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Fourier  transform  ion  cyclotron  resonance  mass  spectrometry  (FTICR/MS) 
equipped  with  various  ionization  techniques  such  as  internal  and  external  electron  impact 
(El)  and  electrospray  ionization  (ESI)  has  been  employed  for  gas  phase  studies  of  two 
different  classes  of  molecules:  organometallic  and  biological. 

One  of  the  organometallic  studies  determined  the  free  energies  of  ionization 
(AGU50o)  for  eight  ruthenium  complexes  with  the  general  formula 
(Ti5-C5Me5)(r|5-C5Me5.nFn)Ru  with  n = 0 - 5.  This  study  was  accomplished  using  an 
internal  El  technique  in  combination  with  electron-transfer  equilibrium  (ETE)  reactions. 
The  additivity  of  fluorine  substituent  effects  in  ruthenocene  ionization  energetics  was  also 
investigated. 

The  other  organometallic  study  investigated  the  reaction  strength  of  Cp2ZrOH+ 

(Cp  = r|5-  cyclopentadienyl)  with  some  selected  alcohols,  esters,  ethers,  amines,  and 
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amides  using  an  external  El  approach.  The  ETE  method  was  also  applied  here  to 
determine  the  free  energies  of  ionization  (AG^  350°)  for  Cp2ZrOR+  complexes,  and  a relative 
AG^  350°  ladder  was  established.  The  external  El  study  produced  improved  results  from  an 
earlier  internal  El  approach,  and  a brief  comparison  of  the  two  experiments  is  given. 

One  of  two  biological  projects  was  to  inspect  the  correlation  between  the  gas 
phase  fragmentation  patterns  of  14  peptides  and  the  locations  within  the  peptides  of  the 
arginine  residue.  Gas-phase  peptide  ions  were  formed  by  ESI  and  infrared  multiphoton 
dissociation  (IRMPD)  produced  by  a C02  laser  was  used  to  fragment  the  ions  in  the  trap. 
The  fragmentation  patterns  obtained  were  compared  with  those  from  sustained  off- 
resonance  irradiation  collisionally  assisted  dissociation  (SORI-CAD). 

The  second  system  of  biological  interest  studied  in  this  work  was  a chemically 
modified  ribonuclease  A (RNase  A).  This  work  showed  that  ESI-FTICR/MS  can  provide 
rapid  scanning  of  sample  purity,  successfully  detect  a single  modification  of  RNase  A with 
N-hydroxysuccimide  acetate  (NHS-acetate),  and  accurately  follow  the  denaturation  and 
chemical  digestion  process  performed  on  RNase  A.  Various  gas-phase  dissociation 
methods  were  attempted  to  dissociate  the  native  and  modified  RNase  A to  pinpoint  the 
exact  location  of  the  modification. 
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CHAPTER  1 
INTRODUCTION 


Fourier  transform  ion  cyclotron  resonance  (FTICR)  mass  spectrometry  (MS)  has 
been  a versatile  and  popular  analytical  tool  since  the  introduction  of  the  concept  of  FTMS 
by  Comisarow  and  Marshall  in  the  early  1970s  .1  The  inventive  combination  of  FT  to  MS 
removed  many  earlier  restrictions  such  as  low  mass  range,  low  mass  resolving  power,  and 
slow  scanning  speed  from  which  MS  had  suffered.  Beside  circumventing  earlier 
limitations,  FTMS  provides  ultrahigh  resolution  mass  measurements,  longer  trapping  time 
for  ions  in  the  cell,  the  ability  to  perform  MS"  experiments,  and  the  flexibility  of 
interfacing  with  many  ionization  techniques  such  as  electrospray  ionization,  laser 
desorption,  matrix  assisted  laser  desorption,  traditional  electron  impact/ionization  and 
chemical  ionization.  For  comparison,  FTICR/MS  can  trap  ions  for  long  periods  of  time 
(up  to  5 x 104  seconds)  compared  to  less  than  a microsecond  ion  lifetime  for  a time-of- 
flight  mass  spectrometer  and  a range  of  10-100  microseconds  for  a quadrupole  mass 
spectrometer.  This  luxury  of  long  ion  trapping  times  conveniently  allows  various  ion 
manipulations  such  as  ion-molecule  reactions  and  multiple  dissociation  experiments  in  ion 
traps.  Although  FTICR/MS  was  originally  used  primarily  for  studies  of  gas-phase 
chemical  reactions,  improvements  made  over  the  past  10  years,  especially  the  successful 
interfacing  of  the  electrospray  ionization  (ESI)  technique  to  MS  in  1984  independently  by 
Aleksandrov  and  co-workers2  and  Yarmashita3  has  revolutionized  the  applications  of  MS. 
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Following  the  successful  application  of  ESI  methods,  FTICR/MS  has  been  extended  into 
the  biological  field,  and  it  has  been  shown  capable  of  studying  large  molecules  of 
biological  interest.  Due  to  its  versatile  applications  in  both  gas  phase  ion/molecule 
reactions  and  biological  areas,  FTICR/MS  is  the  ideal  method  for  investigating 
organometallic  and  biological  systems. 

The  FTICR/MS  characteristic  of  long  ion  trapping  time  at  low  pressure  (~  10'9 
torr)  and  the  ability  to  repetitively  empty  the  cell  of  previous  reactants  and  products  led  to 
many  gas  phase  reaction  studies  including  ion-molecule  reaction  rate  determination4  and 
measuring  physical  quantities  such  as  collision  complex  lifetimes,5  equilibrium  constants,6 
electron  affinities,7  and  ionization  potentials.8  At  low  pressure,  ion  clustering  interferences 
are  eliminated,  and  reactions  are  exclusively  bimolecular. 

The  improvement  of  high  mass  range  detection  since  the  introduction  ofFTMS  has 
enabled  mass  spectrometric  studies  of  ion-molecule  reactions  involving  many 
organometallic  compounds.  For  example,  the  positive  ion/molecule  condensation 
chemistry  of  three  organometallic  compounds  was  reported  by  Parisod  and  Comisarow,9 
the  formation  of  organometallic  ions  in  the  gas  phase  by  reaction  of  Rh+  with  alkanes  was 
studied  by  Byrd  and  Freiser,10  and  metal- switching  reactions  were  discovered  in  which  a 
metal  ion  displaces  another  metal  bound  in  a metallocene  structure  (MCp2).n  In  addition 
to  this  list,  the  thermodynamic  properties  for  many  organometallic  compounds  have  been 
extensively  investigated  and  reported  by  Ryan  and  Richardson.1214 

Thermodynamic  data  are  essential  for  truly  understanding  a chemical  system  and 
its  reactivity.  Although  a large  amount  of  thermochemical  information  has  been  obtained 
for  many  small  organic  and  inorganic  compounds,  this  type  of  information  is  limited  for 
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larger  and  more  complex  organometallic  compounds.  Therefore,  determination  of 
thermodynamic  data  for  organometallic  compounds  is  needed  to  further  their 
understanding. 

The  thermochemistry  of  organometallic  species  has  recently  gained  increasing 
attention  due  to  the  possible  importance  of  these  compounds  in  catalysis.  For  example, 
copper,  zinc,  and  many  transition  metal  organometallic  complexes  play  important  catalytic 
roles  in  biological  systems.15'19  Understanding  these  systems  would  help  to  explain  their 
reaction  pathways  and  catalytic  mechanisms.  Combining  the  electron-transfer  equilibrium 
(ETE)  method12  to  FTICR/MS  can  lead  to  a powerful  technique  for  the  study  of 
organometallic  systems  in  the  gas  phase  and  the  derivation  of  useful  thermodynamic 
information  such  as  free  energies  of  ionization  and  reaction  equilibrium  constants. 

Although  the  development  of  the  FTMS  approach  freed  ICR  mass  spectrometers 
from  many  earlier  limitations,  expanding  the  available  mass  range  from  small  molecules  to 
larger  ones,  FTICR/MS  could  not  be  applied  to  solve  biological  problems  until  the 
introduction  of  new  ionization  techniques  such  as  ESI  and  matrix  assisted  laser  desorption 
ionization  (MALDI).20  One  of  the  reasons  that  MS  was  not  able  to  be  applied  to  biological 
systems  prior  to  the  development  of  ESI  and  MALDI  was  the  problem  of  introducing 
large  ions  into  the  gas  phase  without  destroying  their  structures.  This  problem  was  not 
solved  until  the  successful  interfacing  of  ESI  to  MS.  ESI  and  MALDI  techniques  are 
usually  referred  to  as  “soft”  ionization  techniques,  which  retain  the  structures  of 
biomolecules  and  have  high  sensitivity,  allowing  detection  of  picomole-to-femtomole 
levels  of  material,  thus  enabling  the  direct  analysis  of  biological  fluids  with  a minimum 
amount  of  sample  preparation.  One  difference  between  ESI  and  MALDI  mass  spectra  is 
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that  ESI  generates  multiple  charge  states  while  MALDI  primarily  generates  singly  charged 
molecular  ions.  The  molecular  weights  of  biomolecules  can  be  read  directly  from  a 
MALDI  spectrum,  while  some  simple  calculations  are  required  for  deriving  molecular 
weights  for  ESI  spectrum.  Since  ESI  produces  multiple  charge  states,  and  MS  detects 
m/z  ratios,  most  of  the  mass-to-charge  (m/z)  ratios  of  large  biomolecules  appear  in  the 
100  - 3000  m/z  region  of  the  mass  spectrum.  Generally,  ESI-FTICR/MS  provides  higher 
mass  resolving  power  than  MALDI-FTICR/MS.  Although  the  detailed  mechanism  of  ion 
formation  in  the  ESI  process  still  remains  under  debate  and  speculation,  a substantial 
amount  of  recent  literature  has  shown  that  many  weak  associations  in  solution  do  survive 
the  ESI  process  and  transfer  to  the  gas  phase.21  Some  of  the  weak  associations  which 
have  been  studied  with  ESI-FTICR/MS  include  noncovalent  complexes  such  as  enzyme- 
substrate,  receptor-ligand,  host-guest,  intact  multimeric  proteins,  DNA  duplex  and 
quadruplex  species,  oligonucleotide  complexes  with  drugs  and  proteins,  and  protein-drug 
complexes.22'24 

Due  to  the  availability  of  various  dissociation  techniques  in  FTICR  mass 
spectrometry,  biomolecular  ions  can  be  fragmented  in  the  gas  phase  for  primary  structure 
analysis.  In  addition  to  collision-assisted  dissociation  and  MS"  analysis,25  the  combination 
of  enzymatic  digestion  followed  by  MS  detection  of  fragments  has  been  successfully 
applied  to  sequence  determinations  of  proteins.26  The  result  of  mass  spectrometric 
analysis  of  proteolytic  peptide  mixtures,  especially  after  high  pressure  liquid 
chromatography  (FIPLC)  separation  of  the  peptides,  is  called  a mass-spectrometric  peptide 
mapping.  Given  the  ability  to  perform  direct  analyses  of  multicomponent  proteolytic 
peptide  mixtures  (or  peptide  mapping),  MS  has  already  been  demonstrated  to  be  a 
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powerful  method  for  the  identification  of  covalent  post-translational  modifications27,28  and 
the  unequivocal  characterization  of  multiple  chemical  modification  sites  in  proteins29  by 
mapping  differences  between  the  expected  and  the  observed  primary  structures.  Such 
peptide  mapping  is  very  efficient  compared  with  classical  biochemical  methods,  in  which 
each  collected  chromatographic  fraction  is  sequenced  by  automated  Edman  degradation.30 

One  problem  with  ESI/FTMS  is  that  ionic  buffers  disturb  the  spraying  process  and 
compete  with  analyte  molecules  for  charges.  Some  involatile  substances  can  form  a solid 
cage  which  prevents  escape  of  the  analyte  from  the  evaporating  droplet.  Upper  limits  for 
common  contaminants  in  protein  solutions  are  about  1 mM  for  ionic  salts  and  0.01%  for 
ionic  detergents.30  Therefore,  desalting  the  sample,  either  by  HPLC  or  dialysis,  or  simple 
dilution  is  usually  needed  prior  to  the  sample  introduction. 

Throughout  the  work  presented  in  this  dissertation,  the  versatility  of  FTICR/MS 
has  been  demonstrated  by  applying  the  technique  to  study  two  series  of  organometallic 
complexes  and  two  biological  systems.  Basically,  the  dissertation  can  be  viewed  as  three 
separate  sections:  i)  discussing  the  theory  and  applications  of  a FTICR  mass  spectrometer 
in  chapter  2;  ii)  using  FTICR/MS  and  ion-molecule  reactions  to  investigate  intrinsic 
properties  of  two  sets  of  organometallic  compounds  reported  in  chapters  3 and  4;  and  iii) 
exploring  FTICR/MS  applications  in  the  biological  area  with  two  biologically  related 
projects  presented  in  chapters  5 and  6. 

Theoretical  aspects  of  a mass  spectrometer,  especially  the  basic  principles  of 
FTICR,  including  derivation  of  the  ion  cyclotron  frequency  equation,  ion  motions  in  an 
ICR  cell,  and  ion  excitation  and  detection,  will  be  briefly  discussed  in  chapter  2.  Essential 
components  of  a FTICR  mass  spectrometer  such  as  the  superconducting  magnet,  ultrahigh 
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vacuum  chamber,  ICR  cell,  and  data  processing  station  are  presented.  Ionization 
techniques  such  as  electron  impact  used  for  organometallic  studies  and  electrospray  used 
for  biological  studies  are  also  discussed  in  chapter  2.  The  two  most  commonly  used 
theories  for  explaining  electrospray  ion  formation  are  summarized  in  the  electrospray 
ionization  section. 

Chapters  3 and  4 deal  with  organometallic  studies.  Chapter  3 summarizes  the 
investigative  results  of  free  energies  of  ionization  (AG")  for  eight  ruthenocene 
organometallic  compounds  with  the  general  formula  (Tf-CjMesXrf-CsH^FJRu  (n  = 0 to 
5),  including  two  isomeric  structures,  obtained  using  the  ETE  method  and  FTICR/MS. 
The  questions  of  ligand  and  substituent  additivity  for  these  eight  metallocenes  are 
addressed  with  the  results  obtained. 

Chapter  4 investigates  the  catalytic  reactivity  of  a zirconium  hydroxide  complex, 
Cp2ZrOFT  (Cp  = cyclopentadienyl),  with  various  organic  families  such  as  alcohols,  esters, 
ethers,  amines,  and  amides.  In  the  ICR  cell  the  organic  reagents,  admitted  to  the  high 
vacuum  through  precision  leak  valves,  reacted  with  zirconium  hydroxide  ions,  which  were 
produced  using  electron  impact  (El)  in  an  external  source  and  transferred  into  the  ICR 
cell.  The  reaction  mechanisms  were  checked  with  isotopically  labeled  Cp2Zr018H+ 
complexes  using  a few  alcohol  and  ester  compounds.  The  relative  energies  of  some 
Cp2ZrOR+  products  were  established  by  the  ETE  method  and  FTICR/MS.  The  results 
obtained  from  external  El  and  internal  El  methods  are  briefly  compared,  and  the  logical 
reasons  that  contribute  to  the  different  reaction  results  from  the  two  approaches  are  also 
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Chapters  5 and  6 demonstrate  the  strength  of  FTICR  mass  spectrometry  in 
biological  studies.  One  well  known  advantage  associated  with  mass  spectrometric  analysis 
of  biological  compounds  is  time  saving  partial  protein  sequencing  with  multiple  MS 
experiments.  In  the  work  reported  in  chapter  6,  modifications  of  a protein  with  a known 
sequence  have  been  sequenced  and  identified,  but  it  is  still  impossible  to  sequence  an 
unknown  protein  with  only  MS  even  though  that  is  the  ultimate  goal  for  a structural  mass 
spectrometrist.  However,  if  more  is  known  about  how  locations  of  amino  acid  residues 
within  the  proteins  affect  peptide  fragmentation  in  the  gas  phase,  it  would  make  the 
ultimate  goal  more  reachable.  This  relationship  between  gas  phase  peptide  fragmentation 
patterns  and  the  locations  of  amino  acid  residues  in  a protein  was  investigated,  and  the 
results  obtained  are  reported  in  chapter  5.  Small  peptides  with  3 to  6 amino  acid  residues 
were  dissociated  in  the  ICR  cell  with  the  infrared  multiphoton  dissociation  (ERMPD) 
technique,  a low  energy  dissociation  method,  after  they  were  introduced  into  the  vacuum 
by  electrospray  ionization.  Each  peptide  contained  at  least  one  basic  amino  acid  residue  in 
its  sequence  to  assure  that  each  peptide  was  protonated  so  it  could  be  delivered  into  the 
gas  phase  from  solution.  Other  dissociation  techniques,  especially  the  sustained  off- 
resonance  collisionally  activated  dissociation  (SORI-CAD)  method,31  another  low  energy 
dissociation  method,  were  also  applied  to  these  peptides  for  comparison  purposes. 

Chapter  6 illustrates  the  power  of  FTICR/MS  in  identifying  modification  of  a 
protein  by  following  the  chemical  modification  of  the  enzyme  ribonuclease  A (RNase  A). 

In  this  chapter,  the  advantage  of  FTICR/MS  for  rapid  screening  of  the  purity  of  a sample 
is  demonstrated.  Interpretation  of  the  mass  spectrum  to  identify  post-modified  proteins  is 
elaborated,  and  the  reasons  that  various  gas  phase  dissociation  techniques  failed  to 
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dissociate  the  protein  are  discussed.  The  development  of  a proper  method  for  trying  to 
identify  the  modification  sites  of  the  modified  protein  is  addressed.  The  initial  results 
obtained  from  a different  approach,  which  combines  the  enzymatic  digestion  technique, 
reverse  phase  HPLC  separation,  and  FTICR/MS,  are  reported. 


CHAPTER  2 

THEORY  AND  BACKGROUND 


Mass  spectrometry  (MS)  has  become  one  of  the  most  important  physical  methods 
in  analytical  chemistry  today32  due  to  its  high  sensitivity  and  versatility.  MS  is  one  of  the 
few  methods  that  is  entirely  suitable  for  the  identification  or  quantitative  measurement  of 
trace  amounts  of  chemicals.33  Three  basic  functions  of  a mass  spectrometer  include 
producing  gas-phase  ions  from  sample  molecules  with  various  ionization  techniques  in  the 
source;  separating  and  detecting  gas-phase  ions  according  to  their  mass-to-charge  (m/z) 
ratio  in  the  analyzer;  and  then  processing  the  separated  ions  with  different  data  processing 
techniques.  Some  selected  topics  which  are  relevant  to  this  thesis  will  be  discussed  in  the 
following  sections:  the  El  and  ESI  ionization  techniques,  the  basic  principles  of 
FTICR/MS,  and  four  dissociation  techniques. 

Ionization  Techniques 

The  initial  step  of  the  MS  method  is  to  produce  ions  which  can  be  later  transferred 
into  the  high  vacuum  to  be  analyzed.  Conventionally,  ions  were  produced  using  internal 
ionization  technique  such  as  electron  ionization  (El)  inside  the  mass  spectrometer,  where 
the  ion  analysis  process  also  occurs.  However,  some  recently  developed  external 
ionization  methods,  such  as  electrospray,  produce  ions  at  atmospheric  pressure,  although 
the  ion  analysis  process  still  takes  place  in  a vacuum.  Some  of  the  common  ionization 
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techniques  are  summarized  in  Table  2-1 33  below.  Since  only  internal  and  external  El  and 
electrospray  ionization  techniques  were  used  in  this  thesis,  a more  detailed  discussion  of 
these  techniques  will  be  given  in  the  following  sections. 


Table  2-1. 

Some  common  ionization  methods 


Ionization 

method 

Sample 

preparation  for 
ionization 

Thermal  input 
associated  with 
ionization 

Method 

category 

El 

As  vapor 

Relatively  high 

Thermal 

Cl 

As  vapor 

Relatively  high 

Thermal 

FAB 

Dissolved  in 
matrix 

Virtually  none 

Energetic 

particle 

bombardment 

ICP 

Inductively 

Coupled 

Plasma 

Usually 
dissolved  in 
liquid,  then 
ionized  in  argon 
plasma 

Relatively  high 

Thermal 

GD 

Glow 

Discharge 

Usually  solid, 
then  ionized  in 
high  voltage 
discharge 

Relatively  high 

Thermal 

MALDI 

Mixed  with 
matrix 

Virtually  none 

Energetic 

particle 

bombardment 

ESI 

Dissolved  in 
solvent 

Virtually  none 

Field  desorption 

Electron  Ionization 

Electron  ionization,  also  known  as  electron  impact  (El),  was  the  first  ionization 
method  to  be  used  routinely  in  mass  spectrometry.  It  is  still  one  of  the  most  popular 
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techniques  for  initiating  the  study  of  gas  phase  ion-molecule  reactions.  The  El  source  is  a 
small  enclosure  traversed  by  an  electron  beam  which  originates  from  a heated  filament  and 
is  then  accelerated  through  a potential  of  30-70  eV  into  the  source.  Gas-phase  molecules 
of  interest  in  the  source  then  interact  with  these  electrons.  As  a result,  some  of  the 
molecules  lose  an  electron  to  form  a positively  charged  ion,  and  some  other  molecules 
gain  an  electron  to  form  a negatively  charged  ion.  Depending  on  the  polarity  of  the 
detector,  either  positive  or  negative  ions  can  be  detected.  With  this  method,  many 
molecular  ions  usually  retain  sufficient  internal  energy  after  ionization  to  decompose 
further  to  form  fragment  ions  that  can  be  used  to  characterize  the  structure  of  the  neutral 
molecules.  The  excess  energy  can  also  be  lost  to  the  environment  by  introducing  a time 
delay  which  allows  the  “hot”  ions  to  react  and  transfer  the  excess  internal  energy  to  the 
background  gases  or  a gas  introduced  into  the  cell  for  this  purpose. 

The  basic  ionization  principles  are  essentially  the  same  for  both  internal  and 
external  El  methods  in  FTICR/MS.  The  only  difference  is  the  initial  location  of  ion 
formation.  An  internal  El  technique  produces  ions  within  the  homogeneous  magnetic  field 
while  an  external  El  technique  produces  ions  in  a source  region  which  is  outside  of  the 
fringe  field  of  the  magnet.  The  ions  produced  with  external  El  have  to  be  transferred  into 
the  analyzer  using  some  type  of  ion  guide  or  a series  of  optics.  One  of  the  advantages  of 
using  external  El  is  elimination  of  neutrals  from  the  cell,  which  provides  cleaner  gas  phase 
reactions,  since  only  ions  and  not  their  parent  neutrals  are  transferred  to  the  cell.  Pictorial 
diagrams  of  an  internal  and  external  El  setup  are  shown  in  Figure  2-1. 
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a). 


PV2 


PV1 


Internal  filament 


PV3,4 

PV5,6 


PV1,  PV2  = trapping  plates 

PV3,  4,  5,  6 = excitation/detection  plates 


b).  Magnet 


El  source  Transfer  region 


Analyzer 


PSH  = repelling  plate 

ETP/DETP  = extraction  plates 
SOH  = source  housing 


Figure  2-1. 

Schematic  diagram  of  Bruker’s  a)  Internal  and  b)  External  El  source. 
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Electrospray  Ionization 

The  first  studies  of  the  ESI  phenomenon  date  back  two  and  one-half  centuries,34 
and  extensive  early  studies  of  electrospray  as  an  ionization  method  were  done  by  Malcom 
Dole  and  co-workers  in  1968. 35,36  However,  not  until  the  successful  interfacing  of  ESI 
(see  Figure  2-2)  to  MS  in  19842,3  as  a soft  ionization  technique  to  introduce  structurally 
sensitive  biological  systems  into  the  gas  phase,  did  ESI  become  popular  in  the  MS 
community. 


Heated 

metal 


Figure  2-2. 

Bruker’s  external  ESI  source  with  hexapole  ion  accumulation,  heated  metal 
capillary,  needle  housing,  sample  delivery  tube,  and  syringe. 

In  the  electrospray  ionization  process,  the  dissolved  analyte  molecule  is  protonated 
by  the  solvent  and  then  delivered  to  the  high  voltage  needle  tip  through  a solution  delivery 
tube.  The  syringe  pump  controls  the  flow  rate  to  be  in  the  range  of  0.3-5  pl/min,  and  the 
analyte  is  usually  sprayed  at  a concentration  about  10'6  M in  49:49:2 
methanol/water/acetic  acid  or  acetonitrile/water/acetic  acid.  Addition  of  acetic  acid  to  the 


ESI  solvent  assures  sufficient  protons  to  protonate  basic  analyte  molecules  and  to  provide 
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a stable  spray  current.  Analyte  ions  emerge  as  a fine  spray  of  charged  droplets  from  the 
high  electric  potential  (3-4  kV)  needle  at  atmospheric  pressure.  As  ions  pass  through  the 
capillary  into  the  high  vacuum  chamber  which  hosts  the  mass  analyzer,  the  solvent  in  the 
charged  droplets  evaporates  with  the  aid  of  the  heated  capillary  or  a countercurrent  flow 
of  a drying  gas.  The  capillary  is  located  about  1 cm  away  from  the  needle  and  kept  around 
100  V.  The  production  of  positive  or  negative  ions  depends  on  the  polarities  of  the  needle 
and  capillary.  If  the  needle  is  positively  charged,  positive  ions  are  produced,  and  the 
capillary  is  held  at  lower  positive  voltages  to  allow  ions  into  the  vacuum.  Conversely, 
negative  needle  voltages  produce  negative  ions,  and  a less  negative  potential  on  the 
capillary  transfers  the  negative  ions  into  the  vacuum. 

Although  much  research  has  been  devoted  to  investigation  of  the  mechanism  of 
gas-phase  ion  formation  in  the  ESI  method  (Figure  2-3), 37,38  it  is  still  very  much  under 
active  study. 


Electrospray 

needle 


Figure  2-3. 


ESI  ionization  mechanism 


Currently,  one  of  two  main  theories  used  to  explain  the  mechanism  of  ion 
formation  in  ESI  is  the  Charged  Residue  Model  (CRM),  originated  by  Malcolm  Dole.35  36 
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The  other  theory  is  the  Ion  Evaporation  Model  (IEM),  proposed  by  Iribarne  and 
Thomson.39,40  Dole’s  CRM  is  based  on  experiments  conducted  using  dilute  solutions  of 
polystyrene  molecules  in  a bath  gas  of  nitrogen  at  atmospheric  pressure.  As  the  solvent 
evaporates  from  the  charged  droplets,  the  charge  density  on  the  droplets’  surface  increases 
until  it  reaches  the  Rayleigh  limit.41  At  this  point,  the  Coulombic  repulsive  forces  between 
the  surface  charges  exceed  the  surface  tension,  resulting  in  division  of  the  original  charged 
drops.  If  the  division  process  continues  until  each  droplet  contains  only  a single  molecule 
which  retains  part  of  the  original  charge,  the  mass-to-charge  ratios  of  these  small  droplets 
will  be  detected.  However,  aqueous  electrolyte  solutions  do  not  easily  lose  dissolved  ions 
due  to  an  energy  barrier,  which  is  the  electrostatic  interaction  between  the  ionic  charge 
and  the  dipole  of  the  surrounding  solvent  molecules.  This  energy  barrier  prevents  the  ion 
from  leaving  the  surface.  With  a large  number  of  experimental  results,  Iribarne  and 
Thomson  suggested  that  at  a certain  stage  of  the  evaporation  the  net  charge  on  the  droplet 
with  its  small  radius  was  sufficient  to  produce  an  electric  field  which  was  strong  enough  to 
give  the  ions  an  adequate  potential  energy  to  evaporate  directly  from  the  surface. 

Mass  Analyzers 

After  the  ions  are  produced  in  the  source  region,  they  are  transferred  into  the 
analyzer  region  through  an  ion  guide  or  a series  of  electrical  lenses;  then  a mass  analyzer 
sorts  the  ions  according  to  their  m/z  ratio  and  detects  them.  There  are  different  mass 
analyzers  which  can  be  used  to  achieve  this  purpose,  and  some  of  the  common  ones  are 
presented  in  Table  2-2.  Since  all  data  reported  in  this  thesis  were  collected  using  FTICR 
mass  spectrometers,  the  mass  analyzer  discussion  will  focus  on  only  FTICR/MS. 
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Table  2-2. 

Some  common  mass  analyzers 


Mass 

analyzer 

name 

Quadrup 
ole  MS 

Magnetic 
sector  MS 

Time  of  Flight 
MS  (TOF/MS) 

Ion  trap 
MS 

(IT/MS) 

Ion 

cyclotron 

resonance 

MS 

(ICR/MS) 

Reference 

Dawson, 

197642 

Duckworth, 
Barber,  and 
Venkatasub- 
ramanian, 
198843 

Cotter,  1992;44 
Vestal,  Juhasz, 
and  Martin, 
199545 

March, 

Hughes, 

198946 

Buchanan, 
1 99 147 

Fourier  Transform  Ion  Cyclotron  Resonance  Mass  Spectrometry 

Brief  history.  The  beginning  of  FTICR/MS  can  be  traced  back  to  the  early  1930s 
when  the  cyclotron  accelerator  was  invented  by  Lawrence  and  Livingston.48  Since  then 
FTICR/MS  has  undergone  numerous  stages  of  development.  In  1950,  Sommer  and 
Hippie  developed  the  Omegatron,49  which  was  the  first  cyclotron  instrument  designed  to 
detect  ions  based  on  their  mass-to-charge  (m/z)  ratio.  In  1960,  detection  of  cyclotron 
motion  by  a marginal  oscillator  which  measured  power  absorption  during  excitation  was 
introduced  by  Wobschall. 50,51  However,  the  major  breakthrough  of  FTICR/MS  came  in 
the  mid  1970s  when  Comisarow  and  Marshall  introduced  Fourier  transform  (FT)  analysis 
with  pulsed  and  frequency  sweep  (chirp)  excitation  followed  by  broadband  detection  of  an 
induced  transient  signal.52,53 

The  introduction  of  FT/MS  improved  upon  the  earlier  limitations  of  scanning 
ICR/MS,  which  included  restricted  mass  range,  low  mass  resolution,  and  slow  scanning 
speeds.  Now  FTICR/MS  is  known  for  its  ultrahigh  mass  resolving  power,  accurate  mass 
measurement,  wide  mass  range,  simultaneous  detection  of  all  ions,  ion  remeasurement, 
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long  ion  storage  time,  high  resolution  ion  isolation,  and  multistage  tandem  mass 
spectrometry  (MSn).54 

Although  the  hardware  and  software  of  the  technique  have  changed  over  the  years, 
the  basic  principles  of  FTICR/MS  have  stayed  the  same.  Detailed  descriptions  of  the  basic 
principles  can  be  found  in  many  ICR  papers,55"59  so  only  brief  explanation  will  be  given 
here. 

Derivation  of  the  ICR  equation  FTICR/MS  operates  on  the  basis  of  the 
classical  motion  of  ions  under  the  influence  of  electric  and  magnetic  fields.  When  a 

particle  of  mass  m and  charge  q with  initial  velocity  u is  placed  in  a magnetic  field  B , it 
is  influenced  by  the  Lorentz  force  Fl 

Fl  = q v x B (Equation  2.1) 

and  a centrifugal  force  Fc . 

m u2 

Fc  = (Equation  2.2) 

r 

where  r is  the  radius  of  the  helical  path  that  the  charged  particle  travels  under  the  influence 
of  a magnetic  field. 

As  Figure  2-4  shows,  the  magnitude  of  the  Lorentz  force  must  be  equal  to  the 
magnitude  of  the  centrifugal  force  in  order  to  keep  the  ions  in  the  cyclotron  orbit.  This 
means  that  equation  2. 1 equals  equation  2.2  to  give  equation  2.3. 

m V2 


q t;B  = 


r 


(Equation  2.3) 
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Figure  2-4. 

Forces  acting  on  a charged  particle  in  a magnetic  field. 


Rearranging  equation  2.3  gives  equation  2.4,  which  gives  an  expression  for  the  radius,  r, 
of  the  cyclotron  orbit, 


myi 

"qB 


(Equation  2.4) 


where  u±  is  the  velocity  component  perpendicular  to  the  magnetic  field.  Equation  2.3  can 
also  be  rearranged  to  give  equation  2.5  below. 


ui  qB 


m 


(Equation  2.5) 


VL 

Since  — is  the  expression  for  angular  frequency,  C0c , equation  2.5  can  be  used  to  derive 

r 

the  famous  and  fundamental  ICR  equation. 


(Oc  = 


qB 


m 


(Equation  2.6) 


It  is  also  known  that  cqc  = 2 ttv  , so  ion  cyclotron  angular  frequency  can  also  be  expressed 


in  linear  frequency,  w . 
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qB 

Vc  = (Equation  2.7) 

2 /mi 

According  to  equations  2-6  and  2-7,  the  value  of  an  ion’s  cyclotron  frequency  in  a known 
uniform  magnetic  (no  electric)  field  is  inversely  related  to  its  mass-to-charge  ratio.  Since 
frequency  is  a physical  property  which  can  be  measured  more  accurately  than  any  other, 
FTICR/MS  provides  ultrahigh  mass  measurement  accuracy. 

Ion  motions  in  the  ICR  cell.  If  ions  are  initially  formed  along  the  magnetic  field 
axis,  they  are  confined  radially  by  their  ion  cyclotron  orbital  motion.  However,  without 
any  further  restraint,  ions  would  travel  along  a helical  path  and  be  lost  from  the 
homogenous  region  of  the  magnet  along  the  direction  of  the  magnetic  field  (usually 
defined  as  the  z axis).  In  order  to  prevent  ions  from  being  lost  in  this  direction,  an 
electrostatic  potential  (or  trapping  potential)  is  applied  axially  to  repel  the  ions  back  to  the 
center  of  the  ion  trap,  which  is  located  in  the  most  homogenous  region  of  the  magnetic 
field.  The  typical  values  of  this  electrostatic  potential  are  1-2  volts.  This  introduction  of 
the  trapping  voltage  achieves  the  purpose  of  keeping  ions  from  escaping  along  the  z-axis, 
but  it  also  introduces  unwanted  effects  including:  the  linear  relationship  between  the  ICR 
orbital  frequency  and  the  m/z  ratio  is  no  longer  true,  ICR  line  broadening  and  m/z  shifting 
phenomena  are  introduced,  and  a new  ion  motion,  the  trapping  ion  motion,  arises. 

Squeezing  ions  axially  is  analogous  to  squeezing  a balloon  from  two  ends.  When 
ions  are  forced  back  into  the  center  of  the  trap  from  z-axis,  they  bulge  out  in  the  middle  or 
radially  in  the  trap.  The  outwardly  pushing  force  that  results  from  the  trapping  potential  is 
exactly  opposite  in  direction  to  the  magnetic  component  of  the  Lorentz  force.  The  net 
result  of  the  new  force  is  to  reduce  the  “effective”  magnetic  field  strength.  Since  the  ion 
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cyclotron  frequency,  ryc , is  proportional  to  the  magnetic  field  strength,  introduction  of  the 
trapping  potential  will  lower  the  observed  C0c . 

The  trapping  voltage  in  any  shape  of  an  ion  trap  has  been  well  approximated  by  a 
“quadrupolar”  potential  as  in  equation  2. 8, 60 


where  VT  is  the  electrostatic  trapping  voltage  applied  to  the  trapping  plates,  a and  y are 
numerical  constants  that  depend  on  the  shape  of  the  ion  trap,  and  a is  the  ion  trap 
dimension.  Equation  2.9  gives  the  description  of  the  Lorentz  force  in  the  presence  of  an 

electric  field,  E , 


where  r is  the  vector  position  of  the  ion  at  time  t.  Substituting  the  resulting  electric 
field,  E = -V  V(x,y,z ) , from  equation  2.8  into  2.9,  three  linear  second-order 

differential  equations  are  derived  to  describe  the  motions  of  ions  in  a given  magnetic  and 
electrostatic  trapping  field.  The  equation  involving  only  the  z-coordinate  can  be  solved 
independently  to  give  the  trapping  frequency,  vT,  in  equation  2.10. 


(Equation  2.8) 


Ut  I U K — __  — 

Fl  = m—y  = m— — = qE  + qv  x B 
dt  dt 


d2r  dv 


(Equation  2.9) 


(Equation  2.10) 


The  solutions  for  the  two  remaining  differential  equations  in  x and  y give  the  other  two 
types  of  ion  motions:  the  ion  cyclotron  orbital  motion,  whose  frequency  is  reduced  to  <o 
(equation  2.1 1),  and  a natural  magnetron  motion  with  a frequency  co.  (equation  2. 12). 
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where 


/ 

1 \ 

m 

C 0c  + COc+ 

i- 

{ \ 

YYlcrit  ) 

TYlcn,  = 


qa2B2 
8a  Vt 


(Equation  2.11) 


(Equation  2.11a) 


CO  - 


f 

i > 

m 

(Oc  CQc-\ 

1- 

\ \ 

Mcr.,  ^ 

(Equation  2.12) 


Some  typical  experimental  values  for  the  three  natural  ICR  motions  (ion  cyclotron 
orbital  motion,  trapping  motion,  and  magnetron  motion)  in  a 1 in.  cubic  cell  are  shown  in 
Table  2-3  to  give  one  a good  feel  about  the  relative  frequency  contribution  from  each 
motion.  The  experimental  conditions  used  to  calculate  the  values  in  the  table  are:  a = 1.0 
in.,  a = 1.386  for  a cubic  cell,  B = 3.0  Tesla,  and  VT  = 1.0  V. 


Table  2-3. 

Natural  ion  motional  frequencies  (Hz)  for  singly  charged  ions  trapped  in  a 
quadrupolar  potential  with  a 3 Tesla  magnetic  field  and  1.0  V electrostatic 
trapping  field.61 


Motion 

m=100 

1000 

10,000 

mrri,  =50,492 

v+  =qB/27un 

460,430 

45,837 

4,366 

456 

v.  =aVT/7ia2B 

228 

229 

240 

456 

vT  ~(aqVT/7i2ma2)1/2 

14,496 

4,584 

1,450 
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ICR  cells.  As  discussed  in  the  last  section,  the  ions  have  to  be  confined  in  a finite 
volume  in  order  to  be  excited  and  detected  in  the  homogenous  region  of  the  magnet.  The 
confined  dimension  is  achieved  by  using  an  ion  trap  which  is  known  as  the  ICR  analyzer 
cell  or  just  ICR  cell.  For  FTICR  mass  spectrometrists,  the  ICR  cell  is  usually  described  as 
the  heart  of  the  FTICR  mass  spectrometer  because  it  provides  the  three  most  important 
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functions:  trapping  or  confining  most  of  the  ions  by  a static  homogeneous  magnetic  field 
and  an  electrostatic  trapping  potential,  exciting  and/or  ejecting  their  cyclotron  motion  by 
an  alternating  (r.f.)  excitation  electric  potential,  and  detecting  the  image  currents  induced 
by  the  coherent  ion  cyclotron  motion.  Over  the  years,  many  different  shapes  of  ICR 
cells54,62  have  been  designed  to  meet  the  researcher’s  specific  need  or  to  eliminate  certain 
unwanted  features  of  the  cell.  In  the  next  three  sections,  I will  discuss  three  different  ICR 
cells:  the  standard  cubic  cell,  the  RF-shimmed  infinity  cell , and  the  dual-cell. 

Standard  Cubic  ICR  cell.  The  cubic  cell63  is  the  ancestor  of  ICR  analyzer  cells, 
and  a diagram  of  a typical  cubic  cell  is  shown  in  Figure  2-5.  As  shown  in  the  diagram,  the 
cell  is  made  of  two  trapping  plates,  two  excitation  plates,  and  two  detection  plates.  The 
trapping  plates  are  perpendicular  to  the  direction  of  the  magnetic  field,  and  the  small  hole 
in  the  middle  of  each  trapping  plate  allows  ions  or  laser  beams  to  enter  the  cell. 

Y 


Figure  2-5. 

A standard  cubic  ICR  cell  consists  of  six  plates:  front  and  back  plates  are 
excitation  plates,  top  and  bottom  are  detection  plates,  and  left  and  right  are 
trapping  plates. 
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Rf-shimmed  infinity  cell.  A well  defined  operation  of  the  excitation  event  in  a 
FTICR  pulse  experiment  is  a crucial  prerequisite  for  obtaining  reliable  relative  signal 
intensities,  good  resolution  for  MS/MS  selection,  and  optimum  phase  coherence  of  the 
excited  ion  packets.  Therefore,  finding  a proper  excitation  method  is  on  the  minds  of 
many  research  scientists,  and  one  of  the  most  successful  examples  is  SWIFT,64  which 
assumes  a linear  dependence  of  ICR  signals  with  the  excitation  power.  However,  the 
linearity  requirement  cannot  be  fulfilled  even  in  the  simple  case  of  an  off-resonance  single- 
frequency excitation  pulse.65  The  non-linearity  character  can  lead  to  ejection  of  ions  along 
the  magnetic  field,  and  it  arises  from  the  coupling  of  axial  and  radial  ion  motions  during 
the  excitation  event  in  standard  ICR  cells.  The  excitation  of  the  axial  ion  motion  is 
especially  effective  for  frequencies  equal  to  2oot66  and  coc  + 2oot,67  where  ooc  is  the  cyclotron 
frequency,  and  co,  is  the  trapping  frequency.  Improvements  of  trying  to  reduce  the 
radiofrequency  (rf)  components  of  the  excitation  power,  which  are  parallel  to  the  magnetic 
field  lines,  were  made  separately  by  Marshall’s68  and  Rusell’s69  groups.  In  their  studies, 
the  guard  ring  principle,  introduced  in  the  early  1950s  by  Sommer  et  al.49  for  the 
Omegatron  detector,  was  adapted.  By  properly  adjusting  the  rf  voltages  applied  to  the 
individual  guard  wire  or  guard  ring  the  rf  excitation  field  can  be  shimmed  to  be  nearly 
homogeneous  within  the  ion  trap.  Although  the  ion  excitation  is  somewhat  optimized  in 
this  way,  the  sensitivity  for  detecting  the  image  currents  is  reduced  because  the  guard 
rings  partially  screen  the  detection  electrodes.  Introduction  of  rod  electrodes69  for  the 
detection  restores  the  peak  current  measurement,  but  the  signal  intensity  is  distributed 
over  a series  of  odd  harmonics  of  the  cyclotron  frequency,  kwc  where  k = 1,  3,  5,  . . . 

The  basic  idea  behind  Infinity  cells,  reported  by  P.  Caravatti  and  M.  Allemann  in 
1991, 70  is  to  modify  the  electric  excitation  rf  field  of  a cell  of  finite  dimensions  and  to  give 
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the  ions  an  “illusion”  of  an  infinitely  long  cell.  The  infinity  cell  eliminates  the  parallel 
component  of  the  electric  excitation  rf  field,  which  arises  from  the  finite  dimensions  of  ion 
traps,  and  achieves  the  “illusion”  effect.  As  the  Bruker  Infinity  cell  shows  in  Figure  2-6, 
the  standard  trapping  electrodes  are  replaced  by  a series  of  metallic  segments,  and  these 
segments  exhibit  the  potential  contours  of  an  imaginary  plane  perpendicular  to  the  axis  of 
an  infinity  long  cylindrical  cell.  The  outermost  segments  are  capacitively  coupled  to  the 
excitation  rf  on  the  adjacent  excitation  electrodes,  and  the  connection  of  the  different 
segments  to  an  appropriate  resistive  voltage-divider  chain  provides  the  correct  gradation 
of  the  rf  potential  across  the  trapping  plate. 

RF-shimmed 
trapped  plate 


PVx3 

PVx2 


PVx4  Infinity™  cell 


Figure  2-6, 

Schematic  diagram  of  the  Bruker  Infinity™  cell.  PVx2,  3,  and  4 are  “side  kick” 
plates. 

Dual  cell.  The  dual  cell  (see  Figure  2-7)  contains  two  cubic  cells  that  share  one 
trapping  plate.  These  two  cubic  cells  are  connected  to  each  other  through  a small  orifice 
with  a typical  diameter  of  2 mm.  This  arrangement  provides  the  differential  pumping 
between  the  high  pressure  cell  (source)  and  the  low  pressure  cell  (analyzer)  to  give  a 
pressure  difference  of  103  which  leads  to  higher  resolution.  Ions  are  formed  in  the  source 
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Figure  2-7. 

A dual  cell  showing  two  excitation  plates,  two  detection  plates,  two  trapping 
plates,  and  the  middle  plate  splitting  two  cells.  The  conductance  limit  connects 
two  separated  cells. 


cell  and  then  are  transferred  into  the  analyzer  through  the  orifice  in  the  common  trapping 
plate. 

In  fact,  each  cell  has  all  the  features  of  a single  cell.  Providing  that  the  common 
plate  acts  as  one  of  the  trapping  plates  of  a single  cell,  each  side  of  a dual  cell  can  be  used 
as  an  independent  cell.  If  the  common  plate  is  grounded  during  ion  formation,  a statistical 
equilibrium  of  the  ion  population  between  the  two  cells  is  expected.47,58 

Excitation  and  detection  of  ion  cyclotron  motion  Excitation  is  a process  to 
accelerate  ions  in  the  cell  to  larger  cyclotron  orbits  by  absorbing  energy  from  an  externally 
applied  oscillating  electric  field  (r.f. ) of  the  same  frequency.  When  ions  are  first 
generated,  they  have  different  kinetic  energies.  Therefore,  when  they  enter  the  cell,  ions 
of  different  size  have  different  cyclotron  radii,  and  ions  with  the  same  mass  have  different 
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phases  of  motion.  When  ions  absorb  the  resonant  r.f.  energies,  their  cyclotron  orbits 
increase,  and  they  spiral  outward  (Figure  2-8  b).  If  ions  of  a given  mass  are  irradiated 
long  enough  at  their  resonant  frequencies,  their  cyclotron  orbits  eventually  become  large 
enough  to  strike  the  cell  plates  and  be  neutralized  by  the  plates.  The  process  of  losing  ions 
to  the  cell  plates  is  referred  as  ion  ejection  or  ion  isolation.  During  the  ion 
ejection/isolation  process,  unwanted  ions  are  eliminated  from  the  cell  to  provide  a cleaner 
environment  for  the  ions  of  interest.  When  the  ions  of  interest  are  isolated,  they  are 
resonantly  excited  to  larger  orbits,  and  coherent  ion  motion  is  achieved  at  this  time  (see 
Figure  2.8  c). 

After  a group  of  coherently  excited  ions  has  been  excited  to  larger  orbits  (but  not 
large  enough  to  be  ejected  from  the  cell),  positive  and  negative  charges  will  be  induced  on 
the  detecting  plates  as  shown  in  Figure  2-8  d.  In  this  process,  negative  charges  would  be 
induced  on  the  bottom  plate  from  the  proximity  of  the  positive  ions,  and  positive  charges 
would  be  induced  on  the  top  detection  plate.  However,  when  the  positive  ions  move  away 
from  the  bottom  detection  plate  and  closer  to  the  top  plate,  negative  charges  will  be 
induced  on  the  top  detection  plate  while  positive  charges  are  induced  on  the  bottom 
detection  plate.  Since  the  two  detection  plates  are  connected  in  a closed  circuit,  the  result 
is  an  induced  alternating  current  which  is  usually  referred  as  the  “image  current”. 

The  image  current  is  directly  proportional  to  the  ion  density  inside  the  ICR  cell  as 
given  by  equation  2.13, 


Nq2rB 
sinty? 


(Equation  2.13) 


ma 
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Figure  2-8. 

Schematic  diagram  of  ion  excitation  and  detection,  a)  A positively  charged 
particle  moves  in  a cyclotron  motion  with  a small  radius,  b)  After  the  charged 
particle  absorbs  r.f.  energy  it  increases  its  cyclotron  radius,  c)  “Packed”  ions  move 
together  in  cyclotron  motion  after  coherence  is  achieved,  d)  Image  current  I(t) 
induced  by  the  alternating  positive  and  negative  charges  is  converted  to  voltage 
and  detected  on  the  detection  plates. 
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where  I s(t)  is  image  current,  N is  number  of  ions  in  the  cell,  and  a is  the  dimension  of 


the  cell.  However,  the  direct  linear  relationship  between  the  ion  density  and  the  image 
current  can  be  distorted  if  too  many  ions  are  inside  of  the  cell  due  to  the  space  charge 
effect.  The  image  current  is  converted  to  a voltage,  then  the  voltage  is  amplified  by  a 
preamplifier  and  digitized  by  an  analog-to-digital  converter.  The  voltage  signal  of  each 
ion  mass  is  given  in  the  expression  of  equation  2. 14, 


K(t)  = 
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(Equation  2. 14) 


where  Vs(l)  is  the  voltage  signal,  and  C is  the  circuit  capacitance.  Due  to  ion-molecule 

collisions  inside  the  cell,  the  signal  recorded  decays  with  time.  The  time  domain  signal  is 
Fourier  transformed  to  give  a frequency  spectrum  which  can  be  converted  to  a mass 
spectrum. 

Ions  are  usually  excited  with  one  of  three  common  types  of  excitation:  (1)  impulse 
excitation,  (2)  chirp  excitation,  or  (3)  stored  waveform  inverse  Fourier  transform 
(SWIFT)  excitation. 

Impulse  excitation.  In  the  impulse  excitation  technique,71  a peak  voltage  with  an 
amplitude  of  more  than  100  V and  a pulse  duration  of  about  0.5  ps  is  applied  to  the 
excitation  plates.  During  this  short  period  of  time,  all  ions  in  the  cell  are  accelerated  to 
larger  cyclotron  orbits.  This  method  is  good  for  exciting  ions  with  a broad  mass  range. 

Chirp  excitation.  In  the  chirp  excitation  method,  all  ions  over  the  mass  range  of 
interest  are  resonantly  excited  by  applying  a fast-sweep  rf  excitation  pulse  to  the  excitation 
plates.  The  r.f.  frequency  is  swept  from  the  lowest  to  the  highest  frequencies  to  meet  the 
requirements  of  the  spectrum.  This  is  a simple  implementation  to  excite  a wide  mass 
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range  without  the  need  for  high  r.f.  amplitudes  as  required  in  impulse  excitation.72  The 
problem  associated  with  this  technique  is  that  the  power  spectrum  is  non-uniform  in 
amplitude.  The  non-uniform  power  spectrum  causes  an  uneven  excitation  for  ions  of 
different  m/z  values,  and  the  relatively  broad  frequency  tailing  on  both  sides  of  the  power 
spectrum  limits  the  frequency  selectivity.59 

SWIFT  waveform  excitation.  The  SWIFT  excitation  technique64  eliminates  the 
problems  associated  with  chirp  excitation.  Any  desired  excitation  magnitude  spectral 
profile  can  be  specified  as  a discrete  magnitude  spectrum,  and  then  inversely  Fourier 
transformed  to  produce  the  corresponding  time  domain  excitation  waveform.  The 
resulting  excitation  provides  a optimally  flat  and  selective  excitation. 

ICR  experimental  pulse  sequence.  A typical  pulse  sequence  used  in  an  FTICR 
experiment  is  given  in  Figure  2-9.  The  sequence  starts  with  a quench  pulse,  which 
removes  residual  ions  from  previous  experiments.  Immediately  following  the  quench  pulse 
is  an  ion  formation  pulse  which  allows  the  ions  to  be  produced  and  trapped  in  the 
analyzer.  After  the  selected  ions  are  isolated  during  an  ion  selection  pulse,  a short  delay  is 
added  to  let  the  isolated  ions  relax  back  to  the  center  of  the  cell.  Then,  either  a collision 
gas  can  be  introduced  to  about  10'7  torr  through  the  pulsed  valve  for  a CID  reaction,  or  a 
laser  pulse  can  be  turned  on  for  a photodissociation  experiment.  Next,  a variable  reaction 
delay  can  be  chosen  which  reduces  the  pressure  in  the  analyzer  to  the  original  pressure  of 
1 O'9  torr.  Finally,  the  ions  of  interest  are  excited  and  detected  sequentially  during  the 
excitation  and  detection  pulses. 
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Figure  2-9. 

Representation  of  a typical  ICR  pulse  sequence. 

Hardware.  A FTICR  mass  spectrometer  (see  Figure  2-10)  consists  of  four 
essential  components:  a solenoidal  superconducting  magnet,  an  ultra-low  pressure  vacuum 
system,  an  analyzer  cell,  and  a data  system. 

Magnetic  field.  The  (usually  superconducting)  magnet  used  in  ICR/MS  has  a 
cylindrical  opening  in  the  middle  where  the  ICR  cell  is  located.  Commercially  available 
magnets  which  are  commonly  seen  in  the  research  laboratories  come  in  various  field 
strengths  of  3 to  9.4  Tesla.  Magnets  with  higher  field  strengths  such  as  14  T and  21  T 
have  also  been  constructed  in  a few  national  laboratories.  In  general  the  higher  the 
magnetic  field  strength,  the  better  the  mass  spectrometer  operating  characteristics.  Nine 
advantages  associated  with  a higher  magnetic  field  for  FTICR/MS73  are: 
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Figure  2-10. 

Schematic  diagram  of  a FTICR  mass  spectrometer  with  source,  transfer,  and  cell 
regions. 


The  mass  resolving  power  increases  linearly  with  increasing  the  magnetic  field,  B. 
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(Equation  2.15) 


The  mass  resolving  power  for  two  dimensional  FTICR  spectra  is  improved  by  B2. 

The  speed  of  acquiring  time-domain  data  with  a given  mass  resolving  power  is 
linearly  proportional  to  B.  T is  the  data  acquisition  period. 

m qBT 

- — = 0.132 (Equation  2. 16) 

m 

Ion  cyclotron  kinetic  energy  increases  quadratically  with  increasing  magnetic  field 
strength. 
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(Equation  2. 17) 


The  upper  mass  limit  increases  quadratically  with  increasing  magnetic  field 
strength.  The  variable  a is  a characteristic  dimension  of  the  trap,  and  a is  a 
geometric  factor  determined  by  the  trap  shape. 
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(Equation  2. 18) 


The  ion  trapping  efficiency  is  increased  so  ions  can  be  held  for  longer  period  of 
time  in  the  trap. 


The  ion  density  in  the  cell  is  also  increased  with  increasing  the  magnetic  field 
strength.  /7max  is  the  maximum  ion  density. 

B 2 

"max=  8 \^n  (Equation  2.19) 

Increasing  magnetic  field  linearly  improves  the  quadrupolar  axialization74 
efficiency. 

The  maximum  number  of  ions  defined  by  onset  of  peak  coalescence  is  proportional 
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to  B2.  (Maximum  number  of  ions,  above  which  two  closely  spaced  resonances 

coalesce,  is  much  less  than  the  maximum  number  of  ions  that  may  be  confined  in  a 

Penning  trap.) 

Ultra-high  vacuum.  In  order  to  provide  ultrahigh  mass  resolving  power,  FT  mass 
spectrometers  require  ultra-low  pressures  (<  10"8  torr)  in  the  cell  region.  To  achieve  ultra- 
low  pressures,  differential  pumping  is  most  often  used  since  a substantial  pressure 
difference  has  to  be  achieved  from  atmospheric  pressure,  in  which  the  ions  are  formed 
externally,  to  < 10'8  torr  in  the  cell  where  the  ions  are  excited  and  detected.  In  the  example 
given  in  Figure  2-10,  the  whole  vacuum  chamber  is  divided  into  three  sections  consisting 
of  the  source  region  (« 10'6  torr),  the  transfer  region  (« 10‘7  torr),  and  the  analyzer  region 
(<  10'8  torr).  The  source  region  hosts  various  external  ionization  techniques.  The  transfer 
region  contains  ion  guides  or  electrostatic  lenses  which  are  used  to  guide  externally 
produced  ions  into  the  ICR  cell.  The  analyzer  region  holds  the  most  important  component 
of  the  instrument,  the  ICR  cell.  The  respective  pressure  in  each  region  can  be  achieved 
with  various  types  of  pumps  such  as  oil  diffusion  pumps,  turbo-molecular  pumps,  or 
cryopumps  (liquid  He).  Although  turbo-molecular  pumps  and  cryopumps  are  commonly 
seen  with  the  newer  FTICR  instruments,  oil  diffusion  pumps  are  still  being  used  on  the 
older  instruments.  Details  of  the  vacuum  pumping  setup  for  each  FTMS  instrument  used 
in  this  thesis  are  provided  in  the  proper  experimental  section  of  the  chapter. 

ICR  cell.  See  discussion  of  ICR  cells  given  earlier. 

Data  system.  Different  FTMS  companies  use  different  data  acquisition  and 
analysis  systems  to  operate  their  mass  spectrometers.  For  example,  the  Xmass  software 
with  a Silicon  Graphics  computer  controls  Bruker’s  FTMS  instruments;  an  Odyssey  data 
system  interfaced  with  a Sun  station  operates  the  Finnigan  FTMS  instruments;  and 
Ionspec  uses  PC  computer  systems  for  its  instruments. 
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As  the  block  diagram  of  a FTICR  mass  spectrometer  in  Figure  2-11  illustrates, 
once  the  induced  image  current  is  detected  in  the  ICR  cell,  it  is  converted  to  voltage,  then 
the  voltage  is  amplified  with  a preamplifier,  digitized  using  an  analog-to-digital  converter 
(ADC),  and  stored  in  a computer  to  be  processed  later. 


Figure  2-10. 

A block  diagram  of  a FTICR  mass  spectrometer. 


Dissociation  Techniques 

The  advantage  of  an  ion  trap  MS  is  its  capability  of  holding  the  ions  in  the  trap  for 
as  long  as  5 x 104  seconds  prior  to  detection.75  During  this  long  trapping  time,  beside 
measuring  molecular  weights  of  ions,  various  gas  phase  reactions  can  be  carried  out  with 
these  trapped  ions,  such  as  MS"  experiments,  collisionally  activated  dissociation  (CAD) 
experiments,  and  ion-molecule  reactions.  The  MSn  experiments  combined  with  different 
CAD  methods  lead  to  extensive  examination  of  the  gas-phase  reactions  and  ion 
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dissociation,  which  reveals  useful  structural  information  for  precursor  ions.  In  general,  an 
MS  experiment  is  carried  out  by  applying  a dissociation  method  to  selected  precursor  ions 
to  increase  the  internal  energy  of  the  ions.  Subsequently,  the  precursor  ions  dissociate, 
and  the  fragment  ions  are  excited  and  detected  using  the  standard  procedures.  Ion 
activation  can  be  achieved  by  several  means  summarized  in  Table  2-4.  Only  four 
dissociation  methods  are  discussed  with  details  below. 

Collisionally  Activated  Dissociation 

Collisionally  activated  dissociation,  also  known  as  collisionally  induced 
dissociation  (CID),  was  first  introduced  independently  by  Jennings76  and  McLafferty77 
in  1968.  In  a CAD  experiment,  the  first  step  is  to  isolate  ions  of  interest  (or  precursor 
ions)  by  ejecting  unwanted  ions  using  one  of  the  excitation  methods  discussed  earlier. 

Then  the  precursor  ions  are  resonantly  excited  to  a certain  kinetic  energy  by  a r.f. 
frequency.  Next  the  activated  ions  react  with  a collision  gas  at  a pressure  10'7  torr  or 
higher,  which  is  introduced  into  the  analyzer  cell  through  a pulsed  valve.  In  this  step,  the 
collision  process  converts  the  translational  energy  of  the  precursor  ions  into  internal 
energy,  which  is  followed  by  a unimolecular  decomposition  to  yield  the  products.  After  a 
short  reaction  delay  time,  all  ions  including  the  undissociated  precursors  and  fragment  ions 
are  excited  and  detected.  The  CAD  method  has  worked  very  well  for  small  organic 
molecules,  but  it  does  not  work  well  for  large  biological  systems  because  the  method  can 
not  convert  enough  translational  energy  to  a desirable  amount  of  internal  energy  to  cause 
sufficient  fragmentation  of  the  biological  precursor  ions  due  to  the  very  many  degrees  of 
freedom  that  a large  biomolecule  possesses.88 
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Sustained  Off-Resonance  Irradiation  Collisionally  Activated  Dissociation 

The  problem  encountered  by  the  CAD  method  as  mentioned  above  can  be 
eliminated  with  the  SORI-CAD  technique.  In  the  SORI-CAD  method,  the  precursor  ions 
are  excited  by  an  off-resonance  electric  field  with  the  irradiation  frequency  500-2000  Hz 
away  from  the  cyclotron  frequency  of  the  precursor  ion.  As  a consequence,  the  ions  can 
be  activated  and  deactivated  through  multiple  cycles  with  a period  corresponding  to  the 
difference  between  the  irradiation  frequency  and  the  natural  cyclotron  frequency  of  the 
ions.  This  multiple  activation  and  deactivation  cycle  technique  allows  the  ions  to  be 
irradiated  or  activated  for  longer  periods  of  time  without  being  ejected  from  the  cell,  and 
the  precursor  ions  can  also  undergo  multiple  collisions  with  the  collision  gas  to  accumulate 
enough  internal  energy  to  cause  sufficient  fragmentation.  The  typical  collision  energy  is 
less  than  10  eV,  so  low  energy  fragmentation  pathways  can  be  accessed  with  SORI-CAD. 
SORI-CAD  has  been  used  to  study  fragmentation  pathways  of  nucleotides78  and  the  gas 
phase  structures  of  proteins  such  as  cytochrome  C.79  Disadvantages  of  SORI-CAD 
include  the  formation  of  off-axis  product  ions  and  the  unwanted  ejection  of  product  ions 
that  may  fall  into  the  region  of  the  irradiation  frequency.80 

Photodissociation 

Photodissociation  is  the  generation  of  fragmentation  following  absorption  of 
photons,  and  it  was  initially  applied  to  FTMS  by  Dunbar.81  In  the  PD  method,  ions  absorb 
one  or  more  photons  until  they  gain  sufficient  energy  to  dissociate  into  fragments,  so  the 
primary  requirement  (or  disadvantage)  of  the  method  is  that  the  precursor  ions  must 
contain  a chromophore  which  absorbs  photons  at  the  irradiating  wavelength.  Depending 
on  the  absorption  wavelength  of  the  chromophore,  irradiation  wavelength  ranges  from  193 
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nm  for  UV  photodissociation  to  10.6  pm  for  IR  photodissociation.  Applications  of  PD 
can  be  found  in  many  fields  of  chemistry,  ranging  from  studying  the  energetics  of  gas- 
phase  dissociations  of  small  ions82  to  the  sequencing  of  large  mass  oligonucleotides.86 


Table  2-4. 

Some  common  dissociation  techniques  and  their  primary  references 


Dissociation  technique 

Relevant  reference(s) 

collisionally  activated  dissociation 
(CAD) 

Jennings,  196876;  Haddon  & McLafferty, 
196877;  McLuckey  199283;  Cooks, 
199584;  Cody,  Burnier  & Freiser  198285. 

multiple  excitation  collisional  activation 
(MECA) 

Williams,  Henry,  & McLafferty,  199086; 
Guan,  Hofstadler,  & Laude,  199387. 

sustained  off-resonance  irradiation  for 
collisionally  activated  dissociation 
(SORI-CAD) 

Gauthier,  Trautman,  & Jacobson, 
199131;  Senko,  Speir,  & McLaffety, 
199479. 

very  low  energy  collision-induced 
dissociation  (VLE-CID) 

Boering,  Rolfe,  & Brauman,  199288. 

capillary-skimmer  dissociation  (CS) 

Smith,  199089;  Little,  199490 

surface-induced  dissociation  (SID) 

Mabud,  Dekrey,  & Cooks,  198591 

electron-induced  dissociation  (EID) 

Cody  & Freiser,  198792;  Gord,  199293. 

photodissociation  (PD) 

Chen,  Hays,  & Dunbar,  198481;  Peiris, 
199394. 

blackbody  infrared  radiative  dissociation 
(BIRD) 

Price,  Schnier,  & Williams,  199695 

Capillary-Skimmer  Dissociation 

Capillary-skimmer  dissociation,  also  known  as  nozzle-skimmer  dissociation,  has 
mainly  been  reported  for  studying  large  biomolecules. 90,96  In  this  technique,  the 
dissociation  process  takes  place  at  relatively  high  pressures  in  the  ESI  source  region  where 
the  high  pressure  allows  a large  number  of  low-energy  collisions,  resulting  in  ion 


38 


activation  and  subsequent  dissociation  of  ions.  The  amount  of  dissociation  of  ions 
obtained  from  this  method  is  influenced  by  the  potential  gradient  difference  between  the 
capillary  and  skimmer,  and  the  effective  dissociation  occurs  within  a certain  potential 
difference.  The  product  ions  can  result  from  various  charge  states  of  the  molecular  ions. 
Multiple  dissociation  experiments  cannot  be  carried  out  using  this  technique  because  once 
the  ions  have  left  the  external  ESI  source  and  are  transferred  into  the  analyzer  cell,  they 
can  only  be  manipulated  in  the  cell.  However,  other  dissociation  techniques  can  be  used 
to  gain  further  useful  information  about  the  fragmented  ions  resulting  from  the  CS 
dissociation  process. 


CHAPTER  3 

ADDITIVITY  OF  FLUORINE  SUBSTITUENT  EFFECTS  IN 
RUTHENOCENE  IONIZATION  ENERGETICS 

Introduction 

Ionization  energy  is  an  important  quantity  for  characterizing  the  thermodynamics 
of  oxidation-reduction  reactions  for  molecules.12,14  The  thermochemistry  of  organometallic 
compounds  has  especially  received  increasing  attention  as  researchers  attempt  to  study 
thermodynamic  contributions  to  reactivity  and  catalysis.97  Studies  have  determined  this 
intrinsic  property  for  many  atoms  and  organic  molecules  in  the  gas  phase.98'100  However, 
ionization  energies  have  been  determined  for  far  fewer  organometallic  compounds. 

Traditionally,  thermochemical  quantities  were  measured  with  combustion  reaction 
calorimetry101  and  ionization  energies  were  derived  using  photoelectron  spectroscopy 
(PES). 102,103  The  vertical  ionization  energies  determined  by  PES  can  only  be  an  accurate 
measurement  of  the  adiabatic  gas-phase  ionization  energy  if  the  differences  in  equilibrium 
geometries  of  the  ion  and  neutral  complexes  are  small,12,14  and  any  structural  changes  can 
result  in  a vertical  ionization  potential  which  is  significantly  different  from  the  adiabatic 
ionization  potential  referenced  at  0 K.104  Therefore,  if  one  is  primarily  interested  in 
determining  thermochemical  quantities  near  room  temperature,  the  electron-transfer 
equilibrium  (ETE)  method  coupled  to  statistical  mechanics  has  proven  to  be  an  attractive 
alternative  to  PES.12,14 
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Mautner’s  group  conducted  early  ETE  studies  of  ion  gas-phase  ionization 
thermochemistry  using  pulsed  high-pressure  mass  spectrometry  (PHPMS).105  In  this 
technique,  the  Keq  for  ETE  can  be  extracted  from  a time  plot  of  ion  currents 
corresponding  to  the  reaction  equilibrium,  which  is  established  in  a high  pressure  source 
following  ion  production  by  a pulsed  high  energy  electron  beam.  The  ETE  reaction  can 
also  be  followed  using  ion  trap  mass  spectrometry  such  as  FTICR/MS.  As  described  in 
chapter  2 of  this  dissertation,  FTICR/MS  can  trap  ions  for  many  seconds  after  the  ions  are 
produced  and  transferred  into  the  trap.  The  ions  in  the  cell  are  assumed  to  be  nearly 
thermalized  through  hundreds  of  collisions  within  several  seconds  of  trapping  time  in  a 
10'8  - 10‘5  torr  pressure  range.  This  thermalizing  process  relaxes  the  ions  back  to  center  of 
the  cell,  and  the  thermalized  ions  can  react  more  efficiently  with  the  reaction  reagents 
introduced  through  the  leak  valves. 

The  fluorine-substituted  ruthenium  complexes  (r)5-C5Me5)(r|5-C5F5)Ru  (l)106  and 
(r|5-C5H5)(r)5-C5F5)Ru  (2)107  were  synthesized  and  reported  by  Hughes  and  co-workers, 
and  their  gas-phase  ionization  free  energies  (AG;0)  were  determined  by  Richardson  and  co- 
workers.13 The  gas-phase  ionization  free  energy  study  found  that  fluorine  is  not  an 
unusual  substituent  as  compared  to  the  other  halogens  with  moderately  electron- 
withdrawing  substituents.  The  y parameter  describes  the  effect  of  a given  modified  Cp 
ligand  on  the  ionization  energies  for  the  d6  and  d8  metallocenes.  A y value  of  1.5  was 
obtained  for  the  (r|5-C5F5)  ligand,  indicating  that  five  fluorine  substituents  raised  the 
ionization  free  energy  relative  to  (t^-CjHs).  However,  neither  the  effect  of  an  individual 
fluorine  Cp  substituent  on  oxidation  energetics,  nor  if  the  effects  of  multiple  fluorine 
substituents  are  additive,  are  known  at  this  time.  Since  a fluorine  substituent  possesses  the 
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unusual  combination  of  a strong  electron  withdrawing  effect  and  significant  electron 
donation  to  the  ring,  it  cannot  be  assumed  that  the  effect  is  additive.  Furthermore,  it  is  not 
known  if  ligand  additivity  can  be  assumed  for  r|5-C5F5  and  any  Cp  derivative  L in  (r|5- 
C5F5)(L)M 0/+  metallocene  ionization  energetics,  as  has  proven  to  be  the  case  for  a few 
other  ruthenocenes  LL'  Ru  0/+  (errors  in  predicted  values  <3  kcal/mol).108 

In  the  present  work,  the  ETE  method  with  FTICR/MS  was  chosen  to  investigate 
the  free  energies  of  ionization  (AG")  in  the  gas  phase  for  fluorine-substituted 
ruthenocenes  with  formulas  (Tf-C5Me5)(Tf-C5Me5.nFn)Ru  (n  = 1 - 5)  and  for  (r)5-C5H5)(r|5- 
C5F5)Ru  (2).  The  questions  of  ligand  and  substituent  additivity  for  these  eight 
metallocenes  will  be  addressed  with  the  results  obtained. 

Throughout  this  thesis,  the  ligand  r)5-  cyclopentadienyl  (r|5-C5H5)  and  its 
permethylated  derivative  r|5-  pentamethylcyclopentadienyl  (r|5-C5Me5)  are  abbreviated  as 
Cp  and  Cp*,  respectively,  and  their  structures  are  given  in  Figure  3-1.  Representative 
structures  of  the  fluorine  substituted  ruthenocenes  are  also  shown  in  Figure  3-1. 

Experimental 


Material 

All  eight  fluorinated  ruthenocene  samples  studied  in  this  project,  with  general 
formula  (ti5-C5Me5)(T]5-C5H5  J?JRu  (n  = 1 to  5)  and  (ri5-C5H5)(ri5-C5F5)Ru,  were 
provided  by  Russell  P.  Hughes  of  the  Department  of  Chemistry,  Dartmouth  College.  The 
compounds  (q5-C5Me5)(r|5-C5F5)Ru  (l),106  Cn5-C5H5)(r)5-C5F5)Ru 
(2),107  and  (r)5-C5Me5)(,n5-C5H5.nFJRu  (n  = 1 to  4)13  were  prepared  and  reported  by 
Hughes  and  coworkers.  Nuclear  magnetic  resonance  (NMR)  spectroscopy  and  mass 
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spectrometry  were  used  to  check  the  purity  of  these  samples.  All  the  electron  transfer 
reagents  were  purchased  from  Aldrich,  and  the  liquid  samples  were  degassed  by  several 
freeze-pump-thaw  cycles  prior  to  use. 


(r|  5-C5Me5)(r]  5-C5F5)Ru  (t15-C5H5)(t15-C5F5)Ru 

(1)  (2) 


Figure  3-1. 

A representative  structure  sketch  for  Cp,  Cp*,  (r)5-C5Me5)(r|5-C5F5)Ru  (1), 
and  (r|5-C5H5)(r|5-C5F5)Ru  (2). 


Mass  Spectrometer 

All  experiments  reported  here  were  performed  with  an  FTICR/MS  instrument  that 
included  a 3.0  tesla  horizontal  bore,  room  temperature  Oxford  superconducting  magnet 
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(Oxford  Instruments,  Cryomagnetic  Systems,  Osney  Mead,  Oxford  0X2  ODX,  England), 
a Nicolet  FTMS-1000  vacuum  system,  and  an  IonSpec  personal  computer  based  data 
system.  A schematic  drawing  of  the  instrument  is  given  in  Figure  3-2.  The  magnet  had 
an  inside  diameter  of  6 inches  cylindrical  opening;  the  vacuum  system  had  an  inside 
diameter  of  2.875  inches;  and  the  experiments  in  this  project  were  done  with  a standard  1 
inch  cubic  stainless  steel  ICR  cell. 

The  vacuum  chamber  indicated  in  Figure  3-2  was  evacuated  by  an  Alcatel  300  L/s 
oil  diffusion  pump  which  was  fore-pumped  by  a mechanical  pump,  and  it  maintained  a low 
10'8torr  pressure.  The  ICR  cell  was  located  in  the  center  of  the  homogenous  magnetic 
field,  and  an  internal  filament  assembly  was  located  next  to  the  ICR  cell  and  on  the 
opposite  side  from  the  solids  probe.  The  inlet  system  was  pumped  by  an  Alcatel  150  L/s 
oil  diffusion  pump  which  was  also  backed  by  a mechanical  pump.  The  inlet  system  was 
heated,  with  its  temperature  controlled  at  about  100°C.  Liquid  samples  and  solid  samples 
with  relatively  high  vapor  pressures  were  introduced  into  the  vacuum  through  two  Varian 
precision  leak  valves  contained  in  the  sample  inlet  system.  Solid  samples  with  relatively 
low  vapor  pressures  were  admitted  using  a stainless  steel  solids  insertion  probe.  The 
temperature  of  the  solids  probe  was  controlled  to  provide  the  proper  sample  pressure  for 
the  reactions  of  interest. 

Due  to  the  high  stray  magnetic  fields,  the  nude  ionization  gauge  used  for 
measuring  the  pressure  of  the  vacuum  chamber  was  shielded  with  p-metal  and  placed  in  a 
position  away  from  the  ICR  cell.  This  configuration  resulted  in  pressure  measurements 
which  were  a factor  of  two  lower  than  the  true  pressure  at  the  FTICR  cell  when  the  gate 
valve  of  the  300  L/s  pump  was  in  a l/4th  open  position.109  Since  all  experiments  were 
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conducted  with  the  gate  valve  in  its  1/4  open  position,  all  pressures  measured  and  used  for 
rate  constant  calculations  throughout  the  project  were  multiplied  by  a factor  of  two.  The 
ionization  gauge  was  connected  to  a Granville-Phillips  digital  ionization  gauge  controller. 

The  mass  spectrometer  was  equipped  with  a bake-out  controller  which  allowed  the 
operator  to  heat  the  vacuum  chamber  between  experiments  to  provide  a clean  vacuum  for 
the  next  experiment.  In  the  studies  presented  in  this  chapter,  the  vacuum  chamber  was 
baked-out  for  ca.  8 hours  between  studies  on  different  samples. 


Figure  3-2. 

Schematic  diagram  of  3 telsa  FTICR  instrument  including  standard  cubic  analyzer 
cell,  nude  ionization  gauge,  solids  probe,  and  inlet  system. 
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Results  and  Discussion 


Sample  Preparation 

In  this  project,  the  ETE  method  was  used  for  investigating  the  thermal  equilibrium 
reactions  of  the  eight  ruthenocene  compounds  with  several  organic  reference  compounds. 
(See  Figure  3-3).  The  FTICR/MS  instrument  was  used  to  follow  the  general  reaction 
represented  by  equation  3.1,  where  LL'Ru(g)  is  the  ruthenocene  compound  of  interest,  R is 
the  reference  compound,  and  L and  L'  are  the  ligands  of  the  ruthenocene  compounds. 

LL'Ru(g)  + R+(g)  **  LL'Ru^j  + R(g)  (Equation  3.1) 

The  ruthenocene  derivatives  were  sublimed  from  the  solids  probe  to  give  partial 
pressures  in  the  1 x 10  '7-  1 x 10  "6  torr  range.  For  all  of  the  ruthenocene  derivatives 
additional  heating  of  the  solids  probe  was  required  to  generate  partial  pressures  of  the 
neutrals  in  the  1 x 10  '7  torr  range.  All  the  reference  compounds  were  admitted  to  the 
vacuum  through  a Varian  precision  leak  valve,  and  the  molecular  ions  were  produced  by 
internal  electron  impact.  The  number  of  ions  generated  was  controlled  by  the  electron 
energy  and  the  electron  current  of  the  filament  to  provide  good  ion  signals  with  minimal 
unwanted  fragment  ions.  Electron-transfer  equilibrium  was  established  for  the 
ruthenocene  derivatives  and  the  reference  compounds  when  the  product  ion  intensities 
reached  a constant  ratio.  Even  though  the  equilibrium  was  usually  established  in  less  than 
2 seconds,  the  reaction  was  monitored  for  another  10-13  seconds  to  ensure  true 
equilibrium. 

Data  Treatment 

The  reaction  free  energies  (AGet°)  were  determined  using  equation  3.2.  The 
temperature  of  the  ICR  cell  was  previously  determined  to  be  350  K with  an  Omega  RTD 
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Figure  3-3. 

Free  energy  ladder.  Reference  compounds  are  on  the  left.  AGet°  and  AGU  350° 
values  are  on  the  right.  All  energies  are  in  kcal  mol"1. 
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thin  film  detector,108  and  this  temperature  was  used  in  the  free  energy  calculations. 

AG°ct>350k  = - RT  In  Keq  = - RT  ln{[I(M+)PR]/[I(R+)PM]}  (Equation  3.2) 

The  reactant  ions  and  neutrals  were  allowed  to  equilibrate  at  the  known  pressures  P^^ 
and  PR,  and  the  ion  populations  (I)  in  the  analyzer  cell  were  determined  from  the  mass 
spectrum.  If  the  free  energy  of  ionization  (AG  i 350k)  for  the  reference  compound  (R),  is 
known,  then  the  corresponding  energy  for  the  ruthenocene  compound  can  be  obtained 
from  equation  3.3. 

AG  ^oieCM)  = AG  p50k  (R)  + AG  eti35ok  (Equation  3.3) 

Previously  reported  AH;0  values  and  estimates  of  AS"  for  the  reference  compounds 
have  been  used  to  derive  AG^350°(R)  values  at  350  k.105,110-112  The  pressures  for  neutrals 
were  obtained  directly  from  the  N2-calibrated  ion-gauge  reading  and  were  likely  to  be  low 
for  both  the  metallocenes  and  the  reference  compounds.  However,  the  assumption  was 
made  that  the  reference  compounds  and  the  metallocenes  have  similar  ion  gauge 
sensitivities  and  system  factors  so  that  the  equilibrium  constants,  which  depend  on  the 
pressure  ratios,  could  be  obtained  with  acceptable  accuracy.  Previous  pressure  calibration 
studies  in  our  laboratory108  showed  that  the  assumption  is  approximately  true  for  the 
reference  compounds  and  for  non-fluorinated  metallocenes  if  the  pumping  speed  in  the 
main  chamber  is  reduced  substantially.  Thus,  the  relative  ionization  free  energies  are 
estimated  to  be  accurate  within  ±0.5  kcal  mol'1.  Due  to  pressure  errors,  it  is  estimated 
that  the  absolute  errors  are  ± 1.5  kcal  mol'1.  Note  that  an  unusually  large  pressure  ratio 
error  of  5/1  would  lead  to  an  error  of  ±1. 1 kcal  mol'1  in  AGct°  at  350  K. 
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Free  Gas-Phase  Ionization  Energies  ( AG/)  of  the  Eight  Fluorine  Substituted 
Ruthenocenes 

Typical  results  for  the  ETE  experiments  with  the  ruthenocenes  and  the  organic 
reference  compounds  are  shown  in  Figure  3-4.  These  data  were  collected  with  a total 
pressure  of  2.7  x 10'7  torr  and  a Cp*RuC6HF4  to  Cp2Ru  ratio  of  1.3:1.  The  AGU50°  value 
of  Cp*RuC6HF4  is  165.2  kcal  mol'1.  The  measured  AGt350°  values  for  the  ETE  reactions 
of  all  eight  fluorine  substituted  ruthenocenes  with  various  reference  compounds  are  also 

T 

shown  in  Figure  3-3  as  a free  energy  ladder.  The  derived  values  of  AG^  350°  are  displayed 
next  to  the  compounds  in  the  ladder.  The  value  obtained  in  this  work  for  1 agrees  within 
± 0.1  kcal/mol  with  that  reported  earlier.13  Typical  calculations  of  Keq,  AGet°  and  AG^,/ 
values  are  illustrated  with  the  reaction  of  Cp*RuC6HF4  and  Cp2Ru  in  equation  3.4.  The 
reference  reagents  were  first  ionized  with  internal  El  to  give  Cp2Ru+.  Then  the  Cp2Ru+ 
ions  reacted  with  Cp*RuC6HF4  sublimed  from  the  solids  probe.  The  reaction  equilibrium 
was  achieved  within  8 seconds  of  reaction  delay  as  shown  in  Figure  3-4. 

Cp*RuC6HF4  + Cp2Ru+  * Cp*RuC6HF4+  + Cp2Ru  (Equation  3.4) 

A plot  of  AG^  350°  vs  the  number  of  fluorines,  n,  in  (r)5-C5Me5)(r|5-C5Me5.nFn)Ru  (n 
= 0 - 5)  is  shown  in  Figure  3-5.  The  data  points  form  a straight  line  with  a slope  of  3.7  ± 
0.5  kcal  mol'1  n'1  (95%  confidence  interval)  and  r2  = 0.96.  The  results  show  that  fluorine 
substituent  effects  are  essentially  additive,  with  each  additional  fluorine  raising  the 
ionization  energy  by  ~ 3.7  kcal  mol'1. 

y Parameter  of  Ligand  Substituent  Effect 

The  parameter  y,  developed  earlier  by  Ryan  et.,108  describes  the  effect  of  a given 
modified  Cp  ligand  on  the  ionization  energetics  for  the  d6  (and  d8)  metallocenes.  This 
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Figure  3-4. 

A typical  ETE  experiment  of  the  reaction  of  Cp*RuC6HF4  and  Cp2Ru  complexes. 


Figure  3-5. 

Plot  of  AG^  350°  values  vs  number  of  fluorines.  The  solid  line  is  the  linear 
regression,  and  the  dotted  lines  are  95%  confidence  intervals. 
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parameter  is  assigned  to  the  entire  ligand  rather  than  any  particular  substituent  since  the 
sum  of  individual  substituent  parameters  does  not  necessarily  describe  the  effect  of  multi- 
substituted  Cp  ligands  accurately.  However,  an  additive  ligand  effect  is  assumed  so  that 
the  effect  of  two  different  Cp  ligands  is  assumed  to  be  described  by  the  sum  of  the 
individual  y parameters. 


Table  3-1. 

Free  Energies  of  Ionization  and  Derived  y(L)  Values 


Ligand 

AGL  350o  (kcal  mol1) 

Y(L) 

fi5-C5H5 

152.3 

0 (defined) 

ti5-C5H4F 

154.6 

0.23 

r)5-l,3-C5H3F2 

157.6 

0.46 

ri5-l,2-C5H3F2 

158.6 

0.54 

t!5-1,2,3-C5H2F3 

161.5 

0.76 

fi5- 1 ,2,4-C5H2F3 

162.5 

0.84 

t!5-C5HF4 

165.2 

1.0 

fi5-C5F5 

170.8 

1.5  (1.7  ± 0.2)a 

See  text  for  discussion 


The  general  equation  formulated  to  derive  y parameters  is  given  in  equation  3.5 
AG;°(L  L'M)  = a ( yL  + y L.)  + AGj°(Cp2M)  (Equation  3 . 5) 

with  L and  L'  as  Cp  derivatives  and  a as  a “sensitivity  constant”  which  depends  on  the 
metal  M.  AGi°(Cp2M)  is  ionization  free  energy  of  Cp2M.  For  ruthenium  metal,  the 
previously  determined  constant  a is  13.0108  and  AG°(Cp2]VI)  is  164.6  kcal  mol'1.  With  the 
experimentally  determined  AGj°  values  given  in  Table  3-1,  y parameters  for  Cp 
derivatives  in  ruthenocenes  can  be  determined  using  equation  3.6  by  defining  y(Cp*)  = -1 
and  y(Cp)  = 0.  The  defining  values  of  y(Cp*)  = -1  and  y(Cp)  = 0 imply  that  a negative 
value  for  yL  means  that  a ligand  L stabilizes  the  ionized  form  of  the  metallocene  relative  to 
the  neutral,  lowering  the  ionization  energy. 
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AGi°(L  L'M)  = 13.0  (yL  + y L.)  + AGi°(Cp2Ru)  (Equation  3.6) 

The  derived  y values  for  the  fluorinated  Cp  derivatives  are  tabulated  in  Table  3-1 
along  with  that  of  Cp  for  comparison.  As  Table  3-1  shows,  the  overall  trend  is  clearly 
additive  with  respect  to  F substituent  effects.  However,  some  small  but  reproducible 
deviations  from  additivity  are  observed.  For  example,  when  n = 2 or  3,  the  two  structural 
isomers  appear  to  have  about  1 kcal  mol'1  difference  in  AG;0  values,  but  no  particular  trend 
is  noticed  with  respect  to  the  isomer  structures.  Further  examination  of  Table  3-1  reveals 
a larger  jump  in  the  ionization  free  energy  values  for  n = 4 to  5 (5.7  kcal  mol'1)  than  the 
rest.  But  the  jump  is  still  within  the  95%  confidence  interval  as  Figure  3-5  shows,  and  it 
may  represent  a “perfluoro”  effect.103  Since  that  effect  is  small  enough,  it  may  be 
accounted  for  by  experimental  error. 

With  the  AGj0  value  obtained  from  this  study  and  the  previous  definition  of  y(Cp) 
being  0,  y^-CsFj)  is  determined  to  be  1.9  for  r|5-C5F5  with  Cp  as  the  other  ligand  using 
equation  3.6.  However,  Y(r|5-C5F5)  is  1.5  if  the  other  ligand  is  Cp*.  Since  C5Me5  is  a 
more  electron-rich  ligand  than  C5H5  it  is  expected  to  reduce  the  ionization  energy-raising 
effect  of  r|5-C5F5  more  effectively  than  Cp  itself.  However,  the  number  of  systematically 
studied  compounds  with  respect  to  ligand  additivity  is  small,  and  the  poorest  additivity 
observed  previously  (L  = indenyl)108  also  has  a large  range  for  its  value  of  y (-  0.4  ± 0.2). 
With  this  small  ligand  additivity  effect  (the  range  of  y(C5F5)  = 1.7  ± 0.2  based  on  the  two 
r)5-C5F5  metallocenes  studied  to  date),  it  can  be  concluded  that  r)5-C5F5  is  not  unusual  with 


respect  to  additivity. 
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Conclusion 

The  free  energies  of  ionization  for  seven  fluorine  substituted  ruthenocenes,  (r)5- 
C5Me5)(r|5-C5Me5.nFn)Ru  (n  = 1 - 5),  and  (r)5-C5H5)(r|5-C5F5)Ru  (2)  were  investigated  and 
determined  with  the  ETE  method,  and  the  values  ranged  between  152.3  and  170.8 
kcal/mole.  r|5-C5Fl5  has  the  lowest  gas-phase  free  ionization  energy  while  r]5-C5F5  has  the 
highest  gas-phase  free  ionization  energy.  The  results  of  this  study  show  that  each 
additional  fluorine  raises  the  ionization  free  energy  by  3.7  ± 0.5  kcal  mol'1,  and  the  energy 
difference  between  the  two  possible  isomers  for  n = 2 and  n = 3 was  found  to  be  less  than 
or  equal  to  1 kcal  mol'1.  This  study  suggests  that  the  intrinsic  effect  of  each  fluorine 
substituent  on  ionization  energy  is  additive  without  any  unusual  effect  observed.  The 
large  jump  from  n = 4 to  5 (5.7  kcal  mol'1)  may  represent  a “perfluoro”  effect,  but  the 
effect  is  small  and  may  be  accounted  for  by  experimental  error. 


CHAPTER  4 

INTRINSIC  REACTIVITY  OF  METAL-HYDROXIDE  COMPLEXES: 
GAS-PHASE  REACTIONS  OF  [Cp2ZrOH]+  WITH  ALCOHOLS, 

AMINES,  ETHERS,  ESTERS,  AND  AMIDES 

Introduction 

Many  transition-metal  compounds  have  been  known  to  catalyze  biopolymer 
cleavage  in  solution.  This  catalytic  property  could  provide  insights  into  methods  for  gas- 
phase  fragmentation  of  peptide  and  ester  bonds  since  the  rapid  cleavage  of  esters  has 
usually  been  a prerequisite  for  catalysis  of  peptide  bonds  and  phosphate  ester  hydrolysis  in 
solution-phase  studies.  Metal-bound  hydroxides  are  found  in  a variety  of  synthetic  and 
biological  catalysis  that  are  active  in  hydrolysis  of  carboxylic  acid  and  phosphoric  acid 
derivatives.  Many  transition-metal  complexes  such  as  those  containing  palladium(II),15 
zinc(II),16  cobalt(III),17  and  copper(II)19  have  been  found  to  be  important  catalysis  in 
biological  systems  for  the  cleavage  of  RNA,  DNA,  amides,  and  carboxylic  and  phosphate 
esters.  In  general,  the  pH  dependence  of  catalysis  can  be  used  to  show  that  a hydroxide 
complex  (LxM-OH"+)  or  its  kinetic  equivalent  is  responsible  for  the  observed  reactivity. 

In  reactions  with  activated  esters  (those  with  good  leaving  groups),  metal 
hydroxides  act  as  typical  O-donor  nucleophiles,  and  the  reactivity  is  generally  predictable 
from  the  pKa  of  the  metal-bound  water  .113  With  unactivated  substrates  (which  are  better 
models  for  biological  substrates),  rates  of  catalyzed  hydrolysis  can  be  unusually  high,  and 
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the  concept  of  bifunctional  catalysis  can  be  invoked  in  which  a nucleophile  (hydroxide) 
and  an  electrophile  (the  metal  center)  are  in  close  proximity  and  can  act  in  tandem  to 
increase  the  catalytic  reaction  rate  (e.g.,  via  a zinc  bound  intermediate  shown  here). 
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The  zirconium- hydroxide  complex  (Cp2ZrOH+)  is  a 14-electron  zirconium(IV)- 
hydroxo  complex  which  could  serve  as  a model  for  bifunctional  catalysis  with  the  potential 
for  gas-phase  studies  to  reduce  the  complexity  of  the  system.  Its  chemical  features  are 
similar  to  known  solution  catalysts  (i.e.,  the  metal  center  is  not  coordinatively  saturated 
and  is  strongly  electrophilic),  and  this  suggests  that  it  would  be  quite  reactive. 

Coordinative  unsaturation  can  lead  to  rapid  binding  of  Lewis  basic  substrates  since  ligand 
displacement  need  not  occur.  Therefore,  the  zirconium-hydroxide  complex  has  been 
chosen  as  a model  to  study  the  reactivities  of  metal-bound  complexes.  In  this  chapter,  the 
results  obtained  from  the  investigation  of  the  intrinsic  reactivity  of  the  Cp2ZrOH+  complex 
with  various  organic  groups  such  as  alcohols,  esters,  ethers,  amines,  and  amides  in  the  gas 
phase  with  FTICR/MS  are  reported  and  discussed. 


Experimental 


Materials 

The  complex  Cp2Zr(CH3)2  was  obtained  from  Professor  Richardson’s  laboratory  in 
the  Inorganic  Division  of  the  Chemistry  Department  at  the  University  of  Florida.  Since 
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the  Cp2Zr(CH3)2  sample  was  sensitive  to  moisture  in  the  air,  it  was  kept  in  a dry  box  until 
it  was  ready  to  be  transferred  into  the  ionization  chamber.  All  organic  reagents  except  4- 
nitrophenyl  acetate,  which  was  recrystallized  by  Celeste  Regino  in  the  Inorganic  Chemistry 
Division  at  the  University  of  Florida,  were  commercially  purchased  and  used  without 
further  purification. 

Mass  Spectrometer 

All  experiments  in  this  project  were  performed  with  a FTICR/MS  instrument 
which  consists  of  a 3.0  Tesla  Oxford  superconducting  magnet  (Oxford  Instruments, 
Cryomagnetic  Systems,  Osney  Mead,  Oxford  0X2  ODX,  England),  a Bruker  external 
source  vacuum  chamber  (Bruker  Analytical  Systems,  Inc.,  Billerica,  MA),  and  an  Aspect 
3000  data  system.  The  schematic  diagram  of  the  instrument  is  given  in  Figure  4-1.  The 
vacuum  chamber  of  the  instrument  is  divided  into  three  separately  pumped  regions: 
source,  ion  transfer,  and  analyzer.  In  the  source  region,  various  types  of  ionization 
techniques  such  as  external  electron  ionization  and  chemical  ionization  take  place,  and  this 
region  is  evacuated  to  about  10"6  torr  with  a 330  1/s  Balzer  turbomolecular  pump  which  is 
forepumped  by  a mechanical  pump.  The  ion  transfer  region  contains  the  electrostatic 
lenses  which  guide  ions  into  the  analyzing  chamber,  and  it  is  pumped  by  a 
70  1/s  Balzer  turbomolecular  pump  (Balzers,  Hudson,  NH).  Finally,  the  analyzer  region, 
which  contains  the  rf-shimmed  Infinity™  cell,  is  evacuated  to  a base  pressure  in  the  low 
10'9  torr  range  using  a 330  1/s  Balzer  molecularturbo  pump.  The  turbomolecular  pumps  of 
the  ion  transfer  region  and  the  analyzer  region  share  one  mechanical  pump.  The  source 
chamber  is  connected  to  the  main  vacuum  with  a small  orifice  which  can  be  blocked  by  a 
manual  gate  valve.  All  sample  inlet  valves  and  the  solids  probe  loading  port  are  pumped 
by  an  Edwards  mechanical  pump. 
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Figure  4-1. 

Schematic  diagram  of  the  FTICR  mass  spectrometer  including  3 tesla 

superconducting  magnet  and  a Bruker  external  source  vacuum  chamber. 

The  mass  spectrometer  is  controlled  by  an  Aspect  3000  minicomputer  with  24  bit 
word  length.  The  computer  is  equipped  with  192  K of  memory,  a FT  co-processor  for 
fast  data  treatment,  and  a fast  ADC.  The  spectrometer  operates  in  both  high  resolution 
and  broad  band  FT  modes. 

Sample  Handling 

The  starting  complex,  Cp2Zr(CH3)2,  which  was  air  sensitive,  was  kept  in  a dry  box 
until  it  was  ready  to  be  transferred  to  the  ionization  chamber.  The  first  step  of  the  sample 
transfer  process  was  to  load  a small  quantity  of  the  Cp2Zr(CH3)2  powder  into  a melting 
point  capillary,  which  was  then  sealed  with  grease  to  prevent  moisture  contamination 
before  it  was  taken  out  of  the  dry  box.  Prior  to  loading  the  sample-  containing  capillary 
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on  a solids  probe  tip,  the  capillary  was  cut  down  to  the  proper  length  so  it  would  fit 
securely  in  the  probe  tip.  The  air  exposure  time  of  the  sample  was  kept  to  a minimum  to 
reduce  the  amount  of  moisture  contamination.  After  the  sample  tube  was  secured  on  the 
solids  probe  tip,  the  probe  was  inserted  into  the  first  region  of  the  probe  loading  port, 
which  was  then  rapidly  evacuated  by  a mechanical  pump  to  about  3.0  x 10'2  torr.  Next, 
the  valve  separating  the  solids  probe  port  from  the  source  chamber  was  opened  to  allow 
the  probe  to  be  inserted  into  the  source  chamber.  This  allowed  the  sample  to  be  located 
right  next  to  the  filament  assembly  of  the  external  El  source.  This  setup  enabled  the 
sample  vapor  to  be  easily  ionized  by  the  electrons  produced  using  the  external  El 
ionization  technique. 

The  temperature  of  the  solids  probe  was  maintained  around  40-60  °C  to  assure  that 
there  was  sufficient  sample  vapor  in  the  source  region  to  be  ionized  so  that  optimal  signals 
for  the  ions  of  interest  could  be  obtained.  The  ions  were  produced  with  a proper  electron 
energy  to  maximize  the  signals  of  the  parent  ions,  Cp2Zr(OH)+,  and  to  minimize  the 
formation  of  the  fragment  ion,  Cp2Zr+.  The  typical  value  of  electron  energy  used 
throughout  this  study  was  30  eV.  H20,  D20,  and  the  organic  reagents  were  introduced 
into  the  high  vacuum  through  two  Varian  precision  leak  valves,  and  the  amount  of  vapor 
leaked  into  the  vacuum  was  adjusted  to  provide  adequate  reactivity  of  the  ions  of  interest. 

Results  and  Discussion 


Preparation  of  Cp2ZrOH+ 

Equations  4. 1 and  4.2  together  illustrate  the  formation  of  Cp2Zr(OH)+  from  the 
neutral  Cp2Zr(CH3)2.  As  equation  4. 1 shows,  the  first  step  was  to  produce  Cp2Zr(CH3)+ 
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in  the  source  region  using  the  external  electron  impact/ionization  technique.  Then  the 
Cp2Zr(CH3)+  ion  was  transferred  into  the  analyzer  cell  from  the  ionization  source  through 
a series  of  electrostatic  lenses.  The  second  step  of  the  ionization  process  was  to  introduce 
H20  or  D20  into  the  high  vacuum  through  a Varian  precision  leak  valve,  so  they  could 
react  with  the  Cp2Zr(CH3)+  ions  in  the  cell  to  produce  zirconium  hydroxide  ions, 
Cp2Zr(OH)+,  as  illustrated  in  equation  4.2.  The  amount  of  H20  vapor  leaked  into  the  cell 
was  adjusted  to  maximize  the  conversion  of  Cp2Zr(CH3)+  to  Cp2Zr(OH)+. 

Cp2Zr(CH3)2  + e"  -*  Cp2Zr(CH3)+ + CH3  • + 2 e'  (Equation  4. 1) 
Cp2Zr(CH3)+  + H20  ->  Cp2Zr(OH)+  + CH4  (Equation  4.2) 

After  the  Cp2Zr(OH)+  ions  were  prepared,  they  were  isolated  using  a combination 
of  chirp  and  single  shot  excitation.  They  then  reacted  with  organic  reagents  which  were 
introduced  into  the  analyzer  cell.  The  amount  of  organic  reagent  leaked  into  the  vacuum 
was  adjusted  to  provide  an  optimal  yield  of  the  product  ions  of  interest.  A variable 
reaction  delay  applied  to  each  experiment  controlled  the  amount  of  the  reaction. 

All  five  isotopes  of  zirconium  were  observable  in  the  spectra,  with  m/z  235  and 
237  being  the  most  abundant  isotopes  for  Cp2ZrCH3+  and  Cp2ZrOH+,  respectively.  The 
Cp2Zr  metallocene  moiety  is  assumed  to  remain  intact  for  all  reactions.  The  gas-phase 
reactions  of  the  isolated  complex  ions,  Cp2ZrOH+,  with  selected  alcohols,  ethers,  esters, 
amines,  and  amides  will  be  summarized  and  discussed  in  the  following  sections. 

Reactions  of  Alcohols  with  Cp2ZrOH+ 

The  results  of  the  reactions  with  alcohols  are  given  in  Table  4-1.  It  can  be  seen 
from  the  Table  that  alkoxide  complexes,  Cp2ZrOR+,  are  formed  exclusively  for  all  the 
alcohol  compounds  after  elimination  of  an  H20  molecule  during  the  reaction,  as  equation 
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4.3  illustrates.  The  product,  Cp2ZrOR+,  of  the  reaction  was  observed  immediately 
following  the  introduction  of  the  alcohol  vapor  into  the  ICR  cell.  Further  experiments 
have  shown  that  the  first  product,  Cp2ZrOR+,  reacted  further  with  excess  ROH  in  the  cell 
to  produce  the  second  product  complex,  [Cp2ZrOR+  ROH]+,  at  longer  reaction  delays. 

Cp2Zr^OlT  + Rqil  ->  Cp2ZrOR+  + H20  (Equation  4.3) 

Cp2Zr-OR+  + ROH  -►  [Cp2ZrOR  + ROH]+  (Equation  4.4) 


Table  4-1. 

Results  obtained  for  the  reactions  of  Cp2ZrOH+  with  selected  alcohols.  Only  the 
main  product  ion  of  the  reaction  is  listed  here. 


Name  of  the  compound 

Product 

m/z 

Neutral  lost 

CH3OH  (methanol) 

Cp2ZrOCH3+ 

251 

H20 

CH3CH2OH  (ethanol) 

Cp2ZrOCH2CH3+ 

265 

h2o 

C6H5OH  (phenol) 

Cp2ZrOC6H5+ 

313 

h20 

(CH3)3COH  (t-butanol) 

Cp2ZrOC(CH3)3+ 

293 

h2o 

(CH3)2CHOH  (i-propanol) 

Cp2ZrOCH(CH3)2+ 

279 

h2o 

CH3(CH,)3OH  (butanol) 

Cp,ZrO(CH7)3CH3+ 

293 

h7o 

Typical  spectra  obtained  for  the  reactions  of  Cp2ZrOH+  with  alcohols  are  shown  in 
Figures  4-2  and  4-3.  Figure  4-2  is  a spectrum  of  the  isolated  Cp2ZrOH+  ions,  and  Figure 
4-3  shows  the  formation  of  the  gas-phase  product,  Cp2ZrOCH2CH3+,  with  a 2.5  s reaction 
delay.  The  set  of  peaks  at  m/z  237. 1 in  Figure  4-2  represents  the  parent  ion  Cp2ZrOH+, 
and  the  set  of  signals  at  m/z  265.2  are  due  to  the  product  ion,  Cp2ZrOCH2CH3+. 
Comparison  of  Figures  4-2  and  4-3  clearly  demonstrates  that  the  Cp2ZrOH+  ions  were 
immediately  converted  to  the  main  gas-phase  product  Cp2ZrOCH2CH3+.  Peaks  due  to  the 
minor  product,  [Cp2ZrOC2H5  + HOC2H5]+,  at  m/z  3 1 1 .2  are  also  seen  in  the  spectrum 


given  in  Figure  4-3. 
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Figure  4-2. 

Mass  spectrum  of  isolated  Cp2ZrOH+  ions  with  no  reaction  delay.  The  spectrum  is 
the  accumulation  of  16  scans  at  a pressure  of  ca.  2 x 10'8  torr. 


Figure  4-3. 

Mass  spectrum  obtained  for  the  reaction  of  Cp2ZrOH+  with  ethanol  following  a 
2.5s  reaction  delay.  The  signal  at  m/z  237. 1 is  the  most  abundant  isotope  of 
Cp2ZrOH+,  and  m/z  265.2  is  the  most  abundant  isotope  of  the  product 
Cp2ZrOCH2CH3+.  The  unlabeled  peak  at  m/z  3 1 1 .2  is  the  minor  product 
[Cp2ZrOCH2CH3  + CH3CH2OH]+. 
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Reactions  of  selected  alcohols  with  Cp2Zr018H+.  Isotopic  labeling  studies  for 
selected  organic  reagents  with  Cp2Zr018H+  were  performed  to  investigate  the  reaction 
mechanism,  and  the  results  obtained  are  summarized  in  Table  4-2.  The  reaction 
mechanism  of  the  alcohols  with  the  zirconium  hydroxide  complex  can  be  illustrated  using 
the  well  accepted  four-center/four-electron  complex  explanation  as  illustrated  in  Figure  4- 
4.  For  example,  the  reactions  of  Cp2Zr018H+  with  ethanol  and  phenol  suggest  that  the 
product  Cp2ZrOR+  was  formed  by  replacing  the  -018H  group  of  Cp2Zr018H+  with  the  -OR 
group  of  the  alcohols  after  the  elimination  of  a H2018  molecule  during  the  reaction. 


Table  4-2. 

Results  for  the  reactions  of  Cp2Zr018H+  with  selected  reaction  reagents 


Name  of  the  compound 

Product 

Neutral  lost 

CH3CH2OH  (ethanol) 

Cp2ZrOC2H5+ 

H2018 

C6H5OH  (phenol) 

Cp2ZrOC6H5+ 

h2o18 

(CH3)2C02  (methyl  acetate) 

(Cp2Zr018H  + C3H602)  + 

CH3C02C2H5  (ethyl  acetate) 

(Cp2Zr018H  + C4H802)  + 

CH3C02H  (acetic  acid) 

Cp2ZrO,CCH3+ 

X 

s) 

o 

00 

Determination  of  relative  Zr-O  bond  energies  (AAG°Zj..0)  for  selected 
Cp2ZrOR+  complexes.  Table  4-3  summarizes  the  free  energy  of  the  reactions  (AG°eq) 
and  the  relative  Zr-0  bond  energies  (AAG0^)  of  selected  Cp2ZrOR+  complexes. 
Although  the  absolute  Zr-0  bond  energies  of  Cp2ZrOR+  complexes  cannot  be  determined 
from  this  study,  the  relative  Zr-0  bond  energies  of  Cp2ZrOR+  complexes  can  be  estimated 
as  illustrated  below.  The  kinetic  reactions  can  be  described  using  equation  4.5, 

[Cp2ZrOR]+  + R'OH  - [Cp2ZrOR']+  + ROH  AG°eq  (Equation  4.5) 
where  R and  R'  are  H,  methyl,  ethyl,  i-propyl,  s-butyl,  or  t-butyl.  Equation  4.5  can  be 
broken  down  into  three  separate  equations: 
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Figure  4-4. 

Reaction  scheme  of  an  alcohol  compound  with  Cp2Zr018H+.  T.S.  is  the  reaction 
transition  state. 


Table  4-3. 

Free  energies  (AG°eq)  determined  for  selected  reactions  described  in  equation  4.5, 
and  the  relative  bond  strength  of  Zr-0  (AAG°Zr.0)  determined  for  each  reaction 
pair. 


Reaction  pair 

AG7M(kcal/mol) 

AAGVn  ± 0.7(kcal/mol) 

Cp2ZrOH7CH3OH 

2.6-3. 2 

12.22 

Cp2ZrOCH37C2H5OH 

0.7-0. 8 

0.49 

Cp2ZrOCH37i-propanol 

1.2-1. 6 

1.88 

Cp2ZrOCH37t-butanol 

1.5-1. 6 

2.17 

Cp2ZrOC2H57i-propanol 

0.7 -0.8 

1.35 

Cp2ZrOC2H57s-butanol 

1.3-1. 4 

2.64 

Cp2ZrOC7H,7t-butanol 

1.0-1. 1 

1.93 
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[Cp2ZrOR]+  * [Cp2Zr  ]+  + -OR 

AG  Zr_OR 

(Equation  4.6) 

[Cp2Zr-]+  + -OR'  * [Cp2ZrOR']+ 

■AG 

(Equation  4.7) 

•OR  + R'OH  * OR'  + ROH 

AG°™ 

(Equation  4.8) 

From  the  analog  of  Hess’s  law,  the  addition  of  equations  4.6,  4.7,  and  4.8  gives  back 
equation  4.5,  so  the  values  of  AG°eqcan  be  expressed  as: 

AG°eq  = AGVor-  AGVor,  + AG0™  (Equation  4.9) 

Let  AG  zt-or_AG  7j--or'  = AAG  Zr.0,  and  rearrange  equation  4.9  to  obtain  the  expression 
given  in  equation  4. 10,  which  can  be  used  to  calculate  the  relative  Zr-0  bond  energy: 
AAGV0  = AG°eq  - AG°™  (Equation  4. 10) 

In  equation  4. 10,  AG°™  is  approximately  equal  to  AH0™  by  assuming  that  under  the 
current  experimental  conditions  the  contribution  of  TAS°™is  insignificant  compared  to 
the  values  of  AH°™.  The  TAS°™  term  is  ignored  for  all  calculations  in  this  chapter. 
AH°™  can  be  determined  using  Hess’s  law  with  the  standard  enthalpy  information  given 
in  the  CRC  handbook,114  and  AG°eq  is  determined  using  equation  4.11, 

AG°eq  = - RT  In  Keq  (Equation  4. 1 1) 

where  R in  equation  4.11  is  the  gas  constant,  and  T is  the  temperature  of  the  ions  in  the 
ion  trap.  For  external  ionization  sources,  ions  are  initially  formed  outside  of  the  cell  and 
have  to  be  transferred  into  the  cell  through  a series  of  electrostatic  lenses.  Ions  in  the  cell 
are  assumed  to  be  at  room  temperature,  298.15  K,  and  this  value  was  used  for  all  AG°eq 
calculations  in  this  chapter.  The  reaction  equilibrium  constant,  Kcq,  can  be  determined 
using  the  expression  given  in  equation  4.12, 


[Cpi'ZrOR  ] PRQHaROH 
r 

[Cp2ZrOR\ 


PR'OHaR'OH 


(Equation  4. 12) 
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where  I is  the  normalized  intensity  of  the  ion,  p is  the  pressure  of  the  neutral  gas  in  torr, 
and  a is  the  pressure  correction  factor  to  correct  the  sensitivity  difference  that  the  ion 
gauges  might  have  toward  different  neutrals  in  the  cell.115 

Typical  equilibrium  data  obtained  for  the  alkoxide  reactions  described  in  equation 
4.5  are  illustrated  using  the  kinetic  reaction  of  CpjZrOQH/  with  ethanol  given  in  Figure 
4-5.  Figure  4-5  is  a kinetic  plot  of  the  natural  logarithm  of  the  normalized  ion  intensity, 
ln(I),  versus  the  reaction  delay  in  seconds.  It  can  be  seen  from  Figure  4-5  that  as  the 
starting  ion  Cp2ZrOC4H9+  disappears  with  the  reaction  delay,  the  product  ion 
Cp2ZrOC2H5+  increases  in  concentration.  After  a 3 second  reaction  delay  is  applied,  both 
products  Cp2ZrOC4H9+  and  Cp2ZrOC2H5+  lose  their  ion  intensities  due  to  the  formation  of 
three  new  adduct  products,  [Cp2ZrOC2H5+C4H9OH]+,  [Cp2Zr0C2H5+C4H90H]+,  and 
[Cp2ZrOC4H9+C4H9OH]+.  As  the  kinetic  plot  in  Figure  4-5  shows,  a reaction  equilibrium 
between  the  CpjZrOCJH/  complex  and  the  Cp2ZrOC2H5+  complex  was  established  after  5 
seconds  of  reaction  delay.  This  reaction  equilibrium  was  used  to  determine  the  Keq  value 
of  the  reaction. 

The  formation  of  three  adduct  products  for  the  reaction  of  t-butoxide  complex 
with  ethanol  is  consistent  with  the  other  alcohols  studied  earlier.  At  longer  reaction 
delays,  some  of  the  product  ions,  Cp2ZrOR+,  were  converted  to  similar  adduct  products 
observed  for  the  reactions  of  the  t-butoxide  complex  and  ethanol.  These  adduct 
complexes  have  one  common  feature  - four  electron  donating  ligands  around  the 
zirconium  metal  center.  This  suggests  that  the  most  stable  structure  of  the  zirconium 
complexes  seen  in  this  study  consists  of  four  electron  donating  ligands  including  two  Cp 
complexes,  one  -OR,  and  one  ROH.  The  zirconium  complexes  with  four  electron 
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donating  substituents  are  also  seen  for  other  classes  of  organic  compounds  such  as  ethers, 
amines,  and  amides  (these  results  are  discussed  later). 


Figure  4-5. 

Plot  of  ln(I)  vs.  reaction  delay  for  the  reaction:  [Cp2ZrOC4H9]+  + C2H5OH  ^ 
[Cp2ZrOC2H5]+  + C^OH.  The  reactant  ions  and  the  product  ions  are  labeled  on 
the  plot. 

The  relative  Zr-0  bond  energies  of  the  Cp2ZrOR+  complexes  are  represented  as  an 
energy  ladder  given  in  Figure  4-6.  The  vertical  axis  of  the  diagram  is  AAG0^  in 
kcal/mol,  and  the  right  column  of  the  diagram  represents  the  -R  substituent  of  the 
Cp2ZrOR+  complex.  All  reactions  were  studied  with  either  the  methoxide  complex  or  the 
ethoxide  complex.  It  can  be  seen  from  Figure  4-6  that  the  AAG°Zr.0  value  obtained  for 
the  reaction  of  Cp2ZrOH+  and  Cp2ZrOCH3+  is  10  times  larger  than  the  rest  of  the  AAG°Zr_0 
values.  In  fact,  the  AAG°Zr.0  values  between  the  alkoxide  complexes  are  relatively  small, 
with  a typical  value  less  than  4 kcal/mol,  as  shown  in  Figure  4-7. 

During  the  earlier  derivation  of  AG°eq  in  this  chapter,  AG°Zr.OR-  AG^or-  was 
defined  as  AAG°Zr.0.  This  definition  relates  AAG°Zr.0  of  the  reactions  to  the  relative  Zr- 
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O bond  energies  of  the  Cp2ZrOR+  complexes.  A small  AAG°Zr_0  value  means  that  the  Zr- 
O bond  strength  is  similar  between  the  two  alkoxide  complexes  while  a large  AAG^q 
value  suggests  a large  difference  in  the  Zr-0  bond  strength  of  the  complexes.  Hence,  the 
small  AAG0^  values  obtained  for  the  alkoxide  complexes  as  shown  in  Figure  4-6 
suggest  that  the  Zr-0  bond  strength  between  the  alkoxide  complexes  is  similar.  Figures 
4-6  and  4-7  show  that  the  AAG0Zr.o  value  increases  as  the  -R  group  gets  larger.  This 
trend  implies  that  an  alkoxide  complex  with  a larger  -R  group  can  form  a stronger  Zr-0 
bond  than  a complex  with  a smaller  -R  substituent.  The  formation  of  the  stronger  Zr-0 
bond  for  the  complex  with  larger  -R  substituent  is  due  to  the  polarizablility  property  of  the 
-R  group.  The  larger  the  -R  group,  the  easier  it  can  be  polarized,  and  the  better  an 
electron  donating  group  it  becomes.  Therefore,  a stronger  Zr-0  bond  can  be  formed  for 
alkoxide  complexes  with  larger  -R  substituents.  Since  Zr  is  an  electron-insufficient  metal 
center  in  the  complexes,  it  binds  more  strongly  with  a larger  R group,  such  as  a butyl 
group,  compared  to  a smaller  methyl  group.  From  the  experimental  results,  the  order  of 
the  Zr-0  bond  strength  for  the  alkoxide  complexes  can  be  obtained  as  follows:  Cp2ZrOH+ 
< Cp2ZrOCH3+  < Cp2ZrOC2H5+  < Cp2ZrOCH(CH3)2+  < Cp2ZrOCHCH2(CH3)2+  < 
Cp2ZrOC(CH3)3+ 

The  AAG°Zr.0  value  determined  for  the  t-butoxide  complex  from  the  reaction  of 
the  methoxide  complex  is  0.25  kcal/mol  lower  than  the  AAG0^  value  obtained  from  the 
reaction  of  ethoxide.  The  AAG0^  value  between  the  methoxide  complex  and  the 
isopropoxide  complex  is  0.04  kcal/mol  lower  than  the  AAG0^  value  estimated  from  the 
reaction  of  ethoxide  complex  and  the  isopropanol  compound.  These  differences  are 


probably  due  to  experimental  errors. 
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Figure  4-6. 

The  relative  Zr-0  bond  strength  (AAG°Zr_0)  in  ± 0.7  kcal/mol  for  selected 
Cp2ZrOR+  complexes  where  R = H,  methyl,  ethyl,  i-propyl,  s-butyl,  and  t-butyl. 
The  diagram  is  not  drawn  to  scale. 
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Figure  4-7. 

The  plot  of  the  AAG°Zr.0  values  (±  0.7  kcal/mol)  determined  from  the  reaction 
pairs  studied. 
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Reactions  of  Cp2ZrOH+  with  Other  Organic  Families:  Esters,  Ethers,  Amines, 
and  Amides 

Reactions  of  Cp2ZrOH+  with  four  other  classes  of  organic  compounds  were 
investigated  using  the  same  method  as  for  the  alcohol  studies.  The  results  obtained  for  the 
reactions  of  these  four  classes  of  organic  compounds  are  summarized  in  Table  4-4.  As 
the  Table  shows,  except  for  acetic  acid,  the  adduct  products,  [Cp2ZrOH  + RC02R', 

ROR',  RNHR',  or  RCONHR']+,  were  exclusively  produced  for  the  reactions  of 
Cp2ZrOH+  with  esters,  ethers,  amines,  or  amides,  respectively.  The  general  reaction  can 
be  illustrated  by  equation  4.13. 

Cp2ZrOH+  + reaction  reagents  (esters,  ethers,  amines,  amides) 

— ► adduct  product  (Equation  4. 13) 

Isotopic  studies  of  the  reactions  of  Cp2Zr018H+  with  acetic  acid,  methyl  acetate, 
and  ethyl  acetate  were  carried  out  under  the  same  experimental  conditions  as  used  for  the 
unlabeled  studies.  The  results  obtained  are  given  in  Table  4-2.  The  isotopic  labeling 
studies  generated  similar  results  to  the  unlabeled  studies.  For  methyl  acetate  and  ethyl 
acetate,  the  adducts  were  still  observed  as  the  main  gas  phase  products,  and  the  reaction 
mechanism  is  given  in  equation  4. 14. 

Cp2Zr018H+  + CH3COOCH3  or  C2H5COOC2H5  — ► adduct  product  (Equation  4. 14) 

The  reaction  of  the  isotopic  labeling  study  for  acetic  acid  is  illustrated  in  equation 
4.15.  The  experimental  results  showed  that  Cp2Zr02CCH3+  was  the  primary  gas  phase 
product.  The  isotopic  studies  confirmed  that  Cp2Zr02CCH3+  was  formed  by  replacing  the 
-018H  substituent  of  Cp2Zr018H+  with  the  -02C2H3  part  of  the  acetic  acid  and  the 
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Table  4-4. 

Data  collected  for  the  reactions  of  Cp2ZrOH+  with  various  esters,  ethers,  amines, 
and  amides. 


Name  of  the  compound 

Product 

m/z 

*CH3C02H  (acetic  acid) 

((Cp2Zr02CCH3) 

279 

CH3C02CH3  (methyl  acetate) 

(Cp2ZrOH  + CH3C02CH3)+ 

311 

CH3C02CH2CH3  (ethyl  ester) 

(Cp2ZrOH  + CH3C  02CH2CH3)+ 

325 

CH3(CH2)2C02CH3  (methyl  butyrate) 

(Cp2ZrOH  + CH3(CH2)2C02CH3)+ 

339 

C6H5C02CH2CH3  (ethyl  benzoate) 

(Cp2ZrOH  + C6H5C02CH2CH3)+ 

387 

C1(F)2C202CH3 

(Cp2ZrOH  + C1(F)2C202CH3)+ 

381 

(F)3C202CH2CH3  (ethyl  trifluoroacetate) 

(Cp7ZrOH  + (F)3C707CH7CH4)+ 

379 

C8H704N  (4-nitrophenyl  acetate) 

(Cp2ZrOH  + C8H704N)+ 

418 

(CH3CH2)20  (diethylether) 

(Cp2ZrOH  + (CH3CH2)20)+ 

311 

((CH3)2CH)20  (diisopropylether) 

(Cp2ZrOH  + ((CH3)2CH)20)+ 

339 

(CH3CH2CH2)2NH  (dipropylamine) 

(Cp2ZrOH  + (CH3CH2CH2)2NH)+ 

338 

(CH3CH2)3NFf  (triethylamine) 

(Cp2ZrOH  + (CH3CH2)3NH)+ 

338 

HCON(CH3)2  (N,N-dimethyl  formamide) 

(Cp2ZrOH  + HCON(CH3)2)+ 

310 

CH3CON(CH3)2  (N,N-dimethylacetamide) 

(Cp7ZrOH  + CH3CON(CH3)7)+ 

324 

* indicate  different  product  is  observed. 


elimination  of  one  H2018  molecule  during  the  reaction.  The  reaction  mechanism  of  the 
acetic  acid  reaction  is  similar  to  the  one  proposed  earlier  for  the  alcohols,  but  it  is  different 
from  the  rest  of  esters  studied  in  this  chapter. 

Cp2Zr018H+  + CH3COOH  — ► Cp2Zr02C2H3+  + H2018  (Equation  4. 1 5) 

One  possible  explanation  for  the  discrepancy  in  results  occurring  within  the  ester 
family  (with  general  formula  RC02R')  relates  to  the  -R'  group  of  the  compound.  When 
comparing  the  formula  of  acetic  acid  (CH3C02H)  to  methyl  acetate  (CH3C02CH3)  and 
ethyl  acetate  (C2H5C02C2H5),  it  can  be  seen  that  the  H atom  is  the  -R'  group  for  acetic 
acid  which  is  same  as  the  -R'  of  the  alcohols  (ROH),  but  the  -R'  groups  of  methyl  acetate 
and  ethyl  acetate  are  -CH3  and  -C2H5,  respectively.  Since  the  ion/molecule  reactions  of 
Cp2ZrOH+  with  alcohols  and  acetic  acid  were  observed  in  this  study,  and  the  similar 
reactions  were  not  observed  for  other  esters,  these  results  suggest  that  the  H atom  as  -R' 
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of  the  organic  compound  plays  an  important  role  for  the  formation  of  the  final  products. 
When  the  -R'  is  an  alkyl  group  other  than  the  H atom,  only  adduct  products  are  formed 
for  the  reactions  of  the  Cp2ZrOH+  complex.  When  the  -R'  of  the  organic  compound  is  a 
H atom,  an  ion/molecule  reaction  occurs. 

In  the  reactions  for  the  rest  of  the  ethers,  esters,  amines  and  amides,  only  adduct 
products  were  observed.  The  results  obtained  are  summarized  in  Table  4-4.  These 
reactions  will  be  examined  thermodynamically  and  kinetically  to  investigate  the  factors 
that  prevent  the  formation  of  the  expected  gas-phase  products  which  are  seen  for  acetic 
acid  and  the  alcohols. 

The  studies  of  the  alcohol  family  have  shown  that  the  Zr-0  bond  of  a Cp2Zr-OR+ 
complex  is  relatively  stronger  than  the  Zr-0  bond  of  the  Cp2Zr-OH+  complex.  If  Zr-02R+ 
was  formed  as  the  final  product  for  the  reaction  of  an  ester  compound,  its  Zr-0  bond 
strength  would  be  stronger  than  one  single  Zr-0  bond  and  weaker  than  two  single  Zr-0 
bonds.  Therefore,  the  formation  of  Zr-02R+  for  esters  would  be  more  favorable  than  the 
formation  of  Zr-OR+  for  alcohols,  and  this  can  be  demonstrated  numerically.  With  the 
values  AAG0^  shown  in  Figure  4-7  and  the  derivation  steps  illustrated  in  equations  4.6 
to  4. 10,  the  AAG°Zr.0  values  for  the  reactions  of  Cp2ZrOH+  with  methyl  acetate  and  ethyl 
acetate  are  estimated  to  be  greater  than  |-12.22|  kcal/mol.  With  the  values  of  the  standard 
enthalpy  of  formation  found  in  CRC  handbook,  AG0^  for  the  methyl  acetate  and  ethyl 
acetate  reactions  is  determined  as  - 19.36  and  - 5.9  kcal/mol,  respectively.  After  AAG°Zr^, 
and  AG°  are  determined,  AG°eq  of  the  reactions  for  methyl  acetate  and  ethyl  acetate  can 
be  estimated  using  equation  4.9.  The  derivations  show  that  AG°eq  is  <-31.58  kcal/mol  for 
the  methyl  acetate  reaction  and  < -25.26  kcal/mol  for  ethyl  acetate.  The  negative  AG°eq 
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values  indicate  that  both  reactions  of  Cp2ZrOH+  with  methyl  and  ethyl  acetates  are 
thermodynamically  favorable  for  producing  the  products  Cp2Zr02CCH3+  and 
Cp2Zr02CC2H5+,  respectively.  However,  there  is  no  evidence  that  such  products  were 
experimentally  observed. 

The  other  probable  explanation  for  why  the  expected  products  Cp2Zr02CCH3+  and 
Cp2Zr02CC2H5+  were  not  observed  for  methyl  and  ethyl  acetates  may  be  due  to  kinetics. 
For  example,  the  adduct  products  are  probably  the  reaction  intermediate,  and  the 
activation  energy  barriers  for  the  proposed  four-center/four-electron  transition  states  (see 
Figure  4-8)  are  too  high  to  be  overcome  under  the  experimental  conditions;  therefore,  the 
favorable  transition  structures  which  could  lead  to  the  expected  gas-phase  product  are 
prohibited. 

Figure  4-9  illustrates  two  possible  activation  energy  barriers  that  could  prevent  the 
conversion  of  the  adducts  to  the  final  expected  gas-phase  products.  One  possible 
activation  energy  structure  is  that  the  transition  energy  is  simply  too  high  for  the 
intermediate  structure  to  get  over  the  energy  barrier  (See  Figure  4-9,  a)  and  lead  to  the 
expected  gas  phase  products.  The  second  possibility  is  that  there  is  an  energy  minimum  at 
the  transition  state  as  shown  in  Figure  4-9.  b.  In  this  case,  the  transition  states  (see 
Figure  4-8)  do  not  have  enough  energy  to  roll  themselves  out  of  this  energy  valley  and  to 
convert  the  transition  state  structures  to  the  expected  final  products. 

The  same  argument  holds  for  the  reactions  of  ethers,  amines,  and  amides.  Even 
though  the  reactions  are  thermodynamically  favorable  for  complete  reactions,  the  kinetic 
data  are  not  favorable  for  complete  reactions  to  form  the  expected  final  products. 
Therefore,  the  expected  gas-phase  products  Cp2Zr02CR+,  Cp2ZrOR+,  Cp2ZrNHR+,  and 


72 


ether 


amine 


T.S. 

ester 


i i 
H 

T.S. 

amide 


Figure  4-8. 

Proposed  transition  states  (T.S.)  for  the  reactions  of  Cp2ZrOH+  with  ester,  ether, 
amine  or  amide  reagents. 


Figure  4-9. 

The  possible  structures  of  the  transition  state  which  could  inhabit  the  formation  of 
the  expected  final  products,  Cp2Zr02CR+,  Cp2ZrOR+,  Cp2ZrNHR+,  and 
Cp2ZrOCNHR+  for  the  reactions  of  Cp2ZrOH+  with  esters,  ethers,  amines,  and 
amides. 
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Cp2ZrOCNHR+  for  the  reactions  of  esters  (except  acetic  acid),  ethers,  amines,  and  amides 
never  have  the  opportunity  to  be  formed. 

Comments  on  the  Results  Obtained  from  Both  External  and  Internal  El  Methods 

The  results  obtained  with  the  external  El  method  have  shown  fewer  reaction 
channels  than  the  results  obtained  from  the  internal  El  method.  For  example,  the  external 
El  approach  eliminated  reaction  channels  involving  Cp2Zr(CH3)2  neutral  which  were 
observed  using  the  internal  El  approach.  Since  the  external  El  technique  generates  ions  in 
the  ionization  source,  which  is  located  some  distance  from  the  analyzer  cell,  and  since  only 
ions  can  be  transferred  into  the  analyzer  cell  through  a series  of  electrostatic  lenses,  the 
unionized  neutrals  in  the  source  region  are  pumped  away  before  they  can  enter  the  ICR 
cell.  However,  for  the  internal  El  technique,  the  filament  assembly  was  located  adjacent 
to  the  ICR  cell,  and  therefore  any  unionized  Cp2Zr(CH3)2  complex  could  easily  enter  the 
reaction  cell  to  complicate  the  desirable  reaction  pathways.  Comparisons  of  the  results 
obtained  from  both  the  external  and  internal  El  techniques  demonstrate  that  the  external 
El  approach  provides  a cleaner  environment  for  the  gas  phase  ion/molecule  reactions  than 
the  internal  El  approach. 


Conclusions 

This  study  has  shown  that  under  the  current  experimental  conditions  a Cp2ZrOR+ 
complex  was  produced  when  the  zirconium-hydroxide  complex  (Cp2ZrOH+ ) reacted  with 
selected  alcohols  (ROH).  The  results  obtained  from  the  isotopic  studies  have  confirmed 
that  the  Cp2ZrOR+  product  was  formed  by  cleaving  the  metal  oxygen  bond  of  Cp2ZrOH+. 
The  -OH  substituent  of  Cp2ZrOH+  was  replaced  with  the  -OR  group  of  the  alcohol. 
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However,  only  adduct  products  were  produced  for  the  reactions  of  Cp2ZrOR+  with  other 
classes  of  organic  compounds  such  as  esters,  ethers,  amines,  and  amides.  Unfavorably 
high  energy  barriers  may  be  responsible  for  the  lack  of  the  expected  final  gas-phase 
products.  Although  the  O18  labeling  studies  show  that  adduct  products  were  observed  for 
the  methyl  acetate  and  ethyl  acetate  reactions,  Cp2Zr02CCH3+  was  seen  as  the  primary 
product  for  acetic  acid  after  a H2018  molecule  was  eliminated  from  the  reaction  complex. 

In  this  study,  the  free  energies  (AG°eq)  for  the  reactions  of  several  Cp2ZrOR+ 
complexes  were  determined  and  then  were  employed  to  estimate  the  relative  Zr-0  bond 
energies  (AAG°Zr.0)  of  the  Cp2ZrOR+  complexes.  The  AAG0^  values  determined 
suggest  that  as  the  -R  group  of  the  Cp2ZrOR+  complex  gets  larger,  the  Zr-0  bond 
becomes  stronger.  The  results  also  show  that  the  Cp2ZrOH+  complex  has  the  weakest  Zr- 
O bond  among  all  Cp2ZrOR+  complexes  studied.  This  work  has  demonstrated  that  the 
zirconium-hydroxide  complex  can  affect  hydrolysis  of  alcohols  in  the  gas  phase,  but  the 
hydrolysis  reactions  were  not  observed  for  ethers,  esters,  amines,  and  amides. 


CHAPTER  5 

INFRARED  MULTIPHOTON  DISSOCIATION  (IRMPD)  STUDIES  OF  SMALL 
PEPTIDES  IN  THE  GAS  PHASE  WITH  ESI-FTICR/MS 

Introduction 

Peptide  sequencing  with  mass  spectrometry  (MS)  has  gained  much  attention  in 
recent  years  since  proper  ionization  techniques  to  introduce  charged  peptides  into  the  gas 
phase  have  been  available.  A major  effort  of  mass  spectrometric  peptide  studies  has 
focused  on  understanding  the  dissociation  patterns  of  peptides.116,117  A better 
understanding  of  peptide  dissociation  patterns  will  obviously  lead  to  better 
characterization  of  larger  biological  systems,  so  successful  peptide  sequencing  using  MS 
can  provide  some  insights  into  the  gas  phase  protein  sequencing.  For  example,  when 
peptide  ions  of  differing  amino  acid  sequences  have  identical  mass,  they  can  not  be  easily 
differentiated.  One  means  of  achieving  this  is  to  perform  MS"  experiments  to  sequence 
their  primary  structures.  If  a general  correlation  exists  between  peptide  fragmentation 
patterns  and  the  location  of  a specific  amino  acid  residue  within  a peptide  sequence,  it 
would  simplify  the  mass  spectral  analysis  and  speed  up  the  sequencing  process. 

A majority  of  the  mass  spectrometric  research  involving  peptide  sequencing  has 
been  done  using  either  the  collisionally  activated  dissociation  (CAD)  method  or  (in  FTICR 
studies)  the  sustained-off-resonance  collisionally  activated  dissociation  (SORI-CAD) 
method.  Although  both  of  these  can  provide  information  useful  in  elucidating  the 
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structure  of  a peptide,  each  has  disadvantages  associated  with  it.  CAD  is  a high-energy 
dissociation  method  which  usually  gives  too  much  fragmentation,  leading  to  very  complex 
spectra  which  can  be  a challenge  for  correct  interpretation.  For  both  CAD  and  SORI- 
CAD  methods,  peptide  fragmentation  is  caused  by  collision  with  an  externally  introduced 
neutral  gas  such  as  argon.  Introduction  of  neutral  gases  into  the  high  vacuum  region 
increases  the  analyzer  cell  pressure  and  reduces  the  optimal  low  pressure  environment  for 
obtaining  high  resolution  spectra. 

Previous  experiments  have  shown  that  infrared  multiphoton  dissociation  (IRMPD) 
is  a low-energy  dissociation  technique  that  works  well  for  structure  elucidation  of  small 
and  large  molecules. 80,118,119  Although  SORI-CAD  is  also  well  accepted  as  a low-energy 
dissociation  technique  yielding  fragmentation  pathways  of  low  activation  energy,  IRMPD 
combined  with  FT/MS  has  a number  of  advantages  over  SORI-CAD,  summarized  by 
Wilkins  and  co-workers.120  First,  the  photodissociation  (PD)  method  does  not  require 
introducing  collision  gas  into  the  trap,  thus  the  ideal  low  pressure  conditions  can  be 
maintained  for  optimal  instrumentation  performance.  Second,  product  ions  are  formed 
“on-axis”  (along  the  z-axis,  close  to  the  center  of  the  xy-plane)  without  any  significant 
magnetron  motion.  Third,  since  the  activation  of  precursor  ions  does  not  depend  upon 
imparting  translational  energy,  the  application  of  IRMPD  is  not  restricted  to  a certain  mass 
range  by  instrumental  dimensions,  which  limit  the  mass  range  for  collision  activation. 
Fourth,  because  the  excitation  of  precursor  ions  is  a stepwise  process,  energy  deposition 
increases  with  irradiation  time.  With  PD,  higher  energy  levels  for  dissociation  of  larger 
biological  systems  with  more  degrees  of  freedom  can  easily  be  accomplished  by  increasing 
the  irradiation  time.  One  disadvantage  of  the  IRMPD  method  is  that  it  requires  systems  to 
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contain  a chromophore  which  absorbs  10.6  pm  C02  laser  wavelength.  For  large 
biological  molecules,  this  requirement  can  be  easily  met  because  they  contain  hundreds  of 
amino  acid  residues  which  can  absorb  the  specific  wavelength,  but  small  peptides  with 
only  a few  residues  have  to  be  chosen  carefully  to  meet  that  requirement.  However,  the 
mass  spectra  generated  from  large  molecules  are  usually  very  complex  to  interpret,  so 
small  peptides  consisting  of  six  amino  acid  residues  or  fewer  were  chosen  in  the  present 
study.  Since  previous  studies  indicated  that  arginine-containing  peptides  absorb  thel0.6 
pm  fundamental  wavelength  of  C02  lasers,  such  peptides  were  chosen  to  be  the  target  of 
this  study.  Another  reason  for  choosing  arginine-containing  peptides  is  that  they  were 
previously  studied  with  other  methods  and  the  results  were  available  for 

* 117  121  122 

comparison. 

In  this  chapter,  the  fragmentation  patterns  of  small  arginine-containing  peptides 
(six  or  fewer  amino  acids)  studied  with  the  IRMPD  method  are  reported.  General 
correlations  between  the  dissociation  patterns  of  small  peptides  and  the  locations  of 
arginine  residues  in  the  peptides  were  investigated;  the  dissociation  spectra  observed  for 
singly  and  doubly  charged  states  of  the  same  peptide  were  analyzed  and  compared. 
Dissociation  of  the  same  peptides  was  also  studied  using  the  SORI-CAD  method,  and 
comparison  of  the  SORI-CAD  results  with  the  IRMPD  results  is  summarized. 

Experimental 


Materials  and  Sample  Preparation 

All  peptide  samples  used  in  this  project  were  purchased  from  Sigma  Chemical  Co.  (St. 
Louis,  MO)  and  used  without  further  purification.  All  the  ESI  samples  were  prepared 
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with  concentrations  of  ca.  1 x 10"6  M or  less  in  49/49/2  (v/v/v)  methanol/water/acetic  acid 
except  for  the  TPRK  and  TKPR  samples,  where  50/50  (v/v)  methanol/water  was  the 
solvent. 

Instrumentation 

The  shielded  9.4  T FTICR/MS  - IRMPD  experiments.  All  IRMPD 
experiments  were  performed  at  the  National  High  Magnetic  Field  Laboratory  in 
Tallahassee,  Florida  with  a shielded  9.4  T FTICR  mass  spectrometer  equipped  with  a 
Finnigan  Odyssey  data  system  as  shown  in  Figure  5-1. 123  The  system  consists  of  a 9.4  T 
superconducting  horizontal  bore  solenoidal  magnet  (Oxford  Instruments)  with  a 225  mm 
room  temperature  bore,  a cylindrical  homogenous  region  of  about  10  cm  diameter  by  10 
cm  axial  length,  and  a 90  mm  thick  iron  shield  which  reduces  the  stray  field  to  less  than  5 
gauss  at  a distance  of  30  cm  from  the  face  of  the  shield.  The  vacuum  of  the  instrument  is 
divided  into  four  sections.  The  first  section  is  the  electrospray  source,  which  is  a home- 
built  Chait-type  source.124  This  section  contains  a 20  cm  long  0.6  mm  i.d.  resistively- 
heated  stainless  steel  capillary  (Alltech,  Deerfield,  IL).  This  stage  is  evacuated  by  a 30 
cf7m  rotary  pump  (Varian,  Lexington,  MA)  to  1 torr.  The  second  section  contains  a 0.8 
mm  diameter  skimmer  (Beam  Dynamics,  Minneapolis,  MN)  and  a 60  cm  long  RF-only 
octopole  ion  guide,  and  this  section  is  pumped  to  a pressure  of  1 mtorr  by  a 1500  L/s 
hybrid  turbo-drag  pump  (Balzers,  Hudson,  NH).  The  third  section  of  the  vacuum  system 
contains  a second  octopole,  2 m in  length,  and  is  differentially  pumped  by  a 1500  L/s 
cryopump  (CTI  Cryogenics,  Waltham,  MA).  Both  octopoles  are  constructed  with  1.6  mm 
titanium  rods  with  a 4.8  mm  i.d.,  and  they  are  driven  at  1.3  MHz  by  high  power  RF 
supplies.  The  last  section  is  the  analyzer  region,  with  a base  pressure  of  1 x 10'8  torr. 


79 


which  contains  a 5 cm  diameter  cylindrical  open-ended  Penning  trap  constructed  with 
copper  plates  and  Macor  insulators.125 


Iron  Shield 


Open  Trap 


Ion 

Gauge 


Octupole  Ion  Guide 
/ Ion  Gauges 

— ESI  Source 


Cryopumps,  each  1500  L/s 


Hybrid  Turba/Drag 
Pump  1500  L/s 


Figure  5-1. 

The  9.4  Tesla  mass  spectrometer  at  NHMFL  used  for  IRMPD  experiments. 

All  ESI  samples  were  sprayed  at  a rate  of  0.3  pl/min.  After  the  ions  were  produced 
with  ESI  and  then  transferred  to  the  ion  trap  using  the  second  octopole,  the  ions  of 
interest  were  isolated  with  a combination  of  a chirp  pulse  and  SWIFT.  Then  a model  E48- 
2-115  Synrad  C02  laser  located  on  the  opposite  side  of  the  magnet  from  the  ESI  source 
was  fired  for  0.3  or  0.4  sec  with  output  power  between  20  to  40  W.  The  laser  firing  time 
and  the  laser  power  used  in  all  IRMPD  experiments  were  adjusted  to  produce  a reasonable 
amount  of  dissociation. 

The  shielded  4.7  T FTICR/MS  for  SORI-CAD  experiments.  All  CAD  and 

Sori-CAD  data  were  collected  using  a shielded  4.7  T FTICR  mass  spectrometer  (Bruker 
Analytical  Systems,  Inc.,  Billerica,  MA,  USA)  at  the  University  of  Florida.  A schematic 
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diagram  of  the  instrument  is  shown  in  Figure  5-2.  The  system  consists  of  a shielded 
horizontal  4.7  T superconducting  solenoidal  magnet  (Magnex  Scientific  Ltd.,  Abingdon, 
England)  with  a 150  mm  room  temperature  bore  and  a 170  mm3  cylindrical  RF-shimmed 
Infinity™  analyzer  cell.70 

The  vacuum  system  of  the  instrument  is  divided  into  three  sections.  The  first 
section  is  an  external  electrospray  source,72  which  utilizes  a hexapole  ion  guide  (Analytica 
of  Branford,  Branford,  CT,  USA).  The  ESI  source  was  modified  in  our  laboratory  with  a 
heated  metal  capillary  for  desolvation  purposes,  which  is  similar  in  design  to  that  reported 
by  Ikonomov  and  Kebarle.126  This  stage  is  evacuated  by  a 500  L/min  mechanical  pump 
and  a turbo-drag  pump  with  a mechanical  forepump  to  ca.  2 x 10'3  torr.  A 800  L/s 
cryopump  (liquid  He)  is  used  to  pump  the  region  between  the  hexapole  and  the  first 
conductance  limit  to  give  a source  pressure  of  about  2x10^  torr.  The  second  section 
contains  all  the  electrostatic  ion  transfer  optics  and  is  pumped  by  a 400  L/s  cryopump 
(liquid  He).  The  third  section  of  the  vacuum  is  the  ultrahigh  vacuum  region  which 
contains  the  analyzer  cell,  and  it  is  pumped  by  a 400  L/s  cryopump  (liquid  He)  to  maintain 
a base  pressure  in  the  low  10'9torr  range.  The  stainless  steel  capillary  (500  pm  i d.  x 170 
mm  length)  was  held  at  ground  potential  and  was  maintained  at  a temperature  of  120  °C. 
After  leaving  the  ESI  source,  ions  were  accelerated  by  a 3 kV  potential  difference  and 
focused  through  the  conductance  limits  before  being  decelerated  and  trapped  in  the  ion 
trap. 

The  same  ESI  samples  used  for  IRMPD  investigations  were  used  for  SORI-CAD 
experiments.  They  were  sprayed  with  a flow  rate  of  1 pl/min.  The  mass  spectra  of  the 
molecular  ions  were  collected  with  a cell  pressure  of  1 x 10'9  torr,  and  the  SORI-CAD 
experiments  were  performed  with  Ar  as  the  collision  gas  at  a pressure  of  1 x 10'7  torr. 
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4.7  T superconducting  magnet 


Figure  5-2. 

4.7  Tesla  mass  spectrometer  used  for  all  SORI-CAD  experiments  at  UF. 

Results  and  Discussion 

Nomenclature  of  Peptide  Fragmentation  in  the  Gas  Phase 

With  the  capability  of  mass  spectrometry  for  sequencing  peptides  and  proteins,  it 
became  necessary  to  develop  a systematic  naming  system  for  peptide  fragments  produced 
in  gas  phase.  The  first  comprehensive  nomenclature  was  developed  by  Roepstorff127  at 
Odense  University  in  Denmark,  and  it  was  derived  from  existing  naming  systems  used  by 
several  groups.128,129  Since  Roepstorff  proposed  his  nomenclature  in  1984,  it  has  been 
refined  and  extended  by  several  groups130,131  to  account  for  all  possible  peptide  fragments 
seen  in  the  gas  phase.  In  the  following  section,  the  nomenclature  used  for  the  present 
work  will  be  briefly  summarized. 

Amino  acids  are  the  building  blocks  of  peptides  and  proteins,  and  they  are  based 
on  linear  arrangements  of  repeating  units,  -NH-CHR-CO-,  which  differ  only  in  the  nature 
of  the  side  chain,  R.  When  peptides  or  proteins  undergo  collisional  dissociation  in  the  gas 
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phase,  at  least  three  types  of  fragmentation  paths  are  observed:  N-terminal  peptide 
fragmentation,  C-terminal  peptide  fragmentation,  and  nonsequential  fragmentation.  The 
different  fragmentation  types  are  denoted  a„,  bn,  cn,  and  dn  when  the  positive  charge  is 
retained  on  the  N-terminal  peptide  fragment  and  vn,  wn,  yn",  and  z„  when  the  positive 
charge  is  retained  by  the  C-terminal  peptide  fragment.  When  the  positive  charge  is  located 
on  a mid-sequence  fragment,  the  fragment  is  referred  as  nonsequential  fragmentations. 

1.  N-terminus  charge  retention: 

a.  The  a„,  bn,  and  cn  ions  are  formed  by  main-chain  fragmentation: 


H2N-CHRrCO-NH-CHR2-CO-NH-CHR3-  CO-NH-  CHR,-...-CHR0  -C02H 


2H 


b.  The  cn"  ions  are  formed  by  attaching  two  H atoms  to  the  corresponding  cn  ions. 
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c.  The  dn  ions  are  formed  by  homolytic  cleavage  of  the  bond  between  the  P and  y 
carbons  of  the  side  chain  of  a„+1  ions. 


dn 


2.  C-terminus  charge  retention: 

a.  The  x„,  yn,  and  z ^ ions  are  formed  by  peptide  backbone  fragmentation: 


H,N-CHR<-. . . -CHRd-CO-NH-CHR, 


y3 


CO-NH 


CHR2-C0-NH-CHR,-C02H 


2H 
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b.  The  yn"  ions  are  formed  from  adding  two  H atoms  to  the  corresponding  yn  ions. 


yn 


// 


c.  The  wn  ions  are  produced  by  homolytic  cleavage  of  the  bond  between  the  P and 
y carbons  of  the  side  chain  of  z n+1  ions: 


1 


wt 
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d.  The  vn  ions  are  produced  by  eliminating  the  side  chain  of  ynions  along  with  an 
H-atom: 


v 


n 


3.  Nonsequential  fragmentation:  internal  acyl  ions  and  immonium  ions: 


H 

H2N-CHR1-...-CHRni-CO 


^Single  a.a. 
Immonium 

NH-CH^,- 


+ 


Acyl 


CO-NH-CHR^-CO-NH. . -CHR,,  -C02H 


Immonium 


+ 


R 

© I 

H-N—  C 

I 

H 


Single  a.a.  immonium  ion 
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IRMPD  Results 

Table  5-1  summarizes  the  abbreviations  used  for  amino  acid  residues,  and  Table 
5-2  lists  the  peptides  used  in  this  study.  A typical  IRMPD  spectrum  collected  with  the  9.4 
T mass  spectrometer  at  the  NHMFL  is  given  in  Figure  5-3,  which  shows  the  IRMPD 
spectrum  of  the  RGD  peptide.  The  important  peaks  on  the  spectrum  are  labeled 
according  to  the  nomenclature  discussed  above.  For  example,  the  molecular  ions  MH+ 
occur  at  m/z  = 347.2,  and  the  fragment  ions  MH+  - NH3  at  m/z  = 330. 1 . The  other  main 
fragments  seen  are  b2  + H + OH  at  m/z  = 232. 1,  b2  at  m/z  = 214.1,  and  b2  - NH3  at  m/z  = 
197. 1.  It  is  interesting  to  notice  that  the  RGD  peptide  only  produces  b type  ions  in  the  gas 
phase. 


mass-to-charge  ratio 


Figure  5-3. 

IRMPD  spectrum  of  RGD  obtained  on  the  9.4  T FTICR  system  at  the  NHFML. 
The  major  peaks  on  the  spectrum  are  labeled. 
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Table  5-1.  Amino  acid  residues 


Amino  acid 

Threetone)  letter  code 

MonoisotoDic  mass 

Averaee  mass 

Glycine 

giy  (G) 

57.021 

57.052 

Alanine 

Ala  (A) 

71.037 

71.079 

Serine 

Ser (S) 

87.032 

87.078 

Proline 

Pro  (P) 

97.053 

97.117 

Valine 

Val  (V) 

99.068 

99.133 

Threonine 

Thr  (T) 

101.048 

101.105 

Cysteine 

Cys  (C) 

103.009 

103.145 

Isoleucine 

lie  (I) 

113.084 

113.160 

Leucine 

Leu  (L) 

113.084 

113.160 

Asparagine 

Asn  (N) 

114.043 

1 14.104 

Aspartic  acid 

Asp  (D) 

115.027 

115.089 

Glutamine 

Gin  (Q) 

128.059 

128.131 

Lysine 

Lys  (K) 

128.095 

128.174 

Glutamic  acid 

Glu  (E) 

129.043 

129.116 

Methionine 

Met  (M) 

131.040 

131.199 

Histidine 

His  (H) 

137.059 

137.142 

Phenylalanine 

Phe (F) 

147.068 

147.177 

Arginine 

Arg  (R) 

156.101 

156.188 

Tyrosine 

Tyr  (Y) 

163.063 

163.176 

Tryptophan 

Trp  (W) 

186.079 

186.213 

Table  5-2.  List  of  peptides  studied 

Peptide  # 

Peptide 

Average  mass  (Da) 

1 

RGD 

346.2 

2 

RGDS 

433.4 

3 

RGES 

447.4 

4 

RHF 

578.6 

5 

GRGD 

403.4 

6 

GRGDS 

490.5 

7 

PQRF-amide 

545.6 

8 

GHRP 

465.5 

9 

TPRK 

500.6 

10 

SDGRG 

490.5 

11 

TKPR 

500.6 

12 

DSDPR 

588.6 

13 

RFDS 

523.5 

14 

YIGSR 

594.7 

15 

GPRP 

425.5 

16 

RGPFPI 

685.8 
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Figure  5-4  displays  the  IRMPD  spectrum  of  RGD  which  was  obtained  using  a 3 T 
mass  spectrometer  at  UF,  and  it  can  be  seen  that  the  results  obtained  with  both 
instruments  are  very  similar.  This  demonstrates  that  the  IRMPD  results  are  reproducible. 
Although  both  mass  spectrometers  provided  essentially  the  same  information,  the  IRMPD 
results  reported  in  this  chapter  were  acquired  using  the  9.4  T spectrometer. 

Summary  of  the  IRMPD  results  for  Arg-containing  peptides.  It  has  been 
noted  by  many  research  groups  that  the  presence  and  location  of  certain  amino  acid 
residues,  especially  Arg,  influence  the  fragmentation  spectra  of  proteins  or  peptides.33,132 
Arginine  has  been  identified  as  the  most  basic  amino  acid  residue  in  the  gas  phase,133  and 
its  side  chain  has  been  pinpointed  as  the  most  favored  site  of  protonation  for  singly 
charged  ions.122  Therefore,  most  of  the  ions  produced  in  a dissociation  spectrum  for  an 
arginine-containing  peptide  contain  the  Arg  residue.  This  can  be  illustrated  by  the  results 


Figure  5-4. 

IRMPD  spectrum  of  collected  with  3 T FTICR  system  at  UF. 
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for  arginine-containing  peptides  summarized  in  Table  5-3.  The  table  lists  the  three  most 
intense  fragment  ions  seen  for  each  peptide  and  is  arranged  in  order  as  the  location  of  Arg 
moves  from  the  N-terminus  to  the  C-terminus  of  the  peptides.  As  shown  in  the  table, 
three  correlations  between  the  fragmentation  pattern  and  the  location  of  the  Arg  amino  acid 
residue  within  the  peptide  sequence  can  be  derived  from  the  results,  (i)  bn  type  ion  series 
occur  as  the  main  fragment  ions  when  the  Arg  residue  is  located  at  N-terminus  of  the 
peptide;  (ii)  y„"  type  ion  series  are  seen  as  the  primary  fragments  when  Arg  is  located  at  the 
C-terminus  of  the  peptide  sequence;  and  (iii)  both  bn  and  yn"  type  ion  series  are  observed 
when  Arg  is  located  in  the  middle  of  a peptide  sequence.  Correlations  (i)  and  (ii)  agree 
well  with  the  results  reported  by  Biemann’s  group.131,134  However,  the  IRMPD  results 
reported  above  are  different  from  the  fragment  ions  produced  from  high-energy  CAD 
experiments.  As  Table  5-3  shows,  the  low  energy  IRMPD  dissociation  technique 
generates  predominately  bn  and  yn"  ion  type  ions  which  result  from  main-chain 
fragmentation  of  peptides.  The  fragmentation  ions  observed  from  CAD  experiments  are 
predominantly  vn  and  wn  ions  which  contain  an  Arg  at  the  C-terminus  and  a„  and  dn  ion 
series  resulting  from  an  N-terminal  Arg.135  As  the  definitions  of  vn,  wn,  and  dn  ions  indicate, 
they  are  produced  by  cleaving  the  side  chains  of  yn",  z ^ and  a,,  ions  respectively.  This 
means  that  interpreting  CAD  spectra  involves  an  extra  step  of  calculation  of  side  chain  loss 
compared  to  the  straightforward  analysis  of  an  IRMPD  spectrum  since  the  IRMPD 
technique  only  produces  peptide  backbone  fragments.  The  general  trend  (iii)  observed  for 
Arg-containing  peptides  also  agrees  well  with  the  results  obtained  by  other  groups.33  The 
IRMPD  studies  show  that  the  three  correlations  enable  one  to  distinguish  Arg-containing 
peptides  of  different  sequences  and  same  mass  by  using  MS  in  gas  phase,  and  the  results 
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also  imply  that  a terminal  Arg  residue  provides  for  optimal  sequencing  of  a peptide  due  to 
the  fact  that  only  two  ion  series,  yn"  and  bn,  are  produced. 


Table  5-3. 

IRMPD  Results  of  Arginine-containing  Peptides 


Peptide 

IRMPD  results 

Peptide 

IRMPD  results 

RGD 

b2 

PQRF-amide 

MH+-  C2H5NHC(NH2)NH 

MH+-  NH3 

b3 

b2  + H + OH 

y2" 

RGDS 

b3 

GHRP* 

b3  - h2o  - nh3 

MH+-  H20 

b3 

b3  -H20 

MH+-  H20  - NH3 

RGES 

MH+-  H20 

TPRK* 

b3  + H + OH 

b3  - h2o 

x2 

b3 

b3 

RHF* 

b2 

SDGRG 

y3" 

y2 "-  H20  (very  weak) 

y2"-  h2o  - nh3 

MH+-  2H20 

GRGD 

MH+-  H20  -NH3 

TKPR* 

MH-  C2H5NHC(NH2)NH 

b2  + H20 

b3 

b3  - NH3 

y2" 

GRGDS 

b4 

DSDPR 

MH+ 

MH+-  H20 

y 2" 

MIT-  H20  -NH3 

y2  - h2o 

* peptide  contains  two  basic  amino  acids 
bn  and  yn"  type  ions  are  shown  in  bold. 


Singly  charged  ions  ([M+H]+)  vs.  doubly  charged  ions  ([M+2H]2+).  Both  singly 
(1+)  and  doubly  (2+)  charged  ions  were  observed  with  strong  intensities  for  the  peptides 
PQRF-amide,  DSDPR,  and  RGPFPI.  After  each  charge  state  of  the  peptides  was 
individually  isolated  in  the  analyzer  cell,  it  was  studied  with  the  IRMPD  method.  The 
IRMPD  spectra  obtained  are  given  in  Figures  5-5,  5-6,  and  5-7,  and  the  results  are 


summarized  in  Table  5-4. 
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The  two  IRMPD  spectra  of  the  [M+H]+  and  [M+2H]2+  ions  of  DSDPR  shown  in 
Figure  5-5  were  collected  under  the  same  conditions  except  stronger  laser  irradiation 
power  was  used  to  dissociate  the  singly  charged  ions  than  the  laser  irradiation  power  used 
for  the  doubly  charged  ions.  The  spectrum  of  the  1+  ions  was  acquired  with  28  watts  of 
laser  output  power  while  the  spectrum  of  2+  ions  was  obtained  with  only  8 watts  of  laser 
output  power.  Figure  5-6  and  Figure  5-7  present  the  IRMPD  spectra  of  1+  and  2+  ions  of 
PQRF-amide  and  RGPFPI.  The  PQRF-amide  spectra  were  collected  with  24  watts  of  laser 
irradiation  power  while  the  spectra  of  RGPFPI  were  collected  with  higher  laser  output 
power  of  40  watts  at  10.6  pm  wavelength. 

After  examining  the  assignments  listed  in  Table  5-4  and  comparing  the  appearance 
of  the  two  spectra  of  the  same  peptide  in  Figures  5-5  to  5-8,  it  is  not  difficult  to  see  that  1 + 
and  2+  molecular  ions  produce  different  fragments  and  generate  different  IRMPD  spectra. 
For  example,  the  singly  charged  ions  of  PQRF-amide  produce  primarily  bn  ion  series  while 
the  doubly  charged  molecular  ions  produce  mainly  yn"  type  ions;  the  singly  charged 
molecular  ions  of  DSDPR  generate  predominately  yn"  ion  series  while  doubly  charged 
molecular  ions  give  a mixture  of  bn  and  yn"  ions;  the  singly  charged  molecular  ions  of 
RGPFPI  dissociate  to  yn"  ion  series  while  the  corresponding  doubly  charged  ions  provide 
basically  bn  type  ions.  However,  there  is  very  little  dissociation  observed  for  the  singly 
charged  RGPFPI  ions.  Although  different  charge  states  of  the  same  peptide  generate 
different  spectra,  the  fragment  ions  observed  in  both  spectra  are  still  dominated  by  only  yn" 
and  bn  fragment  series,  as  observed  for  other  peptides  in  this  study,  and  reported  by  Hunt’s 

130 


group. 
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Figure  5-5. 

IRMPD  spectra  of  1+  and  2+  charge  states  of  DSDPR  with  both  spectra  collected 
under  the  same  conditions  except  the  laser  output  power.  The  1+  spectrum  was 
collected  with  28  watts  of  laser  output  power;  the  V spectrum  was  collected  with  8 
watts  of  laser  output  power  at  1064  pm  wavelength. 
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Figure  5-6. 

IRMPD  spectra  of  T and  2+  charge  states  of  PQRF-amide  with  both  spectra 
collected  under  the  same  conditions:  24  watts  of  laser  output  power  at  1064 
wavelength. 
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Figure  5-7. 

IRMPD  spectra  of  1+  and  2+  charge  states  of  RGPFPI  with  both  spectra  collected 
under  the  same  conditions:  40  watts  of  laser  output  power  at  1064  pm  wavelength. 
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Table  5-4. 

Assignments  for  ions  produced  from  both  1+  and  2+  charge  states  of  PQRF-amide, 
DSDPR,  and  RGPFPI  peptides  as  shown  in  Figure  5-4. 


Peptide 

IRMPD  fragments 

PQRF-amide 

546.3 

MH' 

529.3 

MH'-  NH3 

(1+  charge  state) 

511.3 

MH*-  NH3  - H20 

459.2 

MH*  -CH3CH2NHC(NH2)(NH)* 

382.2 

b3 

365.2 

b3  - NH3 

364.2 

b3  - H20 

348.2 

b3  - NH3  - H20 

347.2 

b3  - H20  - NH3 

321.2 

Yz" 

PQRF-amide 

546.3 

MH' 

449.3 

y3" 

(2+  charge  state) 

432.2 

y3"  - nh3 

321.2 

y2" 

273.7 

(M  + 2H)2+ 

265.1 

(M  - NH3  + 2H)2+ 

DSDPR 

589.3 

MH* 

(1+  charge  state) 

272.2 

Yz" 

254.1 

y2"-  H20 

197.7 

noise 

DSDPR 

589.3 

MH* 

(2+  charge  state) 

474.2 

y/' 

387.2 

y3" 

318.1 

b3  (very  weak) 

300.1 

b3  - H20 

295.1 

(M+2H)2+ 

286.1 

(MH2  - H20)2+ 

272.2 

y2" 

175.1 

y." 

*CH3CH2NHC(NH2)(NH)  is  the  side  chain  of  arginine  residue. 
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Table  5-4— continued: 


Peptide 

IRMPD  fragments 

RGPFPI 

686.4 

MFT 

573.3 

c5" 

(1+  charge  state) 

530.3 

y5" 

495.2 

y5"  - H20  - NH3 

458.2 

b4 

376.2 

Ys" 

302.1 

c2"  +H 

RGPFPI 

555.3 

b5 

526.3 

b5  - HCO  from  Proline  at  C-terminal 

(2+  charge  state) 

509.3 

b5  - NH2CH2NH2  from  Arg  side  chain 

495.2 

b5  - H20  - NH=C=NH  (60) 

458.2 

b4 

441.2 

b4  - NH3 

430.2 

c r 

i 

O 

O 

413.2 

b4  - NH3-CO 
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By  comparing  the  results  generated  from  1 + and  T molecular  ions  of  the  same 
peptide,  it  can  also  be  concluded  that  the  doubly  charged  molecular  ions  experience  a 
greater  extent  of  dissociation  than  the  corresponding  singly  charged  ions.  This  can  be  seen 
for  all  three  peptides  whose  spectra  are  given  in  Figures  5-5  to  5-8.  The  greater  extent  of 
fragmentation  for  the  higher  charge  state  than  the  lower  charge  state  is  consistent  with  the 
charge/charge  repulsion  explanation  of  why  multiply-charged  proteins  readily  fragment  in 
the  gas  phase. 

Peptide  fragmentation  vs.  laser  power.  The  extent  of  fragmentation  of  a peptide 
versus  the  laser  irradiation  power  was  studied  with  the  randomly  chosen  GRGDS,  RGES, 
and  RGDS  peptides.  The  relative  intensities  of  the  main  fragment  ions  observed  for  each 
peptide  are  plotted  against  the  laser  output  power  used  for  producing  the  spectra  in 
Figures  5-9,  5-10,  and  5-11.  The  plots  demonstrate  that  the  relative  intensities  of  the 
major  fragment  ions  of  each  peptide  increase  with  the  laser  power.  Although  all  peptides 
studied  show  a greater  extent  of  dissociation  with  higher  laser  power,  the  fragment  ions 
remain  as  the  main-chain  cleavage  ions.  As  the  laser  power  increases,  the  losses  of  NH3 
and  H20  from  fragment  ions  are  common  to  see  as  shown  by  the  plot  of  RGDS  given  in 
Figure  5-9.  For  example,  both  b3  - H20  and  b3-H20-NH3  ions  are  seen  along  with  the  b3 
ion  at  higher  laser  power.  Since  the  intensities  of  two  newly  formed  ions  from  b3  ions 
increase  with  the  laser  power,  the  growth  of  intensity  of  b3  ions  is  slowed.  Therefore,  in 
IRMPD  experiments  the  proper  laser  power  needs  to  be  controlled  for  producing  maximal 
amount  of  low-energy  dissociation  fragments. 

TPRK  vs.  TKPR.  When  two  strong  basic  amino  acid  residues,  K and  R,  are 
present  in  the  same  peptide,  such  as  TPRK  and  TKPR,  the  fragmentation  pattern  does  not 
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Fragmentation  vs.  laser  power 


b2 

■ o b< 


Figure  5-9. 

Plots  of  relative  intensities  for  major  fragment  ions  b2  and  b4  of  GRGDS  observed 
in  the  IRMPD  spectra  vs.  laser  output  power  of  the  experiments. 
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Fragmentation  vs.  laser  power 


Figure  5-10. 

Plots  of  relative  intensities  for  major  fragment  ions  b2,  b3,  and  b3-H20  of  RGES 
observed  in  the  ERMPD  spectra  vs.  laser  output  power  of  the  experiments. 
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Fragmentation  vs.  laser  power 


laser  output  power  (watt) 
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b3  - H20  - NH3 

mh*-h2o 

Figure  5-11. 

Plots  of  relative  intensities  for  major  fragment  ions  b2,  b3  b3-H20,  b3-H20-NH3,  and 
MH+-  H20  of  RGDS  observed  in  the  IRMPD  spectra  vs.  laser  output  power  of  the 
experiments. 
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always  strictly  follow  the  predictions  from  earlier  in  this  chapter:  bn  type  of  fragment  ions 
are  observed  for  N-terminal  Arg-containing  peptides  and  yn"  type  of  ions  are  generated  for 
C-terminal  Arg-containing  peptides.  Discrepancies  can  be  seen  from  the  results  for  TPRK 
and  TKPR  peptides  summarized  in  Table  5-5.  The  table  shows  that  TKPR  produces  both 
bn  and  yn"  ion  series  rather  than  only  yn"  type  of  ions  which  are  seen  for  C-terminal  Arg- 
containing  peptides.  For  the  TPRK  peptide,  x„  ions  are  observed  along  with  bn  ions, 
instead  of  the  mixture  of  yn"  and  bn  ions  obtained  for  peptides  with  Arg  in  the  mid- 
sequence. This  observation  implies  that  the  fragmentation  mechanism  of  these  peptides 
with  more  than  one  strong  basic  residue  is  different  from  that  of  peptides  containing  only 
one  basic  amino  acid.  The  different  fragmentation  mechanism  is  probably  due  to  the 
availability  of  more  than  one  protonation  site  in  the  amino  acid  sequence. 

Although  TPRK  and  TKPR  contain  exactly  the  same  amino  acid  residues  with  only 
different  orders  of  arrangement,  the  appearances  of  their  IRMPD  spectra  are  completely 
different,  as  shown  in  Figure  5-12.  Even  though  both  spectra  were  collected  with  20 
watts  of  laser  output  power  and  the  same  experimental  conditions,  TPRK  obviously 
dissociates  to  a greater  extent  than  TKPR.  Since  TPRK  contains  two  strong  basic  amino 
acid  residues,  arginine  and  lysine,  right  next  to  each  other  in  the  peptide,  stronger 
protonation  competition  between  the  basic  R-groups  of  the  amino  acid  residues  could 
occur.  Since  the  charge/charge  repulsion  mechanism  guides  the  peptide  fragmentation  in 
the  gas  phase,  it  is  not  surprising  that  a greater  extent  of  fragmentation  occurred  for 
peptides  containing  higher  number  of  basic  amino  acid  residues  and  for  peptides  with  their 
basic  amino  acid  residues  located  close  to  each  other  within  the  sequence. 
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Table  5-5. 

List  of  assignments  for  ions  observed  in  the  IRMPD  spectra  of  TPRK  and  TKPR 
peptides. 


Peptide 

IRMPD  fragments 

TPRK 

501.3 

MH+ 

483.3 

MH+  -H20 

457.3 

MH+  -C02 

439.3 

MH+  -C02  -H20 

373.2 

b3  + H + OH 

355.2 

b3 

337.2 

b3  - H20 

329.2 

*2 

320.2 

b3  - H20  - NH3 

311.2 

x2  - H20 

293.2 

x2  - 2H20 

276.1 

x2  - 2H20  - NH3 

TKPR 

501.3 

MH+ 

483.3 

MH+  -H20 

457.3 

MH+  -C02 

439.3 

MH+  -C02  -H20 

379.2 

unable  to  assign 

345.2 

b3  + H + OH 

327.2 

b3 

309.2 

b3  - H20 

272.2 

y2" 

254.2 

y2"  - h2o 

230.1 

b2 

212.1 

b2  - H20 
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Figure  5-12. 

The  IRMPD  spectra  of  TPRK  and  TKPR  peptides.  Both  spectra  were  collected 
with  same  laser  output  power  (20  watts)  and  exactly  the  same  experimental 
conditions:  a)  TPRK,  b)  TKPR. 
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SORI-CAD  vs.  IRMPD 

The  sustained-off-resonance  collisionally  activated  dissociation  (SORI-CAD) 
method  is  another  well  accepted  low-energy  dissociation  technique  for  protein  or  peptide 
sequencing  in  gas  phase.  Dissociation  of  some  peptides,  such  as  SDGRG,  GPRP,  and 
YEGSR,  was  studied  with  a pressure  of  10"7  torr  of  Ar  as  the  collision  gas  using  the  SORI- 
CAD  method  on  the  4.7  T instrument  at  UF,  and  the  results  obtained  were  compared  with 
the  corresponding  ERMPD  results.  Although  most  of  the  results  obtained  from  both 
approaches  are  compatible  (see  Figures  5-13,  5-14,  and  5-15),  the  IRMPD  method  can 
provide  more  fragmentation  information  for  some  peptides  such  as  the  GPRP  peptide  given 
in  Figure  5-13.  The  IRMPD  method  not  only  provides  sufficient  fragment  ions,  it  also 
provides  a cleaner  spectrum  than  the  SORI-CAD  method,  especially  in  the  low  mass 
region,  as  shown  in  the  spectrum  of  RGES  in  Figure  5-14. 

It  has  been  demonstrated  that  if  experimental  conditions,  including  neutral  gas 
pressure,  reaction  delay,  and  excitation  power  are  well  controlled,  the  SORI-CAD  method 
can  produce  a very  similar  spectrum  to  the  one  generated  using  the  IRMPD  approach. 
However,  if  one  of  the  parameters  is  not  optimized  properly,  the  resultant  spectrum  can  be 
very  different  from  the  ERMPD  spectrum.  When  the  collision  gas  pressure  is  too  high,  the 
reaction  delay  is  too  long,  or  the  excitation  power  is  too  strong,  complicated  spectra  can 
arise  from  excessive  dissociation.  For  example.  Figure  5-15  shows  a complicated 
spectrum  obtained  for  the  peptide  RGES  using  SORI-CAD  technique.  It  is  obvious  that 
the  many  peaks  seen,  especially  in  the  low  mass  region,  are  much  more  difficult  to  analyze 
than  a simpler  ERMPD  spectrum  (see  Figure  5-15). 
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a).  S ORI-C AD  spectrum  of  SDGRG 


mass-to-charge  ratio 


b).  IRMPD  spectrum  of  SDGRG 


Figure  5-13. 

Comparison  of  an  IRMPD  spectrum  vs.  a SORI-CAD  spectrum  obtained  for  the 
SDGRG  peptide,  a).  ERMPD  spectrum  obtained  with  the  9.4  T instrument  at 
NHMFL,  b).  SORI-CAD  spectrum  obtained  with  the  4.7  T instrument  at  UF. 
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a).  SORI-CAD  spectrum  of  GPRP 


b).  IRMPD  spectrum  of  GPRP 
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Figure  5-14. 

Comparison  of  SORI-CAD  and  IRMPD  spectra  of  the  GPRP  peptide, 
a).  SORI-CAD  spectrum,  b).  IRMPD  spectrum. 
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a).  SORI-CAD  spectrum  of  RGES,  -1000  Hz  off 


b).  IRMPD  spectrum  of  RGES 


Figure  5-15. 

Comparison  of  IRMPD  and  SORI-CAD  spectra  of  the  RGES  peptide, 
a).  SORI-CAD,  b).  IRMPD 
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This  study  has  shown  that  the  SORI-CAD  and  IRMPD  methods  can  provide 
essentially  the  same  information  if  the  experimental  parameters  of  both  methods  are 
properly  controlled.  However,  once  one  has  bought  the  laser  and  has  properly  aligned  it, 
the  only  experimental  parameters  that  need  to  be  adjusted  are  the  irradiation  time  and  the 
laser  power.  For  SOR1,  every  dissociation  experiment  requires  that  the  neutral  gas 
pressure  and  excitation  power  be  properly  controlled.  Therefore,  the  IRMPD  approach 
seems  to  be  the  ideal  choice  for  peptide  sequencing  in  the  gas  phase. 

Conclusions 

Studies  of  IRMPD  of  Arg-containing  peptides  in  the  gas  phase  reveal  the 
importance  of  gas-phase  fragmentation,  and  the  correlations  can  be  used  to  distinguish  the 
location(s)  of  Arg  in  a peptide,  which  can  lead  to  rapid  peptide  sequencing.  The  general 
relationships  between  the  peptide  fragmentation  patterns  in  the  gas  phase  and  the  location 
of  an  arginine  residue  in  a peptide  can  be  summarized  as:  bn  type  ions  are  observed  as  the 
primary  fragmentation  series  when  the  Arg  residue  is  located  at  the  N-terminus  of  the 
peptide;  yn"type  ions  are  seen  as  the  principal  fragmentation  series  when  the  Arg  residue  is 
located  at  the  C-terminus  of  the  peptide;  and  a mixture  of  mainly  bn  and  yn"  type  ions  is 
observed  in  the  dissociation  spectra  when  the  Arg  residue  is  located  in  the  middle  of  the 
peptide  sequence.  A second  major  observation  for  the  Arg-containing  peptides  is  that 
losses  of  H20  and  NH3  lead  to  prominent  fragment  peaks  in  the  spectrum. 

The  results  of  this  study  also  show  that  a greater  extent  of  fragmentation  is 
observed  for  higher  charge  state  molecular  ions  than  lower  charge  state  molecular  ions 
when  the  dissociation  spectra  of  1+  and  T charge  states  for  several  peptides  are  examined. 
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It  was  also  seen  that  the  fragmentation  pattern  is  not  always  strictly  followed  by  the 
peptides  that  contain  more  than  one  basic  amino  acid  residue  in  the  sequence.  When  the 
peptide  fragmentation  is  studied  as  a function  of  the  irradiating  laser  power,  the  results 
show  that  increasing  the  laser  power  generally  increases  the  amount  of  peptide  dissociation. 
In  conclusion,  this  study  has  shown  that  IRMPD  is  a true  low-energy  dissociation 
technique  when  compared  to  SORI-CAD  techniques,  producing  only  peptide  backbone 
fragmentation.  It  may  be  the  ideal  method  for  peptide  or  protein  sequencing  with  MS  in 
the  gas  phase. 


CHAPTER  6 

INVESTIGATING  CHEMICALLY  MODIFIED  RNASE  A WITH  ESI-FTICR/MS 

Introduction 

Over  the  past  ten  years,  developments  in  MS  have  made  it  possible  to  study  many 
biological  systems,  including  complexes  arising  from  noncovalent  interactions  between 
ligands  and  receptors,136,137  subunits  of  protein  complexes,138,139  and  oligonucleotides.24 
Other  compelling  biological  problems  lie  in  the  world  of  modified  enzymes.140  Although 
these  modified  enzymes  are  somewhat  useful,  their  rate  enhancements  do  not  match  those 
of  any  naturally  occurring  enzymes.  Thus,  improving  the  turnover  rate  of  a modified 
enzyme  has  been  an  important  goal.  One  means  of  achieving  this  goal  is  covalent 
modification  of  a naturally  occurring  enzyme  with  a transition  metal  cofactor.  This 
improves  the  binding  properties  of  the  enzyme  with  the  catalytic  enhancement  of  the 
cofactor. 

Due  to  its  feasibility  for  the  direct  analysis  of  multicomponent  proteolytic  peptide 
mixtures  (peptide  mapping),  MS  has  been  demonstrated  as  a powerful  method  in  the 
identification  of  covalent  post-translational  modifications141,142  and  of  multiple  chemical 
modification  sites.29  The  mass  spectrometric  peptide  mapping  method  has  been 
successfully  employed  to  investigate  specific  modification  reactions  such  as  amino- 
acylation,  carboxylate-amidation,  and  the  bifunctional  cysteine-modification  by  phenyl 
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Ill 


arsinoxide.143'145  The  experiments  reported  in  this  chapter  were  directed  at  following  the 
chemical  modification  process  of  ribonuclease. 

The  enzyme  ribonuclease,  chosen  for  this  study,  has  been  known  as  RNase,  RNase 
A,  RNase  I,  pancreatic  ribonuclease,  and  a pancreatic  enzyme.  It  is  the  most  stable  and 
well  characterized  enzyme  at  the  molecular  level.  The  124  amino  acid  residues  give  the 
protein  a molecular  weight  of  13,680  Da.  RNase  A contains  four  disulfide  bonds,  which 
occur  between  residues  26  and  84,  40  and  96,  58  and  1 10,  and  65  and  72,  as  shown  in 
Figure  6-1.  When  the  bridges  are  reduced,  the  protein  is  denatured  and  becomes  inactive. 
However,  it  is  possible  to  reduce  the  bridges  partially  and  retain  the  activity  of 
the  enzyme.  Removal  of  four  peptides  at  the  carboxyl  terminus  also  destroys  enzyme 
activity.  RNase  A is  an  endoribonuclease  which  cleaves  RNA  to  both  mononucleotides  and 
oligonucleotides  in  a two  step  process146  as  shown  in  Figure  6-2.  In  the  first  step,  a cyclic 
2',3'-phosphate  intermediate  is  produced,  due  to  intramolecular  attack  by  the  2'-OH  group 
of  the  ribose  moiety  on  the  phosphodiester  linkage.  The  second  step  involves  hydrolysis  of 
the  cyclic  2',  3 '-phosphate,  resulting  in  the  formation  of  a terminal  3-monophosphate  ester. 
The  enzyme  has  its  highest  activity  with  single  stranded  RNA,  and  it  has  very  low  activity 
with  double  stranded  RNA. 

The  His- 12,  His-1 19,  and  Lys-41  residues  have  been  shown  to  be  associated  with 
the  active  site  of  RNase  A,147  and  the  N-terminal  lysine  is  located  adjacent  to  Lys-41 . As 
the  3-dimensional  structure  in  Figure  6-3  shows,  the  N-terminal  Lys  is  more  exposed  to 
the  environment,  so  it  can  be  easily  accessed  and  acetylated.148  Therefore,  RNase  A seems 
to  be  an  ideal  choice  for  enzyme  modification  with  the  N-terminal  lysine  as  the  target  for 
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Figure  6-1. 

Amino  acid  structure  of  RNase  A.  The  four  disulfide  bonds  formed  from  eight 
cysteine  residues  are  shown  in  black. 
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Figure  6-2. 

Two-step  reaction  of  RNase  A with  RNA:  in  the  first  step,  a cyclic  2',3'-phosphate 
intermediate  is  produced,  due  to  intramolecular  attack  by  the  2'-0H  group  of  the 
ribose  moiety  on  the  phosphodiester  linkage;  the  second  step  involves  hydrolysis  of 
the  cyclic  2',  3 '-phosphate,  resulting  in  the  formation  of  a terminal  3'- 
monophosphate  ester 
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Figure  6-3. 

Ribbon  structure  of  RNase  A with  N-terminal  lysine  labeled  as  N;  the  four  disulfide 
bridges  are  shown  in  black  balls  and  sticks. 
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modification.  In  this  chapter,  the  initial  results  of  a study  of  purification,  modification, 
denaturation,  and  digestion  obtained  using  ESI-FTICR/MS  are  reported. 

Experimental 


Materials 

The  native  and  oxidized  RNase  A samples  were  purchased  from  Sigma  Chemical 
Co.  (St.  Louis,  MO.)  The  native  RNase  A was  purified  and  chemically  modified  with  the 
proper  modifiers,  and  all  ESI  samples  were  prepared  in  concentrations  of  about  1 x 10'6  M 
or  less  in  49/49/2  (v/v/v)  methanol/water/acetic  acid  by  Trisha  Joy  Mendez  in  the  Inorganic 
Division  of  Chemistry  Department  at  the  University  of  Florida.  The  HPLC  peptide 
standard  mixture,  used  for  identifying  the  HPLC  separated  peptides  with  MS,  included  five 
peptides:  Gly-Try  (238  Da),  Val-Try-Val  (379  Da),  methionine  enkephalin  acetate  (Try- 
Gly-Gly-Phe-Met,  573  Da),  leucine  enkephalin  (Try-Gly-Gly-Phe-Leu,  555  Da),  and 
angiotensin  II  acetate  (Asp-Arg-Val-Tyr-Ile-His-Pro-Phe,  1046  Da).  The  peptide  mixture 
was  also  purchased  from  Sigma  Chemical  Co.  (St.  Louis,  MO.) 

Mass  Spectrometer 

All  mass  spectrometric  experiments  were  performed  with  a shielded  4.7-Tesla  FT- 
ICR  mass  spectrometer  (Bruker  Analytical  Systems  Instruments,  Inc.,  Billerica,  MA)  at  the 
University  of  Florida.  A detailed  description  of  the  instrument  can  be  found  in  the 
experimental  section  of  chapter  5 of  this  dissertation.  The  capillary  is  a heated  metal 
capillary  which  was  modified  from  an  Analytica  electrospray  source  (Analytica  of 


Branford,  Branford,  CT). 
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All  data  were  collected  with  a cell  pressure  of  1 x 10'9  torr,  and  the  SORI-CAD 
experiments  were  performed  with  either  Kr  or  Ar  as  the  collision  gas  at  a pressure  of  1 x 
10"7  torr.  The  samples  were  sprayed  at  a rate  of  1 pl/min,  and  the  stainless  steel  capillary 
was  maintained  at  a temperature  of  120°C.  The  electrospray  conditions  were  tuned  to 
optimize  the  signals  of  interest.  The  electrospray  needle  was  located  at  a distance  typically 
about  1 cm  away  from  the  capillary,  which  provided  maximal  ion  current  and  optimal 
signals.  The  capillary  and  skimmer  voltage  were  adjusted  to  values  which  maximized  the 
signal  intensities  for  the  ions  of  interest,  and  typical  values  were  around  100  and  10  volts 
for  the  capillary  and  skimmer,  respectively. 

Molecular  Weight  (Mr)  Determination  for  Large  Biomolecules 

ESI/MS  has  been  well  known  for  its  ability  to  provide  accurate  molecular  weight 
determination  for  many  biological  systems,  including  bovine  insulin,  bovine  trypsin,  and 
bovine  carbonic  anhydrase.  Large  biomolecules  sprayed  with  the  electrospray  ionization 
technique  have  a distinctive  pattern  of  peaks  in  which  adjacent  peaks  always  appear  to  vary 
by  addition  or  subtraction  of  one  charge  for  the  ions  observed  in  the  spectrum,  and  the 
average  charge  state  is  observed  to  increase  in  an  approximately  linear  fashion  with  the 
average  molecular  weight  of  the  biomolecule  (M,).89 

The  M,  of  a macromolecule  can  be  easily  determined  from  mass  spectra  with  the 
general  assumption  that  adjacent  peaks  of  a charge  state  distribution  differ  by  only  one 
charge,  and  that  ions  are  produced  by  proton  attachment  to  the  molecular  ion.89  Although 
determination  of  IV^has  been  automated  by  many  data  systems,  the  basis  used  to  calculate 
the  M,  of  a protein  has  stayed  the  same.  In  order  to  understand  the  spectra  of  multiply 
charged  ions  in  ESI/MS,  it  is  important  to  understand  the  basis  for  determining  the 
molecular  weight  with  MS. 
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A typical  positive  ion  spectrum  of  a single  protein  consists  of  a series  of  peaks  with 
each  one  representing  the  molecular  ion  plus  a specific  number  of  protons.  The  ions  have 
the  general  form: 

(M  + nH)n+ 

where  M is  the  M,.  of  the  protein,  n is  an  integer  represents  the  number  of  protons  attached 
to  the  protein,  and  H is  the  proton  with  mass  of  1 .00794  Da.  Since  MS  measures  mass-to- 
charge  ratios  (m/z),  each  peak  seen  on  a typical  spectrum  can  be  expressed  as: 

m/z  = (M  + nH)  / n 

If  both  m/z  and  n are  known,  M can  be  readily  calculated  as: 

M = n (m/z  - H) 

If  either  M or  n is  known,  M and  n can  be  determined  from  two  equations  derived  from 
two  consecutive  peaks  (m/z),  and  (m/z)2  differing  by  one  proton  in  the  series, 

(m/z),  = (M  + n,H)/  n, 

(m/z)2  = (M  + n2H)  / n2 

where  (m/z),  < (m/z)2,  and  n2  = n,  - 1 . Solving  (m/z), and  (m /z)2  equations  yields: 

M = n2  ((m/z),  - 1) 
and 

n2  = ((m/z),  - 1)  / ((m/z)2  - (m/z),) 

For  example,  with  two  adjacent  peaks  in  a charge  state  distribution  of  a protein 
with  (m/z),  = 1223  and  (m/z)2=  1426, 

n2  = ((m/z),  - 1)  / ((m/z)2  - (m/z),)  = (1223  - 1)  / (1426  - 1223)  = 6.02 
the  charge  state  n2  is  calculated  as  6,  therefore  the  molecular  weight  of  the  protein  can  be 
determined  with  this  newly  determined  n2  value. 

M = n2  ((m/z),  - 1)  = 6 (1223  - 1)  = 7332  Da. 
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Results  and  Discussion 

The  overview  of  this  project  and  the  general  procedure  followed  are  illustrated  with 
the  flow  chart  in  Figure  6-4.  The  very  first  step  of  the  project  was  to  check  the  purity  of 
the  protein  of  interest  to  assure  that  it  was  initially  clean.  Given  that  the  sample  size  in  this 
project  is  quite  small  (100  pi  sample  with  concentration  of  10'6  M),  a quite  small 
contamination  of  the  RNase  A samples  could  be  a potential  hazard  which  could  block  the 
effectiveness  of  the  planned  chemical  modification.  Impurity  of  the  sample  can  also  lower 
the  effective  mass  resolving  power  of  the  mass  spectrometer  and  decrease  the  instrumental 
detection  limit.  Hence,  the  goal  of  the  very  first  step  was  to  examine  the  purity  of  the 
native  RNase  A sample  purchased  from  the  Sigma  Chemical  Co.  After  the  native  RNase  A 
was  purified  using  the  purification  steps  illustrated  in  the  flow  chart  of  Figure  6-5,  the 
second  step  was  to  chemically  modify  the  protein  with  proper  modifiers  and  confirm  the 
modification  by  measuring  the  molecular  weight  for  both  native  and  modified  RNase  A 
with  ESI-FTICR/MS.  At  this  point  it  is  important  to  know  that  the  protein  only  increases 
its  molecular  weight  by  the  molecular  mass  difference  of  one  modifier  molecule  from  its 
native  form  since  modifying  only  the  N-terminal  lysine  is  the  goal  of  the  project.  After  the 
native  protein  was  confirmed  to  be  singly  modified,  the  third  step  was  to  identify  the  exact 
modification  site  using  either  MS  alone  or  the  peptide  mapping  MS  technique.  In  the 
peptide  mapping  MS  approach,  the  protein  of  interest  is  first  enzymatically  digested  to 
yield  a peptide  mixture;  then  the  HPLC  technique  is  used  to  separate  the  peptide  mixture; 
next  MS  is  used  to  identify  the  HPLC  aliquots  and  to  sequence  the  peptides;  and  finally  the 
peptide  sequences  are  examined  carefully  and  are  used  to  derive  the  complete  sequence  of 
the  protein.  Although  different  dissociation  methods  such  as  Sori-CAD,  Capillary- 
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skimmer(CS),  and  photodissociation  (PD)  were  tried  to  fragment  the  gas  phase  proteins, 
no  method  provided  enough  useful  information  for  pinpointing  the  modification  site.  Since 
MS  alone  was  not  successful  in  sequencing  the  protein,  the  peptide  mapping  technique, 
which  was  demonstrated  as  a powerful  method  for  identifying  partial  protein  modifications, 
was  employed.  The  results  obtained  will  be  categorized  and  discussed  in  more  detail  below 
in  their  appropriate  section. 
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Figure  6-4. 

Overview  of  the  project  with  processes  leading  to  various  forms  of  the  protein 
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Figure  6-5. 

Process  used  for  purifying  the  RNase  A protein. 
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Native  RNase  A 

The  spectrum  of  the  unpurified  RNase  A sample  is  given  on  the  top  of  Figure  6-6 
and  the  spectrum  of  the  purified  sample  is  given  on  the  bottom  of  Figure  6-6.  Each 
spectrum  presented  in  Figure  6-6  is  the  sum  of  16  individual  transients  with  128k  data 
points  collected.  In  this  chapter,  while  the  horizontal  axis  represents  the  mass-to-charge 
ratio  for  every  spectrum,  the  vertical  axis  is  plotted  as  absolute  intensity  unless  otherwise  is 
specified.  The  impure  spectrum  clearly  shows  that  the  native  RNase  A has  contamination 
in  the  region  of  m/z  100-1500.  Comparison  of  the  two  spectra  indicates  that  a much 
cleaner  sample  emerged  following  the  purification  process.  Due  to  the  impurity  present  in 
the  sample,  the  top  spectrum  barely  shows  three  charge  states  while  the  purified  sample 
provides  five  charge  states.  The  purified  sample  shows  about  a factor  of  100  enhancement 
in  its  signal-to-noise  ratio.  The  spectrum  of  purified  RNase  A in  Figure  6-6  was  used  to 
determine  the  average  molecular  weight  of  RNase  A,  and  the  calculated  value  was  13,686 
kDa. 

One  distinct  feature  of  the  spectra  shown  in  Figure  6-6  is  that  every  charge  state  of 
the  protein  has  a few  peaks  which  correspond  to  phosphate  or  sulfate  adducts  with  a mass 
ca.  98  Da  higher  than  the  unaddicted  peaks.  It  is  common  to  see  complexes  with  1,  2,  and 
3 phosphate  or  sulfate  groups,  which  are  noncovalently  bound  to  the  protein  in  ESI/MS 
spectra  such  as  these.  For  example,  the  8+  charge  state  of  the  molecular  ion  corresponds  to 
m/z  171 1.8149,  and  the  molecular  species  containing  one  and  two  phosphate  or  sulfate 
groups  are  observed  at  m/z  171 1 .8149  and  1724.0463,  respectively.  Assignments  of  other 
signals  seen  in  the  spectra  are  listed  in  Table  6-1. 
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Figure  6-6. 

Spectra  of  impure  and  purified  RNase  A collected  under  the  same  experimental 
conditions  with  a background  pressure  of  10'9  torr  and  a sample  concentration  of 
10'6  M in  49/49/2  methanol/H20/acetic  acid.  Each  spectrum  is  the  result  of 
accumulation  of  16  individual  transient  with  128k  data  points.  The  top  spectrum 
shows  the  impure  sample  purchased  from  Sigma  Chemical  Co.,  and  the  bottom  one 
shows  the  purified  RNase  A sample. 
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Table  6-1. 

Mass  spectrometric  signals  observed  in  purified  native  RNase  A and  their 
assignments. 


m/z 

Assignment 

1245.1952 

(M  + 1 1 H)11+ 

1369.5558 

(M  + 10  H)10+ 

1521.7966 

(M  + 9 H)9+ 

1532.6397 

*(M  + 9 H + P)9+ 

1711.8149 

(M  + 8 H)8+ 

1724.0463 

(M  + 8H+Pf 

1956.8958 

(M  + 7 H)7+ 

1970.8195 

(M  + 7 H + P)7+ 

* P is  the  abbreviation  for  phosphate 


Discussion  of  Phosphate  Adducts 

Multiple  additions  of  roughly  98  Da  have  been  observed  accompanying  the 
molecular  ions  of  large  biomolecules,  as  are  seen  for  RNase  A in  Figure  6-6.  The  mass 
increase  of  roughly  98  Da  could  correspond  to  addition  of  a phosphate  or  sulfate  group 
since  both  of  them  have  similar  molecular  mass.  Because  most  of  the  published  literature 
reports  the  98  Da  adducts  as  phosphate  adducts,149  the  adducts  will  be  referred  as 
phosphate  group  throughout  this  dissertation.  As  previously  reported  by  many  groups,150 
the  number  of  phosphate  adducts  increases  as  the  charge  state  decreases  or  the  m/z  value 
increases.  Chait’s  group  has  reported  that  the  phosphate  ion  can  be  dissociated  from  the 
protein  by  treating  the  sample  with  barium  acetate,124  and  studies  by  Bruce  and  et  al.151 
show  that  the  phosphate  can  be  readily  dissociated  in  the  heated  capillary  inlet.  However, 
dissociation  by  collision  activation  in  the  capillary-skimmer  region  is  less  effective,  and  the 
dissociation  is  significantly  reduced  when  the  complexes  are  collisionally  activated  in  the 
ICR  cell  where  other  dissociation  processes  dominate. 
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Dissociation  Studies  of  the  Isolated  8+  Charge  State 

One  of  the  advantages  of  FT/MS  is  that  the  ions  of  interest  can  be  isolated  and 
stored  for  many  minutes,  if  desired,  while  other  reactions  with  the  isolated  ions  can  be 
carried  out  in  the  trap.  This  feature  assures  that  the  reactions  of  the  ions  of  interest  are  not 
influenced  by  the  contaminants  in  the  reaction  cell.  For  example,  the  8+  charge  state  with 
various  phosphate  adducts  of  native  RNase  A at  m/z  1711  can  be  isolated  (see  the  top 
spectrum  of  Figure  6-7),  and  the  complete  isotopic  distribution  of  the  8+  charge  state 
without  the  phosphate  adducts  can  also  be  isolated  (shown  in  the  bottom  spectrum  in 
Figure  6-7)  due  to  the  ultrahigh  resolving  power  of  the  FT/MS  instrument. 

After  the  8+  ions  are  isolated,  various  dissociation  methods  such  as  SOR1-CAD, 
photo-dissociation  with  a UV-visible  laser  (355  nm  and  532  nm),  and  capillary-skimmer 
dissociation  (CS-CAD)  were  tried  in  order  to  dissociate  the  8+  ions.  However,  the  amount 
of  fragmentation  of  the  protein  ions  was  not  sufficient  to  be  an  effective  aid  in  partial  gas 
phase  sequencing  (see  Figure  6-8).  The  low  extent  of  fragmentation  is  probably  due  to  the 
large  number  of  degrees  of  freedom  which  exist  in  RNase  A,  so  the  energy  absorbed  by  the 
system  from  the  dissociation  methods  is  not  enough  to  raise  the  internal  energy  high 
enough  to  cause  the  precursor  ion  to  dissociate. 

Singly  Modified  RNase  A 

Once  satisfactory  purity  of  the  RNase  A sample  was  achieved,  the  next  step  was  to 
chemically  modify  the  protein  at  the  N-terminal  lysine.  In  this  step,  it  is  important  to  insure 
that  the  modification  takes  place  at  only  one  site  on  the  protein;  in  other  words,  the 
molecular  mass  of  the  modified  protein  can  only  be  different  by  one  molecular  mass  of  the 


modifier  from  the  native  protein. 
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Figure  6-7. 

The  8+  charge  state  is  isolated  with  correlated  sweep  in  the  analyzer  cell.  a),  the  8+ 
ions  including  single  phosphate  adduct;  b).  the  8+  ions  with  no  phosphate  adducts. 
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Figure  6-8. 

A SORI-CAD  dissociation  spectrum  of  the  8+  charge  state  obtained  from  8 scans 
with  128k  data  size. 
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Modification  with  2-iminotheolane.  One  of  the  chemical  modifiers  tried  by  Ms. 
Mendez  was  2-iminotheolane  (C4H7NS-HC1,  mw  = 137.6  Da)  whose  reaction  mechanism 
with  a lysine  amino  residue  is  illustrated  in  Figure  6-9.  The  mechanism  indicates  that  2- 
iminotheolane  is  added  to  the 


Lysine  2-iminotheolane 

Figure  6-9. 

Reaction  mechanism  of  2-iminotheolane  with  a lysine  amino  acid  residue  of  a 
protein.  The  middle  structure  of  the  reaction  is  the  intermediate,  and  the  final 
product  is  shown  on  the  right. 

N-terminal  nitrogen  atom  of  the  lysine  residue.  If  only  the  N-terminal  lysine  is  modified, 
the  protein  should  increase  its  molecular  weight  by  approximately  137.6  Da, 
corresponding  to  one  2-iminotheolane  molecule.  However,  when  the  modified  solution 
was  examined  by  ESI-FTICR/MS,  the  mass  spectrum  collected  (see  Figure  6-10)  implied 
that  more  complicated  modifications  had  occurred  in  solution  than  were  expected.  Figure 
6-10  gives  the  mass  spectrometric  results  for  one  of  the  2-iminotheolane  modifications. 
The  top  spectrum  of  the  figure  shows  that  three  charge  states  of  9+,  8+,  and  T are  clearly 
seen  around  1500,  1700,  and  2000,  respectively,  with  many  extra  peaks  occurring  within 
each  charge  state  compared  to  the  spectrum  of  the  native  protein  given  in  Figure  6-6. 
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Since  the  8+  ions  are  the  most  intense  peaks  in  the  spectrum,  the  m/z  1700-1900  mass 
range  was  expanded  for  a closer  examination  of  the  individual  peaks.  This  is  shown  in  the 
bottom  two  spectra,  with  the  left  spectrum  being  the  expansion  for  all  three  charge  states 
and  the  right  one  being  the  expansion  for  only  8+  ions.  When  the  extra  peaks  were  closely 
analyzed,  they  actually  corresponded  to  different  combinations  of  2-iminotheolane  and 
phosphate  adduct  additions  to  the  native  RNase  A.  The  peaks  of  1725.28  could  be  the  8+ 
molecular  ion  with  one  phosphate  adduct,  1742.84  could  correspond  to  the  8+  ions  with 
addition  of  one  2-iminotheolane  molecule.  Although  the  2-iminotheolane  modification  did 
not  give  the  expected  single  modification  results,  it  would  provide  confidence  that  the 
RNase  A protein  could  be  modified,  even  singly,  with  proper  modifiers  under  the  proper 
conditions.  Since  2-iminotheolane  gave  more  modification  than  expected,  the  search  for 
the  right  chemical  modifier  continued. 

Modification  with  N-hydroxysuccimide  acetate  (NHS-acetate).  The  next 
modifying  reagent  tested  was  N-hydroxysuccimide  acetate  (NHS-acetate)  whose  reaction 
with  a N-terminal  lysine  amino  acid  residue  is  given  in  Figure  6-11. 152  As  the  reaction 
scheme  depicts,  the  NHS-acetate  reagent  works  differently  from  the  2-iminotheolane 
reagent.  NHS-acetate  adds  only  41  Da  (-OCCH  group  from  the  NHS-acetate  reagent)  to 
thee  nitrogen  atom  of  the  N-terminal  lysine  while  a whole  molecule  of  2-iminotheolane  is 
attached  to  the  N-terminal  nitrogen  atom  of  a N-terminal  lysine  residue  of  the  protein. 

The  addition  of  -OCCH  “neutralizes”  the  side  chain  of  the  Lys,  and  the  positive  charge  on 
the  nitrogen  atom  relocates  to  another  part  of  the  protein. 

The  41  Da  mass  addition,  implying  a single  modification,  was  detected  when  the 
NHS-acetate  modified  protein  sample  was  examined  with  ESI-FTICR/MS.  A typical  mass 
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Figure  6-10. 

A typical  mass  spectrum  of  the  2-iminotheolane  modification  on  the  RNase  A 
protein.  The  top  spectrum  shows  the  complete  mass-to-charge  region  between 
200  to  3000;  the  bottom  left  spectrum  gives  the  expansion  of  all  9+,  8+,  and 
7+charge  states;  and  the  bottom  right  spectrum  provides  a greater  expansion  for 
only  the  8+  charge  state. 
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Figure  6-11. 

Reaction  of  N-hydroxysuccimide  acetate  (NHS-acetate)  with  the  N-terminal  lysine 

amino  acid  residue  of  a protein. 

spectrum  collected  for  the  NHS-acetate  modified  RNase  is  shown  in  the  bottom  spectrum 
of  Figure  6-12,  and  its  corresponding  native  RNase  A spectrum  is  given  on  the  top  of  the 
figure.  Both  spectra  display  the  region  of  m/z  1300-2100,  where  the  protein  signals 
appear.  The  bottom  spectrum  of  the  modified  protein  clearly  exhibits  four  sets  of  signals 
which  correspond  to  the  charge  states  7+,  8+,  9+,  and  10+  at  m/z  around  1900,  1700,  1500, 
and  1300,  respectively.  When  the  spectrum,  especially  the  8+  charge  state,  is  carefully 
examined,  two  distinctive  series  of  signals  can  be  revealed  on  the  same  spectrum.  One  ion 
series  is  more  intense  than  the  other,  and  the  signals  of  one  series  occur  in  between  the 
signals  of  the  other  series,  with  the  observation  especially  pronounced  for  the  8+  ions.  The 
set  of  stronger  signals  at  m/z  1716.64,  1728.92,  and  1741.01  represent  the  modified 
proteins,  and  the  set  of  weaker  signals  occurring  at  m/z  1711.51  and  1723.52  correspond 
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to  the  native  proteins.  The  mass  difference  between  1716.64  of  the  modified  protein  and 
171 1.5  lof  the  native  protein  is  exactly  41  Da,  which  corresponds  to  the  exact  mass  of  a 
single  acetate  group.  The  addition  of  41  Da  implies  that  the  native  protein  is  singly  and 
successfully  modified  with  the  NHS-acetate  reagent.  Table  6-2  summarizes  the 
assignments  for  all  the  peaks  seen  in  the  spectrum  of  the  modified  peptide  shown  in 
Figure  6-12.  The  signals  of  the  charge  state  T are  not  listed  in  the  table  due  to  their  low 
intensities.  The  mass  spectrometric  results  show  that  RNase  A can  be  successfully  and 
singly  modified  with  the  NHS-acetate  reagent  given  the  proper  solution  environment, 
including  the  pH  value  of  the  solution  and  correct  ratio  of  protein  to  modifier. 


Table  6-2. 

List  of  signals  observed  on  the  modified  RNase  A spectrum  shown  in  Figure  6-13, 
and  their  assignments. 


m/z 

Assignment 

1373.58 

(M  + 10  H + NHS-acetate)10+ 

1521.60 

(M  + 9 H)9+ 

1526.08 

(M  + 9H  + NHS-acetate)9+ 

1536.86 

(M  + 9H  + P + NHS-acetate)9+ 

1711.51 

(M  + 8 H)8+ 

1716.64 

(M  + 8 H + NHS-acetate)8+ 

1723.57 

(M  + 8H  + P)8+ 

1728.92 

(M  + 8 H + P + NHS-acetate)8+ 

1741.01 

(M  + 8 H + 2P  + NHS-acetate)8+ 

Methodology  Development  for  Locating  the  Protein  Modification  Site 

After  the  protein  was  proven  to  be  successfully  and  routinely  modified  with  NHS- 
acetate  reagent,  the  next  step  was  to  identify  the  exact  site  of  the  modification.  One 
question  which  needed  to  be  answered  was:  Does  the  modifying  reagent  really  modify  the 
N-terminal  lysine  as  it  was  originally  anticipated?  To  answer  this  question  is  not  as  easy  as 
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it  seems.  In  the  following  few  sections,  the  methodology  developments  of  searching  for 
the  proper  approach  will  be  described. 
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Figure  6-12. 

Comparison  of  the  spectrum  between  the  native  RNase  A sample  (see  the  top 
spectrum)  and  the  singly  modified  RNase  A with  NHS-acetate  (given  in  the 
bottom  spectrum). 


In  the  beginning,  pure  gas  phase  dissociation  techniques  such  as  sustained  off- 
resonance  irradiation  collisionally  activate  dissociation  (SORI-CAD),  photodissociation 
(PD),  and  skimmer  capillary  collisionally  activated  dissociation  (SC-CAD)  were  tried  to 
achieve  the  goal  of  partially  sequencing  the  protein  in  gas  phase.  However,  as  the  SORI- 
CAD  spectrum  in  Figure  6-8  shows,  these  techniques  did  not  generate  any  useful 
fragments  which  could  be  used  to  aid  the  partial  peptide  sequencing  in  the  gas  phase  to 
identify  the  peptide  modification  site. 
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The  next  method  checked  was  the  peptide  mapping  method  with  mass 
spectrometry.  As  described  earlier,  the  first  step  of  the  method  is  to  digest  the  RNase  A 
with  a proper  protease  (trypsin  in  this  case);  then  to  apply  the  reverse  phase  high  pressure 
liquid  chromatography  (HPLC)  technique  to  separate  the  small  peptide  mixture  resulting 
from  the  digestion  process;  and  finally  to  use  ESI-FTICR/MS  to  examine  each  HPLC 
aliquot  and  sequence  the  peptide  using  the  MS"  method.  Figure  6-13  illustrates  the 
general  guidelines  used  in  the  present  project.  The  first  step  is  to  denature  the  RNase  A 
with  guanidine  HC1,  which  cleaves  the  hydrogen  bonds  and  eliminates  the  Coulombic 
interactions  within  the  protein.  The  four  disulfide  bonds  are  more  accessible  to  the 
environment  following  the  protein  denaturation,  so  the  second  step  is  to  reduce  these 
disulfide  bonds  with  a proper  reducing  reagent  such  as  dithiothreitol  (DTT).  Since  the 
process  of  reducing  the  disulfide  bonds  opens  four  sulfur-sulfur  bonds,  it  provides  eight 
highly  reactive  -SH  groups,  and  the  immediate  following  step  is  to  protect  or  block  these 
active  -SH  groups  with  4-vinyl  pyridine  to  prevent  reformation  of  the  disulfide  bonds.  In 
the  last  step,  the  protected  protein  is  ready  to  be  digested  with  the  trypsin  enzyme,  which 
only  cleaves  the  carbon  side  of  peptide  bond  at  arginine  and  lysine  except  when  lysine  is 
located  adjacent  to  a proline  amino  acid  residue,  to  yield  a peptide  mixture.  Before  MS  is 
used  to  identify  the  peptides  in  the  mixture,  reverse  phase  HPLC  is  used  to  separate  the 
peptide  mixture. 

With  this  method,  both  native  and  modified  RNase  A are  digested  with  trypsin, 
which  should  cleave  both  proteins  at  the  similar  sites.  The  peptide  mixtures  resulting  from 
both  proteins  are  essentially  the  same  except  one  or  more  of  the  digested  peptides  from 
the  modified  protein  should  be  41  Da  heavier  than  the  corresponding  peptide  from  the 


133 


1 


Denature  w/ 
Guanidine  HC1 


I 


Disulfide  bonds 
reduction  with 
DTT 


I 


Protect  S with 
4-VP 


S s S s 

I I I I 

4-VP  4-VP  4-VP  4_VP 


1 


Trypsin  digestion 


r 


Peptide  mixture 


Figure  6-13. 

Pictorial  representation  of  the  peptide  mapping  method  including  the  processes 
such  as  protein  denaturation,  disulfide  bonds  denaturation,  and  protein  digestion. 
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native  RNase  A.  When  the  peptides  of  the  native  RNase  A are  identified,  they  then  are 
compared  with  the  digested  peptides  obtained  from  the  modified  RNase  A to  fish  out  the 
modified  fragment.  Once  the  modified  fragment  is  identified,  MS"  experiments  will  be 
performed  to  pinpoint  the  specific  site  of  the  RNase  A modification. 

Denatured  and  Reduced  Forms  of  RNase  A 

Figure  6-14  a)  and  b)  display  the  spectral  results  of  the  native  RNase  A and  its 
reduced  form.  After  examining  the  two  spectra  given  in  Figure  6-14,  it  can  be  seen  that 
the  multiple  charge  distribution  envelope  of  the  reduced  RNase  A is  shifted  to  a lower 
mass-to-charge  range  (meaning  higher  charge  states)  than  the  native  form.  This 
observation  can  be  explained  by  noting  that  the  reduced  protein  provides  more  possible 
basic  sites  for  protonation,  leading  to  higher  charge  states  following  protonation.  Such  a 
charge  envelope  shifting  phenomenon  has  also  been  observed  by  Smith’s  group.153  For 
example,  the  charge  states  observed  for  the  reduced  form  of  RNase  A range  from  10+  to 
19+  while  the  native  RNase  A only  produces  charge  states  ranging  from  T to  1 1+.  The 
reduced  RNase  A reported  in  Figure  6-14  b is  about  850  Da  heavier  than  its  native  form 
(see  Figure  6-14  a),  and  the  850  Da  is  due  to  the  addition  of  the  eight  protecting  groups 
(4-vinyl  pyridine)  to  the  eight  formerly  bridging  sulfur  atoms  of  the  reduced  RNaseA. 

Identifying  the  HPLC  Separated  Peptides  with  MS 

The  fourth  step  of  the  general  approach  of  the  peptide  mapping  method  shown  in 
Figure  6-13  is  to  digest  the  disulfide-reduced  proteins  with  trypsin.  The  peptide  mixture 
generated  from  the  digestion  process  is  then  diagnosed  by  MS  after  the  peptides  are 
separated  using  the  reverse  phase  HPLC  technique.  Before  MS  was  used  to  examine  the 
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Figure  6-14. 

The  spectrum  of  the  native  RNase  A is  given  in  a),  and  the  spectrum  of  the 
reduced  form  of  RNase  A is  shown  in  b). 
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peptide  mixture  from  the  RNase  A samples,  a HPLC  peptide  standard  mixture  purchased 
from  Sigma  was  used  to  test  the  reliability  of  HPLC  and  MS.  The  HPLC  peptide  standard 
mixture  contains  five  peptides  which  are  listed  with  their  expected  m/z  values  of  1+  and 
their  available  2+  charge  states  in  Table  6-3. 


Table  6-3. 

List  of  peptides,  with  their  expected  m/z  values  of  1+  and  the  available  2+  charge 
states,  as  contained  in  the  HPLC  standard  peptide  mixture  purchased  from  Sigma 
Chemical  Co.  This  peptide  mixture  is  used  to  test  the  reliability  of  the  reverse 
phase  HPLC  method  and  MS. 


Peptide 

m/z 

Assignment 

Gly-Try  (1) 

239.2 

(1  + H)+ 

Val-Try-Val(2) 

380.2 

(2  + H)+ 

methionine  enkephalin  acetate  (3) 

574.3 

(3  + H)+ 

leucine  enkephalin  (4) 

556.8 

(4  + H)+ 

279.2 

(4  + 2H)2+ 

angiotensin  II  acetate  (5) 

1047.8 

(5  + H)+ 

524.3 

(5  + 2H)2+ 

Figure  6-15  gives  the  complete  spectrum  of  the  peptide  mixture  with  all  five 
peptides  simultaneously  observed  plus  some  unidentified  peaks  of  low  m/z.  Peptide  1 at 
m/z  239.2  is  seen  with  very  weak  intensity  compared  to  others;  the  2+  charge  states  of  4 
and  5 are  also  seen  with  strong  intensities  along  with  their  1+ molecular  ions; 
and  the  1 + charge  states  of  peptides  2 and  3 are  also  clearly  seen.  The  results  obtained 
from  this  standard  peptide  mixture  again  illustrate  the  phenomenon  that  as  a peptide  gets 
larger,  its  charge  states  detected  by  MS  also  grow  higher  because  larger  proteins  have 
more  available  sites  for  protonation. 

The  spectra  in  Figures  6-16,  6-17,  and  6-18  correspond  to  four  HPLC  fractions 
collected  from  the  standard  peptide  mixture.  In  the  first  HPLC  fraction,  the  MS  spectrum 
a)  reveals  that  the  1 + charge  state  of  the  peptide  2 and  the  2+  charge  state  of  the  peptide  4 
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Figure  6-15. 

Mass  spectrum  of  the  HPLC  standard  peptide  mixture  purchased  from  Sigma 
Chemical  Co. 
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Figure  6-16. 

Mass  spectrum  of  the  first  HPLC  aliquot  collected  from  the  HPLC  standard 
peptide  mixture. 
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Figure  6-17. 

Mass  spectrum  of  the  second  HPLC  aliquot  collected  from  the  HPLC  standard 
peptide  mixture. 


mass-to-charge  ratio 


Figure  6-18. 

Mass  spectrum  of  the  third  HPLC  aliquot  collected  from  the  HPLC  standard 
peptide  mixture. 
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are  eluted  together.  The  similar  unidentified  low  m/z  values  seen  in  the  peptide  mixture 
are  also  seen  in  this  HPLC  fraction.  One  new  peak  at  150  is  also  seen  in  the  first  HPLC 
fraction,  and  it  is  due  to  some  contaminations  from  the  HPLC  column  since  it  is  present  in 
all  other  HPLC  fractions  but  is  absent  from  the  peptide  mixture  spectrum  given  in  Figure 
6-15.  The  spectrum  given  in  Figure  6-17  shows  that  the  majority  of  the  1+  ions  of  the 
peptide  2 are  collected  in  the  second  HPLC  fraction.  The  spectrum  given  in  Figure  6-18 
gives  mainly  the  1+  ions  of  the  peptide  3. 

Obviously,  the  results  shown  in  Figures  6-16  to  6-18  demonstrate  that  a HPLC 
separated  peptide  mixture  can  be  unmistakably  identified  using  ESI-FTICR/MS.  These 
preliminary  experiments  confirm  the  capability  of  ESI-FTICR/MS  to  detect  the  HPLC 
aliquots,  and  show  that  the  technique  is  ready  to  be  used  for  investigating  the  digested 
native  RNase  A and  its  modified  form  to  locate  the  exact  modification  site  of  the  RNase 
A. 


Conclusions  and  Future  Work 

Work  performed  on  this  project  has  demonstrated  that  ESI-FTICR/MS  can  be  an 
important  analytical  technique  for  following  chemical  modifications  of  enzymes.  The 
method  can  provide  rapid  sample  purity  screening  for  samples  commercially  bought  or 
prepared  in  the  laboratory.  For  example,  the  early  detection  of  many  contaminants  at 
different  stages  of  the  project  pointed  out  the  need  for  sample  cleaning  before  proceeding 
to  the  next  stage.  Due  to  its  high  mass  resolving  power  characteristic,  the  ESI- 
FTICR/MS  technique  has  been  demonstrated  to  successfully  and  routinely  identify  the 
singly  modified  RNase  A with  N-hydroxysuccinimide  acetate. 
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Studies  of  gas  phase  dissociation  techniques  with  SORI-CAD,  CS-CAD,  and  UV- 
visible  photodissociation  did  not  generate  enough  fragments  for  peptide  sequencing  in  the 
gas  phase.  It  seems  likely  that  the  combination  of  peptide  mapping  with  MS  detection  will 
still  be  the  best  way  to  attack  the  problem  of  pinpointing  the  exact  site  of  the  protein 
modification.  A proper  method  used  to  locate  the  exact  modification  site  is  still  being 
intensively  explored. 

Much  has  been  achieved  on  this  project  so  far;  however,  considerable  work  still 
needs  to  be  done  in  future.  An  outline  of  the  future  work  needed  to  be  done  is  shown  in 
the  flow  chart  of  Figure  6-19.  Once  the  exact  modification  site  is  accurately  identified,  a 
metal  cofactor  will  be  covalently  bound  to  the  modified  protein,  and  then  the  reaction  rate 
of  the  newly  modified  enzyme  will  be  investigated.  Of  course,  the  FTMS  technique  will 
be  one  of  the  essential  tools  for  every  step  of  future  studies. 
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Figure  6-19. 

An  outline  for  the  future  work  needed  to  be  done  on  this  project. 


CHAPTER  7 

CONCLUSIONS  AND  FUTURE  WORK 


As  can  be  seen  from  the  projects  reported  throughout  this  dissertation,  Fourier 
transform  ion  cyclotron  resonance  mass  spectrometry  (FTICR/MS)  is  one  of  the  most 
versatile  analytical  techniques  for  studying  different  chemical  systems  or  different  aspects 
of  one  particular  system.  Due  to  its  ability  to  be  interfaced  with  a variety  of  ionization 
sources,  many  classes  of  molecules  which  were  previously  not  able  to  be  introduced  into 
the  gas  phase  can  now  be  studied  with  this  technique.  At  present,  FTICR/MS  is  being 
used  to  study  organic,  inorganic,  organometallic,  and  large  biological  systems.  The 
versatility  of  the  technique  is  also  illustrated  by  the  results  which  can  be  obtained  from  one 
class  of  compounds,  such  as  determination  of  the  free  energies  of  ionization  for  one 
organometallic  complex  and  investigations  of  the  reactivity  of  another  organometallic 
complex.  The  same  technique  can  be  used  to  elucidate  structural  features  of  large  or  small 
biomolecules. 

In  chapter  3,  the  free  energies  of  ionization  of  eight  ruthenocene  organometallic 
compounds  with  the  general  formula  (ri5-C5Me5)(r|5-C5H5_nFn)Ru  (n  = 0 to  5),  including 
two  isomeric  structures,  were  measured.  The  values  determined  ranged  from  152.3  to 
170.8  kcal/mol.  Although  the  addition  of  each  successive  fluorine  ligand  to  the 
ruthenocene  complexes  did  not  dramatically  increase  the  free  energy  of  ionization  of  the 
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compound,  the  ligand  additivity  effect  is  still  seen  in  the  results  obtained.  The  studies  also 
disproved  the  expectation  that  five  fluorine  ligands  might  cause  the  “perfluoro”  effect 
within  the  ruthenocene  compounds. 

Investigations  of  the  reactivity  of  a zirconium  hydroxide  complex,  Cp2ZrOH+  (Cp 
= cyclopentadienyl),  with  selected  organic  compounds  such  as  alcohols,  esters,  ethers, 
amines,  and  amides  using  external  El  for  ion  formation  showed  that  only  the  alcohol 
family  gave  a true  ion/molecule  reaction,  producing  the  alkoxide  complex,  Cp2ZrOR+. 
When  the  results  from  external  El  studies  were  compared  to  these  obtained  using  internal 
El  for  ion  formation,  different  reaction  pathways  were  seen.  For  the  alcohol  family,  both 
approaches  for  generating  the  ions  gave  rise  to  Cp2ZrOR+  after  elimination  of  a water 
molecule  during  the  reaction.  The  relative  free  energies  of  ionization  for  these  Cp2ZrOR+ 
complexes  were  determined.  The  order  of  relative  Zr-0  bond  strength  for  the  alkoxide 
complexes  studied  was  determined  as:  Cp2ZrOH+<  Cp2ZrOCH3+  < Cp2ZrOC2H5+  < 
Cp2ZrOCH(CH3)2+  < Cp2ZrOCHCH2(CH3)2+  < Cp2ZrOC(CH3)3+.  Results  obtained  with 
the  external  El  technique  showed  that  other  families  of  organic  compounds  gave  only 
adducts  as  primary  products,  while  earlier  internal  El  experiments  had  suggested  specific 
ion/molecule  reactions. 

FTICR/MS  studies  of  biological  molecules  presented  in  this  dissertation  have 
shown  that  the  technique  can  be  used  to  study  large  biomolecules.  MS  sequencing  of 
peptides  or  proteins  in  the  gas  phase  has  been  demonstrated  as  a possible  alternative 
method  for  traditional  protein  sequencing.  Although  the  peptide  mapping  MS  approach 
can  be  time  saving  for  identifying  an  acetylated  protein  when  compared  to  traditional 
chemical  methods,  MS  alone  is  still  not  commonly  used  to  sequence  an  unknown  protein. 
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However,  a more  complete  understanding  of  the  general  correlations  between  gas  phase 
protein  fragmentation  and  the  location  of  an  amino  acid  residue  in  the  protein  will  be  very 
helpful  for  sequencing  an  unknown  protein.  Since  protein  dissociation  in  the  gas  phase 
can  be  very  complex,  small  peptides  (6  or  fewer  amino  acids)  were  studied  to  serve  as 
models  for  larger  biomolecules.  The  results  of  this  study,  reported  in  chapter  5,  yielded 
three  general  correlations:  bn  type  ions  are  observed  as  the  primary  fragmentation  series 
when  an  arginine(Arg)  residue  is  located  at  the  N-terminus  of  the  peptide,  yn"type  ions  are 
seen  as  the  principal  fragmentation  series  when  the  Arg  residue  is  located  at  the  C- 
terminus  of  the  peptide,  and  a mixture  of  mainly  bn  and  yn"  type  ions  is  observed  in  the 
dissociation  spectra  when  the  Arg  residue  is  located  in  the  middle  of  the  peptide  sequence. 
The  results  obtained  have  also  shown  that  the  precursor  ions  of  higher  charge  states 
provide  more  fragmentation  information  than  the  precursor  ions  of  lower  charge  states. 
The  IRMPD  results  obtained  for  all  peptides  studied  are  compatible  with  SORI-CAD 
results.  However,  the  results  of  this  work  suggest  that  IRMPD  is  the  better  choice  for 
peptide  sequencing  in  the  gas  phase  since  it  produces  almost  exclusively  main  peptide- 
backbone  fragments. 

The  studies  of  ribonuclease  A (RNase  A)  reported  in  chapter  6 illustrate  another 
biological  application  of  FTICR/MS.  This  work  has  demonstrated  that  FTICR/MS  can 
provide  rapid  screening  of  the  purity  of  samples.  In  addition,  successful  identification  of  a 
single  modification  of  RNase  A with  N-hydroxysuccinimide  acetate  shows  that  ESI- 
FTICR/MS  can  be  a novel  tool  for  following  chemical  modifications  of  proteins.  It  can 
quickly  measure  molecular  weights  of  unmodified  and  modified  proteins  to  determine  if 
modification  has  been  accomplished.  The  work  has  indicated  that  SORI-CAD,  skimmer 
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capillary  (SC)  dissociation,  and  photodissociation  (PD)  are  less  effective  methods  for 
sequencing  large  proteins  in  the  gas  phase.  However,  the  peptide  mapping  mass 
spectrometric  method  has  shown  promising  results.  For  future  work,  the  next  step  is  to 
examine  digested  peptide  mixtures  of  both  native  and  modified  RNase  A with  ESI- 
FTICR/MS  to  identify  the  peptide  fragment  which  contains  the  chemical  modifier, 
followed  by  MS"  experiments  to  sequence  the  modified  fragment  and  to  pinpoint  the  exact 
modification  site.  Once  the  exact  modification  site  is  accurately  located,  the  catalytic 
strength  of  the  newly  modified  enzyme  will  be  examined. 
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